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Capítulo 1 


Qual é o assunto 


Começamos por estabelecer o que são as Equações Diferenciais Ordiná¬ 
rias e por que devemos estudá-las. Apresentaremos também conceitos e 
definições úteis por todo o curso. 

Não só neste capítulo, observaremos passagens que podem parecer es¬ 
tranhas ou errôneas, por exemplo, reescrever ^ = y' como dy = y' dx. Há 
muitos outros “deslizes” que podem A devem ser notados. Recorremos ao 
bom senso e boa vontade do estudante (que aprenderá muito mais ao tentar 
justificar cada passagem) ou, onde =§p aplica, a esta observação em [FiNe]: 

O rigor que a Análise ganhava no decorrer do século XIX 
começou a pôr em dúvida èprtos métodos... Os teoremas de 
existência e unicidade êe solução surgem nessa fase. ... sua 
busca através de procéàsüsAnformais se torna justificável e pro¬ 
missora, uma vez que a “solução” assim obtida pode ser verifi¬ 
cada a posteriori. 

É com esse espírito que algumas técnicas que estudaremos são relaxadas, 
porque a corretude de suas conclusões pode ser estabelecida depois. Isso 
não significa que deveremos abandonar o rigor matemático totalmente! Os 
“deslizes” deverão sempre ser justificáveis ou formalizáveis, ficando todo o 
raciocínio “sob controle”. 

Em resumo, usaremos nossos conhecimentos de Cálculo para descobrir 
coisas e, depois, prová-las. 

A título de exemplo, visitaremos uma passagem da História da Mate¬ 
mática (veja [BoMe]) sem relação direta com EDOs: 

£ertos testes de convergência, aqueles usando integrais impróprias em 
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particular, mostram que a série 

111 1 
12 + 22 + 3 ^ + ''' + tf + • • • 

converge a algum número real, mas esses testes não dizem qual é esse limite. 
Gente graúda como Leibniz e Jacques Bernoulli não soube determiná-lo. 

Então Leonhard Euler encontrou um caminho: Queremos somar inver¬ 
sos de quadrados ou, afinal, somar inversos. Naquela época, eram muito 
populares as relações de Girard para polinómios de grau arbitrário e algu¬ 
mas de suas conseqüências; Euler pensou na soma dos inversos das raízes de 
um polinómio, que é igual ao oposto (“menos”) do coeficiente do termo linear 
dividido pelo coeficiente do termo constante. (Cheque isso para ax 2 +bx+c: 
aprendemos que as raízes Xi,X 2 satisfazem X\ + X 2 4o/a e X\X 2 = c/a\ 
dividindo uma relação pela outra, vem (x\ + av/Síj >'í '2 = (— b/a)/(c/d)\ 
simplificando, obtemos (a^) -1 + ( X\ ) _1 = —b/c, como enunciado.) 

No nosso caso, Euler precisava determinar um polinómio com raízes 
l 2 , 2 2 ,3 2 ,..., n 2 ,..., de modo que a soma de seus inversos (precisamente 
a série procurada) seria obtida pela relação de Girard. Claro, nenhum 
polinómio tem um número infinito de raízes, então deveremos operar com 
um “polinómio de grau infinito”. Por outro lado, funções trigonométricas 
têm raízes que se distanciam regularmente; por exemplo, as raízes do seno 


0, ±7r(it27r,..., ±nn ,... 

Não queremos 0 em nossa lista de raízes, porque não sabemos tomar seu 
inverso; então podemos considerar a função (sen 9)/O, cujas raízes são as 
mesmas excetuando-se 0. 

Euler foi um grande manipulador de séries de potências e provavelmente 
lhe era mais natural lidar com a série em vez da própria função. Neste caso, 
podemos obtê-la dividindo-se a série de sen0 pelo próprio 6: 


= 1 - 


Eis um “polinómio de grau infinito”, mas não exatamente aquele que pro¬ 
curamos. Mais interessante do que as raízes a intervalos de comprimento 7 r 
seriam os seus quadrados 

(±^) 2 ,(±2^) 2 ,...,(±mr) 2 ,..., 

mais parecidos com os números originais. Substituindo O 2 = x, temos 


„ x x 
- 3! + 5Í “ 
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cujas raízes são tais quadrados; esse é outro “polinómio de grau infinito”. A 
soma dos inversos de suas raízes é, portanto, 

11 1 /-l/3!\ - 1 

7T 2 + ( 27 t) 2 (nn) 2 + V 1 ) 6’ 

o que leva a 

111 1 _ 7T 2 

l 2+ ¥ + 3 2+ '" + rf + '"-~6' 

Nesse raciocínio, Euler aplicou regras de polinómios de grau finito a 
séries, sem qualquer justificativa. Contudo, havendo um “candidato” (7r 2 /6) 
a ser o limite da série, provar que a série converge a esse mesmo valor é 
uma tarefa mais simples. 

v 

1.1 O que são as equações diferenciais 

O conceito de equação diferencial é melhor explicado por partes: 

1: Uma equação envolve coeficientes e funções conhecidos e uma ou mais 
incógnitas cujos valores devemos determinar. Por exemplo, 2x + 3 = 5 é 
uma equação, e x = 1 é sua única solução ou raiz no conjunto IR. Outros 
exemplos são t 2 + 6 = 5t e a; 3 sen(a: 4 ) = 1, enquanto podemos construir 
inúmeras equações fáceis ou difíceis de resolver, com muitas, poucas, uma 
única ou nenhuma solução. 

2 : Diferencial denota um conceito matemático muito próximo ao de deri¬ 
vada; portanto, espera-se que em uma equação diferencial constem deriva¬ 
das. (As expressões dy e dx são exemplos de diferenciais; [Suv] introduz 
o Cálculo por meio delas.) Como derivada de constante é zero, as equa¬ 
ções diferenciais só teriam valor prático se, de algum modo, funções fossem 
derivadas e, mais propriamente, essas funções não estivessem pré-fixadas. 
Assim, trabalharemos com incógnitas que são funções. 

3: Em cursos introdutórios de Cálculo, já trabalhamos com problemas as¬ 
sim, quando procuramos responder à pergunta: “Qual é a função cuja deri¬ 
vada é f(x) ?” Desenvolvem-se mesmo diversas técnicas de primitivização 
para solucionar esse problema. 

Outra questão familiar é, em Física, determinar leis de movimento de 
corpos submetidos a configurações variadas de forças, e tais leis são funções, 
geralmente de uma variável real (tempo) a variáveis reais (coordenadas no 
espaço). 
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4: Por exemplo, y' — 2y é uma equação diferencial. Uma solução possível é 
y[x) = e 2x : de fato, y'{x) = e 2x • 2 = 2 y(x), ou seja, satisfizemos a equação 
ao substituir y por e 2x . 

5: As equações diferenciais ordinárias, conhecidas por EDOs, são aquelas 
em que não ocorrem derivadas parciais, ou seja, buscam-se funções de ape¬ 
nas uma variável como raízes das equações. Muitas vezes, porém, operar 
com derivadas parciais é parte central da resolução de uma equação. Quanto 
ao co-domínio, há equações ordinárias (ou ainda, sistemas de EDOs) que 
se referem a funções vetoriais em vez de funções de valores reais. 

Desse modo, as seguintes equações são diferenciais ordinárias: yy’ = 0, 
y' + y" = 3x, seny' = x 2 + 2, y" = \n\y + x\, onde lie, buscam soluções 
da forma y{x), e o sistema Y = [f _g]U, em que a ingognita tem a forma 
vetorial Y(t) = [£$]. 

6 Notações para derivadas: Dada uma função suficientemente diferen- 
ciável /, escrevemos /', /", /"' para indicar as funções que são sua primeira, 
segunda e terceira derivadas. Para uma derivada de ordem superior, diga¬ 
mos n, escrevemos /(”); quando a derivada é feita em relação à variável x, 
temos /(”) = (Note que /é a própria /.) Freqüentemente, quando 
a variável t é usada como tempo, usamos pontos em vez de linhas: y, jj; 
cuidado para enxergá-los! Enfim, é comum encontrarmos equações diferen¬ 
ciais na forma adx + bdy = 0, que significa a + by' = 0 após a divisão por 
dx. 

Alguns textos antigos usam y vnl em vez de y ®. Há autores que adotam 
a convenção de indicar jpf por y x ou y' x (sendo preciso tomar cuidado com 
o que se lê), especialmente quando são muitas as variáveis envolvidas. 

Ainda com respeito a notação, explicamos que exp(;r) e e x são a mesma 
coisa, a função exponencial com base e = 2,71828... A notação “exp” é 
vantajosa quando o argumento envolve frações e integrais, o que acontecerá 
bastante em nosso estudo. 

Determine se as funções indicadas são soluções das equações dadas: 

7: y' + ytjfx = 0, y = cosa;. 

Solução: A função dada é y = cos x e sua derivada éy r = — sen x. Para 
tal %, portanto, temos 

y' + y tg x = — sen x + cos x tg x = 0, 

ou seja, a função proposta satisfaz a equação dada, ou seja, é uma solução. 
8 : (y') 2 = 5x 2 y, y = x 3 . 
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Solução: Dada y = a: 3 , temos y’ = 3a; 2 . Substituindo-se, o lado es¬ 
querdo da equação fica 9a; 4 e o direito iguala 5a; 5 . Esssa expressões não são 
iguais e concluímos que tal função não é solução da equação dada. 

Detalhe: Poderíamos obter expressões diferentes em cada lado, mas 
ainda assim elas serem equivalentes e a função ser realmente solução da 
equação. Para certificar-se, mostre que 9a; 4 e 5a; 5 são iguais apenas para 
um número finito de valores de x e, portanto, não são equivalentes. 

9: y" — y = 0, y = 5e x . 12: 2xy' = 6y, y = —x 2 . 

10: y" — y = 0, y = xe x . 13: (2x + y)dx + xdy = 0, y = —x. 

11: 2xy' = 6y, y = 7x 3 . 14: x = —uj 2 x, x = > 54cos (ut + 9 q). 

Trabalharemos preferencialmente com soluções y(x.): I —> IR. onde o 
domínio I é um intervalo aberto. Contudo, considerar equações e soluções 
complexas, ou de outras formas, é válido e importante campo de estudos. 
Vamos formalizar as definições, a título de fixação? 

15 Definição: Por uma equação diferencial ordinária de ordem n, enten¬ 
demos uma relação na “forma geral” 

F(x,y,y’,...,^ ) ) =0, 

que envolve uma função F de n + 2 variáveis definida em um conjunto D, 
uma variável independente x, uma função incógnita y = y(x ) e as primeiras 
n derivadas de y. (F pode estar “disfarçada” como uma simples adição ou 
não depender de algumas derivadas!) Muitas vezes usaremos t no lugar de 
x e, por sua vez, x no lugar de yl 

16: Resolver essa equação significa encontrar um intervalo aberto I no eixo 
Ox e uma função <p: I —y IR tal que: 

(i) <p(x), tp'(lc)FxSylp( n )(x) existam para cada x G /, isto é, ip seja n vezes 
derivável em todo o seu domínio; 

(ii) (xftipfxf,pp\x ),..., ip^ n \x)) € D para cada x € /, de modo que pos¬ 
samos calcular F nessa (n + 2)-upla; 

(iii) F(x , <p{x), p'{x), ..., ip^ (x)) = 0 também para cada x € I. 

Uma tal função <p é chamada solução da equação, ou ainda “curva integral”, 
e geralmente indicada apenas y. 
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17: Equações da forma apresentada são bastante difíceis de serem estuda¬ 
das, como observaremos na Seção 1.4 (“funções implícitas”). Portanto, no 
que segue, estudaremos equações 

yW = f{x,y, y ',...,y( n -V), 

ditas na forma normal (de ordem n). Em geral, suporemos que / é contínua 
em uma região U, isto é, um subconjunto aberto e conexo de IR™" 1 " 1 . 

18: Assim, uma equação de primeira ordem na forma normal é 

V’ = f(x,y), 

onde f:U—>ReU é uma região em IR 2 . 

A mais simples de todas é quando / depende somente de x: 

y'(x) = f(x) para todo x em um certo intervalo. 

Suas soluções são as primitivas de / ! 

1.2 De onde elas vêm 

Devemos nos perguntar qual significado ou a razão de ser do conceito 
de derivada que, afinal, constitui o material das equações diferenciais. 

Para isso, reveremos dois exemplos tradicionais de problemas que podem 
ser expressos em termos de equações diferenciais. Para resolvê-los, não são 
necessárias técnicas sofisticadas e, por isso mesmo, eles figuram em textos 
de introdução ao Cálculo. '(Veremos o método utilizado na Seção 2.1, deta¬ 
lhadamente.) Nos outros pontos, preocuparemo-nos apenas em determinar 
as equações que regem os fenômenos descritos. 

19: Por definição, 


l/ w \ dy / \ 1- ^!/(ema) ,. 

y (a) = ——(a)= hm —- = hm 

dx Ax( em a ) 


ema) .. y(a + h) - y{a) 
- = lim -,-. 


Isso é um número, enquanto y' ou y'(x) é a função que associa a cada a o 
valor y'(/&). 

Mostra-se que a função 


L(x) = y(a) + y'(a) ■ (x - a) 


é a melhor aproximação linear à função y(x) em torno do ponto a. 
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20 Aspecto geométrico: ([Suv]) Considere esta propriedade: Cada ponto 
de uma certa curva no plano Oxy é o ponto médio do segmento tangente à 
curva por esse mesmo ponto e limitado pelos eixos cartesianos. 

Neste momento, pode ser difícil esquematizar essa propriedade em um 
gráfico. (Tente!) Experimentaremos conhecer essas curvas abstratamente: 

Seja (a;, y) um ponto de uma dessas curvas e suponha que a tangente à 
curva traçada por ele intersecte os eixos nos pontos A=(a,0)eB = (0, b). 
Da hipótese de (x, y) ser o ponto médio entre A e B, deduzimos que A = 
(2x, 0) e B = (0,2 y). Para que A e B estejam ambos bem determinados, 
sempre devemos ter x,y ^ 0. 

O coeficiente angular da tangente AB é ou, simplesmente, —y/x. 
Lembramos, agora: “A derivada da função em um poqítKç o coeficiente 
angular da tangente ao gráfico por aquele ponto.” Em Viosso caso, isso 
significa que y 1 = —y/x, ou seja, 

xy' + y = 0; 

essa é a equação diferencial ordinária que descreve as curvas com a propri¬ 
edade enunciada. 

21: Vamos já resolver essa equação, retornando-a à forma 

dy . '"?/ 
dx x 

Separamos as variáveis e integramos: 

dj, = _d* ^ rdy = _ j\dx ^ t = _i„ w + c 2 . 

y x J y ./ x 

Podemos chamar C 2 — Ci de uma nova constante C: 


cy\ = C, 


e então xy = ±e c . Novamente, =te c é um número real 


ou seja, \xy\ — ê e então xy = ±e . JNovarm 
arbitrário, que chamamos de K, e a solução é 

y = K/x. 


vV 


Portanto, as curvas com a propriedade descrita são hipérboles equiláte¬ 
ras cujas assintotas são os próprios eixos cartesianos. Agora, sabendo disso, 
esquematize a propriedade enunciada. 

Atenção: Note que K pode ser negativa; isso é resultado das primitivas 
calculadas e ocorrerá com freqüência em nossos estudos. 


7 


B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1° quad. 2012 



22 : Mostramos que as curvas com a propriedade enunciada são todas hi¬ 
pérboles de uma certa forma. Devemos mostrar que, reciprocamente, todas 
essas hipérboles têm aquela propriedade? Como fazê-lo? 

Nos próximos exercícios (23-25), forme as equações diferenciais que des¬ 
crevem as curvas planas em Oxy com as propriedades descritas. Basta 
obter as equações, não é preciso resolvê-las. 

23: Em cada ponto, a distância ao longo da tangente até o eixo das abscissas 
é constante 1. 


24: Em cada ponto, a distância ao longo da tangente até o eixo das abscissas 
é igual à distância do ponto à origem. 

25: Em cada ponto, a reta tangente passa pela origem. (Em outras palavras: 
a melhor aproximação linear em cada ponto sempre “vale zero no zero”.) 

26 Aspecto mecânico: Descobriu-se, laboratorialmente, que a velocidade 
de desintegração radioativa é proporcional à quantidade (em uma medida 
adequada) da substância radioativa existente. 


Suponhamos que nos instantes t e t + Aí a quantidade de material 
radioativo era Q e Q+AQ, respectivamente. O quociente é a velocidade 
média de desintegração da substância, durante o intervalo de tempo Aí, e 
seu limite ^éa velocidade instantânea de desintegração no instante í. Por 
hipótese, Q diminui quando t aumenta e, então, ^ é negativa; também por 
hipótese, essa derivada é proporcional a Q. Representando-se o coeficiente 
de proporcionalidade por ouQ> 0 — essa é a escolha de sinal comum em 
Física —, temos 

d Q 

iit = - aQ 


ou, simplesmente, Q + aQ = 0, que é a equação diferencial do decaimento. 

Note bem que foi preciso considerar a velocidade instantânea de desin¬ 
tegração, dada pela derivada, e não a velocidade média entre o começo e o 
fim do processo. 

27: A quantidade de uma substância específica enquanto se desintegra, em 
função do tempo, será solução da equação diferencial obtida. Essa “lei 
de desintegração” deverá ter como parâmetros a quantidade inicial Q 0 da 
amostra“meia-vida”, ou seja, o período durante o qual se desintegra a 


metade da quantidade inicial do material. 


Separamos e integramos, simplificando como com a equação anterior 
anterior: 

^ = - a dt =» ln Q = -at + C => Q = Ce c ~ at . 

Q 
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Partindo do valor inicial Q o para t = 0, determinamos C por substituição: 

Qo = e c ~ a0 => C = InQo =► Q = Q 0 e~ at . 

Onde entra a meia-vida, nessa conta? Indicando a duração desse período 
por t e substituindo na lei determinada, temos Q o/2 = Qoe~ ar , donde 
t = (ln2)/a; isso requer o conhecimento de a. Em geral, quem é conhecido 
é o próprio r, porque basta observar a amostra radioativa até que reste 
metade do material original. Então temos a = (ln2)/r e finalmente 

Q = Q 0 e~ t ^/ T . 

28: Um piscicultor deseja baixar a concentração de sal, inicialmente 36g/L, 
de um dique de criação de 5.000m 3 . Ele abrirá umâ^çõmporta que escoa 
100L/s do dique e, simultaneamente, abastecerá o dique com água doce (de 
concentração salina desprezível) a partir de uma fonte com a mesma vazão. 
Forme a equação diferencial que descreve a concentração x de sal no dique 
ao longo do tempo, assumindo homogeneização imediata da água. 

Indicação: O volume de água no dique permanece constante! Durante 
um intervalo de tempo Aí, se supusermos? a concentração x constante, per- 
dem-se lOOxAí gramas de sal, de modo que Ax « — lOOxAí/5.10 6 em 
gramas por litro. 

Note: Essas variações são aproximadas porque x varia durante Aí: a 
água salobra e a água doce são trocadas simultaneamente, enquanto que 
a variação calculada corresponde a primeiro escoar toda a água salobra e 
então repor o volume com água doce. Tomando Aí —> 0, a aproximação 
torna-se correta, mas a equação fica diferencial. 

29: Com a mesma técnica usada nos exemplos desta seção, resolva a equação 
obtida no exercício anterior e determine quanto tempo a operação tomará 
para a concentração baixar a um terço da original. 

Nos próximos problemas (30-31), não é preciso resolver as equações 
diferenciais, apenas obtê-las: 

30: Um ponto material move-se ao longo de uma reta sob ação exclusiva da 
atração gravitacional exercida por outro ponto material, fixado na origem. 
Determine a equação diferencial que rege seu movimento. Ela é válida para 
descrever jb movimento na origem? Por quê? 

Indicação: Sendo x a posição do ponto móvel, a segunda lei de Newton 
escreve-se F = mx e a lei gravitacional diz \F\ = GMm/x 2 . Por que esse 
módulo é importante? Como podemos eliminá-lo? 
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31: Um poço fino é perfurado ao longo do diâmetro de um planeta esférico e 
sólido de raio R e densidade uniforme ó, sem rotação. Determine a equação 
diferencial que rege o movimento de uma sonda ao longo do eixo Ox do 
poço, com origem no centro do planeta. O que acontece se a sonda passar 
pelo centro do planeta ou pela superfície? 

Lembrete: Se a sonda estiver no poço, sofre ação gravitacional da parte 
esférica do planeta mais próxima do centro, não da casca esférica exterior. 

32: Como em outras áreas do Cálculo e da Matemática, alguns problemas 
sugeridos (sejam geométricos ou mecânicos) aparentemente são esdrúxulos 
e sem interesse. Isso é apenas aparência! Enunciados muito surreais podem 
ser oriundos de diversas situações de engenharia e projeto, mas estar apre¬ 
sentados em forma “fechada”. Nosso objetivo é formar, resolver ou estudar 
equações diferenciais, independentemente de motivação, embora freqüente- 
mente o problema possa ser esclarecido por sua própria origem. 

1.3 O que são PVIs e TEUs 

Em algumas aplicações do estudo de Equações Diferenciais, não estamos 
preocupados em encontrar todas as soluções de uma dada equação, mas em 
determinar qual delas se ajusta aos nossos dados experimentais, como vimos 
nos problemas 27 e 29. 

Consideremos esse assunto mais detalhadamente, começando por recor¬ 
dar este problema de Física colegial: 

33: Determine a lei do Movimento Retilínio Uniformemente Variado de 
um ponto material com aceleração escalar 7 constante, conhecendo-se sua 
posição e sua velocidade iniciais. 

Solução: Utilizemos as funções s(t),V(f) para designar a posição e a 
velocidade no tempo. Como a velocidade é a derivada da posição, s — V, 
desejamos escrever s = J V(t) dt. 

Mas a integral indefinida é uma família de funções, enquanto s é uma 
única função. Um modo correto de expressar aquela idéia é através do ajuste 
da constante de integração, fixando-se uma primitiva qualquer. Porém, qual 
primitiva escolher? (Veja que \/{x —V) e x/(x— 1) são diferentes, mas têm 
a mesma derivada.) Usamos, para isso, uma integral definida: 

V(u)du + A 

para alguma constante A G IR. A escolha do limite inferior da integral é 
arbitrária e feita para facilitar os cálculos seguintes; o limite superior é a 
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variável independente estudada; a variável de integração deve ser uma nova 
letra “muda”. Lembre que, de fato, essa integral definida é uma função de 
t com derivada V e deverá diferir de s, que tem a mesma derivada, apenas 
por uma constante (o valor A). 

Do mesmo modo, a aceleração é a derivada da velocidade: V = 7 , de 
modo que 


V(t) = / 7 du + B = yt + B 


o 


com B € IR. Então também 

U f 2 

s(t) = / {l u + B) du + A = 7 — + Bt + A. 

J 0 2 

Resta determinar A e B. Costuma-se, porém, fixar que s(0) = so e 
V(0) = Vo> para facilitar a inclusão de valores experimentais ou a interpre¬ 
tação da fórmula. Das expressões acima para s e 1%. substituindo-se t = 0, 
obtemos A = sq e B = Vq e podemos reescrever 


s (t) — 7^ + V 0 t + sq. 

34: Determine A e B quando s(í 0 ) = Sq e Wío) = Ui para um instante í 0 
arbitrário. 

O problema que pede s(t) de modo que s = 7 , s(to) = so, è(to) = Vo é 
um exemplo de PVI: 

35: Em geral, para uma equação de ordem n, é preciso informar valores para 
as derivadas de ordens 0 a n — 1 da solução, de modo a fixar uma função 
específica (em vez de uma família de funções) como solução da equação. 

Quando esses n valores são informados no mesmo ponto, isto é, sob o 
mesmo valor da variável independente, então se trata de um problema de 
valor inicial, ou simplesmente “PVI”. Nos termos da definição em 15-16, o 
PVI é dado por um vetor ao € IR” e um ponto inicial xg G I; a solução ip 
deve satisfazer também (<p( xo),ip'(xo ),..., <^ n_ 1 )(£o)) = «o- 

(Quando se especifica, para uma equação de ordem n, os valores da solu¬ 
ção em n pontos distintos, trata-se de um problema de valores de contorno, 
ou de fronteira. Isso é usual para equações de segunda ordem com valores 
especificados nas duas extremidades de um intervalo de definição.) 


36: ([Gui]) Determine as funções s,c tais que [*7 

. r 0 li ro 


seja solução do PVI 




Y e V(0) = 
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Solução: Multiplicando as matrizes do sistema, obtemos as equações 
s'(x) = c{x) e c'(x) = —s(x), assim como as condições iniciais s( 0 ) = 0 e 
c(0) = 1. Em algumas situações, fazer isso explicitamente pode não ajudar 
em nada, mas aqui não temos outras opções. 

Para começar, observamos que 

(s 2 + c 2 )' = 2ss' + 2 cd = 2 sc + 2 c(—s) = 0 , 

de modo que s 2 + c 2 é constante. Para determinar o valor dessa constante, 
calculamos a nova função onde temos informação: s( 0) 2 -vc(0) 2 = 1 , de 
modo que s 2 + c 2 = 1 . De onde veio essa função s 2 + e 2 , ou melhor, 
como determinar uma função semelhante para um novo problema em mãos, 
somente muita prática e inúmeras tentativas frustradas poderão explicar. 

Agora, tentemos determinar essas funções desconhecidas comparando-as 
com outras que tenham as mesmas propriedades. Aqui, pensa-se em seno e 
cosseno por causa de s 2 +c 2 = 1 e as derivadas de*y, c serem o que são! Mas 
como mostraremos que devem necessariamente ser seno e cosseno? Tome 

h( x) = (s(x) — sena ;) 2 + (c(x) — cosa:) 2 ; 

obtemos 

h!{x) = 2 (s(a;) — sena:)(c(a:) — cosa) -f 2(c(x) — cosa;)(—s(a;) + sena:) = 0 . 

Então h é constante e calculamos h(x) = h( 0) = 0. Por h ser real e ser 
uma soma de quadradosfhada termo deverá ser zero: s(a:) — sena; = 0 e 
c(a;) — cosa: = 0 , então s(a;) = sena; e c(x) = cosa:. 

37: Identifique ou formule um PVI para cada um dos problemas 26-31. 
Quais são as interpretações das condições iniciais? 

38: Dê o significado de um PVI que fosse formulado para os problemas 
estudados em 20-25. A respeito de 20, fixe condições literais e resolva-o de 
acordo com a solução apresentada para a equação diferencial. 

39: Uma cultura de bactérias continha inicialmente um milhão de germes. 
Após uma hora, a população chegou a oito milhões. Se a taxa de cres¬ 
cimento da cultura for proporcional ao número de germes, determine em 
quanto tempo a população já havia dobrado. 

Sugestão: Releia 26 e 27. 

40: Via de regra, a solução geral de uma equação diferencial ordinária tem 
tantas constantes arbitrárias quanto é a ordem da equação. Derivando-se 
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essa solução geral o quanto for necessário e substituindo-se os dados do PVI, 
obtém-se um sistema algébrico que determina os valores dessas constantes 
univocamente. 

Supondo já obtida a solução geral, como indicada, obtenha os valores dos 
parâmetros a partir das condições iniciais e encontre a solução específica: 

41: y = A cosa:, y(n) = 3. 42: y = A sen a;, y( tt) = 4. 

43: y = (C — x 2 )/2x, y(x 0 ) = y 0 . 

44: y = Asen(a: - 9 0 ), y(n) = 1, y'( tt) = 0. 

45: y = Asen(a; - 6» 0 ), y{ 0) = 0, y'{ 0) = 2. 

46: y = Ae 2x + Bxe~ x + cosa:, y(0) = 0, y'( 0) = 1. 

Este é um bom momento para esclarecer a. que serão os TEUs de que 
tanto falaremos: 

47: Um Teorema de Existência e Unicidade, ou simplesmente “TEU”, enun¬ 
cia sob quais hipóteses um PVI tem solução (daí o termo “existência”) e 
essa solução é única. (Para uma equação algébrica ax 2 + bx + c = 0, um 
enunciado análogo é este: “Se b 2 = 4pd, então existe uma única solução 
real”) 

Nesse caso, a variável dependente regida pela equação diferencial tem 
seu comportamento totalmente determinado pelas condições iniciais. Cos¬ 
tuma-se dizer que “o futuro é totalmente previsível e o passado é totalmente 
conhecido”, especialmente quando a variável independente mede tempo. 
Provaremos tais teoremas em várias situações. 

48 Curvas não se cruzam: Suponha que valha um TEU para uma equa¬ 
ção y' = f(x,y) com certa /: IR 2 —> IR. Chame os gráficos das soluções 
dessa equação de “curvas”. Mostre e lembre sempre que: 

(i) por cada ponto do plano, passa alguma curva; 

(ii) duas curvas jamais se cruzam, juntam ou separam; portanto, 

(iii) as curvas são paralelas e formam uma “foliação” (partição) do plano. 

Esquematize isso em um diagrama. Adapte essa descrição a equações de 
graus mais altos. 
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1.4 Coisas que já sabemos fazer 

Conhecendo o Cálculo básico e usando a manipulação algébrica escolar, 
já podemos resolver ou estudar algumas equações diferenciais e, também, 
aprender raciocínios e detalhes que serão úteis posteriormente. Esta seção 
explora essas aplicações e técnicas. 

Na Seção ??, apresentaremos engenhos que “resolvem” equações diferen¬ 
ciais ordinárias. O leitor com interesse em máquinas e mecanismos já pode 
estudar essa seção e interpretar visualmente como relacionam-se os termos 
de uma equação diferencial. 


49 TFC: Dada /: ]a, 6 [ —> IR contínua e um ponto (xo, yo) £ ]a, ò[ x IR, 
existe uma única y: ]a,b[ —> IR derivável tal que y'{x) = f(x) para todo 
x £ ]a,ò[ e y(x 0 ) = yo- 

Solução: Primeiramente, note que esse enunciado é um TEU; o PVI 
em questão é y' = /, y(x 0 ) = yo- 

Pelo Teorema Fundamental do Cálctílo, 

UL$) = Vo + f(s) ds 

J Xo 

é solução tanto da equação y' = f como da condição y{x o) = yo- Para saber 
porque escrevemos tal integral assim, reveja a resolução do problema 33. 

Essa solução ao PVI é única: suponha w: ]a, b[ —> IR tal que também 
w' = f. Então (t&sr- y)' = w 1 — f = 0, ou seja, w — y = k constante real. Se 
ainda w(xyo, temos k = (w — y){x o) ■ 0 , donde w = y. 


50: O mesmo raciocínio já prova mais: Se, analogamente, z satisfaz z'(x) = 
/(i).e z(io) = zq, novamente temos: ( y — z)' = 0 , então y — z é constante, 
donde y — z = y(x o) — z(x o) = yo — Zo- Assim, se as condições iniciais yo, Zo 
estão próximas, as soluções y, z estão uniformemente próximas e, melhor 
ainda, diferem por uma constante — ou seja, as soluções para diferentes 
valores iniciais são todas translações verticais umas das outras: 
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51 TVM: Se /: IR —► IR for contínua em [a, 6 ] e derivável em ]a,b[, então 
existe c £ ]a, b[ tal que f(b) — f(a ) = f(c)-(b—a). Essé é o Teorema do Valor 
Médio. Convença-se disso por uma ilustração e estude a demonstração em 
seu texto favorito de Cálculo. 

Ele já é usado em 49, formalmente, na demonstração da unicidade. 

52: O mesmo problema de existência .£ unicidade proposto em 49, conside- 
rando-se a equação y' = ky para k £ IR. 

Solução: Verifique que e kx é uma solução de y' = ky. (Na Seção 2.1, 
veremos como encontrar essa primeira solução.) Se y é solução, 

(ye~ kx y = y' e ~ k f \- y> Y( k) = kye~ kx - kye~ kx = 0, 

isto é, ye~ kx = C para alguma constante C. Logo, y(x) = Ce kx . Como 
y(x o) = yo, vemos que Ç = y 0 e~ kx ° e então y(x) = y 0 e k( - x ~ x °\ unicamente 
determinada. Se 3xy(x) = 0, então yo = 0 e portanto y = 0. 

Gráficos: Trace os gráficos correspondentes a alguns (xo,yo), em cada 
caso k > 0, k = CU® k < 0 separadamente. Observe que, pela existência 
demonstrada, uma curva passa por cada ponto de ]a, ò[ x IR; pela unicidade, 
essas curvas não se intersectam. 

Vejamos mais alguns usos do ferramental de Cálculo: 

53: Seja /: IR 2 — > IR contínua. Prove que se a equação x = f(x, x) não tem 
soluções constantes, então ela também não tem soluções periódicas. 

Solução: Já que x = f e f é contínua, deduzimos que x deve ser 
diferenciável duas vezes e x deve ser contínua, ou seja, qualquer solução x(t) 
da equação é de classe C 2 . Como a equação não admite solução constante, 
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então Vc/(c, c) 7 ^ c, isto é, /(c, 0) 7 I 0. Porque / é contínua e /(•, 0) nunca 
se anula, então f(x(t), 0 ), como função de t, também é contínua e é sempre 
positiva ou sempre negativa. Assuma Vt f(x(t),0) > 0 ; o outro caso pode 
ser tratado analogamente. 

Suponha que x(t) é uma solução periódica, isto é, existe T > 0 tal que 
\/tx{t + T) = x(t). Em [0, T], por ser de classe C 2 , a função 'X tem um 
ponto de máximo t\\ então, x(t\) = 0. Agora, usamos a própria equação 
para calcular x(t\) > 0 . 

Pelo “teste da segunda derivada”, t\ também é ponto de mínimo de x. 
Assim, x é constante em [0, T]. Como x é periódica, concluímos que ela 
deve ser constante em todo IR, contrariando a hipótese inicial. 

54 Continuidade: Vimos em 53 que a solução da equação diferencial deve 
ser de classe C 2 e fizemos uso fundamental desse fato. Vamos, então, rever 
o raciocínio: 

(i) A equação diferencial dada é de segunda ordem, então qualquer solu¬ 
ção x(i) dessa equação, para realmente ser solução (por substituição), 
deverá ser derivável duas vezes. 


(ii) Isso requer que x seja derivável e, portanto, contínua; analogamete, x 
deve ser derivável e contínuçfev 

(iii) Como /, x, x são todas funções contínuas, então também a composição 
f(x(t), x(t)) é contínua como função de t. Mas, de acordo com a equa¬ 
ção dada, essa função é xXty. Concluímos que x também é contínua, 
ou seja, que x é de classe C 2 . 

Conclusão: Os detalhes variam de equação para equação, assim como 
os domínios das funções envolvidas, mas, geralmente, uma solução de uma 
equação diferencial ordinária deverá ser derivável em algum domínio e, as¬ 
sim, contínua nele. 

Algumas soluções de equações diferenciais são descritas por expressões 
diferentes em subconjuntos de seus domínios. Os parâmetros que surgem 
nessas descrições podem ser conhecidos através do valor inicial dado e im- 
pondo-se continuidade; determine-os a seguir: 

55: Considere 


com j/(l) = 3. 


x + A 

se x < 0, 

Bx 2 

se 0 < x < 2 . 

C/x 

se x > 2 , 
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Solução: A condição y( 1) = 3 aplica-se à expressão y(x) = Bx 2 , donde 
3 = B • l 2 e obtemos B = 3. Para que y seja contínua em x = 0, devemos 
ter 0 + A = B ■ O 2 , donde A = 0. Em x = 2, devemos ter B • 2 2 = C/2, 


donde C = 24. 


56: Considere 


V 


Í A — x 2 se x < —1, 

B(x + 2) se — 1 < x < 1, 
e Kx se x > 1, 


O 


com //(O) = 6 . 
57: Considere 


com 2 /( 1 ) = —3. 


O'") 

rsr 


ÍAcos(Tnr) se x < f, 

^ 1 Bx + B/x se x > I.-- 


Na resolução em 21, transformamos as constantes de integração a fim 
de simplificar as expressões resultantes. Mesmo no Cálculo básico, quando 
do cálculo de primitivas, operamos assim algumas vezes: é a chamada “ab¬ 
sorção de constantes”. Mostre que estas substituições são válidas: 

58: —K/3, com K £ IR, equivale à/7, com C £ IR. 

Solução: De fato, C = —K/3 e K = —3C definem funções IR —► IR 
inversas uma da outra. 

59: e K , com K £ IR, equivale a C, com C > 0. 

60: In K, com K > 0, equivale a C, com C £ IR. 

61: ; com A, B,C, D £ IR, equivale a k, com k ^ 0. E ' \ ' ? 

62: A cosí? + Bsen9, com A, B £ IR, equivale a kcos(9 — 9 q), com k ^ 0 e 
0 ^ 0o < 27 t. 

Sugestão: Teste seu raciocínio com — cos0 = cos(0 ± tt). 

Há simplificações análogas que dependem da continuidade de funções: 

63: \f/x,y)\ = K, com K ^ 0, equivale a f(x,y) = C, com C € IR, desde 
que / seja contínua. Dê um contra-exemplo com / descontínua. 

64: l/l = e K g, com K £ IR. equivale a / = Cg, com C £ IR, desde que / e 
g sejam contínuas e não se anulem. 
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Note: Isso significa que / não pode “saltar” entre Cg e —Cg, exceto 
onde / e g tiverem raízes. 

65: (/) 2 = g equivale a / = ^fg ou / = — yfg, desde que / ou í^sejam 
contínuas e não se anulem. 

66 : Uma função contínua que satisfaz uma equação algébrica quadrática 
é igual a uma das duas expressões dadas pela fórmula da Bhaskara, onde 
estas não se anularem. 

Aplicaremos estas últimas análises, ou outras similares, às soluções de 
equações diferenciais que devem ser contínuas, como discutimos acima. Elas 
também valem para o estudo de funções e soluções implícitas, que veremos 
a seguir. Por exemplo: 

67 Solução singular: Considere x(y') 2 + 2 xy' — 0. 

Essa equação é quadrática em y'\ isolando esta derivada, obtemos (por 
continuidade): 

y' = -1 ± 

que requer j; / 0 e 1 + j//a: ^ 0. Essas são duas equações homogêneas (na 
nomenclatura da Seção 2.1) que estudaremos em ?? e ??. Porém, há mais 
a ser discutido aqui: 

Em x = 0, qualquer solução deverá satisfazer y = 0, como se deduz por 
substituição na equação original. 

Quando 1 + y/x < 0, O discriminante da equação quadrática é negativo 
e não há valor real de y' possível para tais valores de x e y, de modo que 
não há função real derivável que possa satisfazer a equação diferencial dada. 
Isso prova que esse caso realmente não admite soluções. 

Se 1 + y/x = ,€^ ou seja, y = — x, então as duas equações homogêneas 
acima são a mesma equação y' = — 1 , de que a mesma função y = —x é 
solução; essa é chamada solução singular da equação original. Note que 
ela não é dada, com nenhum valor de C, pelas soluções C ± 2 \fCx (se 
x > 0 ) e -C%.2y/-Cx (se x < 0 ) das equações homogêneas, coletivamente 
chamadas solução geral da equação dada. 

Sabemos ainda que não podemos usar a fórmula de Bhaskara onde se 
anular. Aqui, isso ocorre quando y = 0 e outras informações deverão ser 
utilizadas se um PVI especificar y(x o) = 0. Note que a solução constante 
0 é dada por C = 0 na solução geral e que, para qualquer valor de C, a 
solução geral anula-se somente quando x —> 0 , permitindo a “troca” de uma 
expressão por outra. 
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68 Verificação: Na Seção 1.1, verificamos se algumas funções eram, ou não, 
soluções das equações dadas. Fazíamos uma “checagem” similar quando cal¬ 
culávamos primitivas, derivando o resultado de nossa conta e comparando-o 
com a função original. 

Esse procedimento é importante e deve ser sempre feito, não apenas por 
si, mas porque pode consistir em um passo adicional na resolução de uma 
equação. 

69 Exemplo: Os métodos que estudaremos na Seção 2.1 permitirão encon¬ 
trar a seguinte solução para a equação y' +1 = \/l + y/x , de que já falamos 
acima: 

y(x) = K 2 — 2K\fx 

para x > 0 (veremos isso em ??, juntamente com o caso x < 0 ). 
Substituindo na equação, assumindo x > 0, vem: 



•> \fx — K = | K — \fx |. 


Essa última identidade — e, portanto, o fato da função proposta ser real¬ 
mente solução da equação dada — somente pode ser verdadeiro se K ^ ^fx. 

Para tanto, ou K ^ 0 pressa condição vale para todo x > 0, ou K > 0 
e essa condição somente vale para x ^ K 2 . Assim, o domínio de validade 
da solução depende da constante envolvida, ou seja, do valor inicial em um 
PVI. 

No primeiro caso, podemos escrever K = —\fC com C ^ 0 e escrever 
y{x) = C + 2 \JCx para x > 0, o que altera significativamente a forma da 
solução. No segundo, escrevemos K = \[C com C > 0 e y[x) = C — 2 \JCx 
apenas para x ti ^ C, 


Essa conclusão chama a atenção para outro assunto: 


70 Intervalo de definição: Dado o PVI y' = y 2 , y{ 0) = yo . descubramos 
o maior domínio conexo (intervalo) que uma solução pode ter. 

Pafa?'resolver essa equação, podemos proceder como na Seção 1.2: inte¬ 
gramos $r 2 dy = dx, obtendo —y _1 = x + C, donde 


<4 


y(x) = 


X + C 
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Substituindo-se o valor inicial dado, temos yo = — 1/C ou C = — 1/yo- 
Se yo / 0, isso resolve o problema; se yo = 0, observamos que a solução 
singular constante 0 é obtida tomando-se yo — t ±0 e C —> =Foo. Ambos os 
casos são contemplados na fórmula final 


y{x) = 


1 - Uox ’ 


obtida multiplicando-se numerador e denominador por —yo- 

Para y 0 = 0, essa função está definida em todo o eixo Ox. Para y 0 / 0, 
seu maior domínio é IR — {1 /í/o } • Procuremos agora o maior intervalo I que 
possa ser usado como domínio de y. Note que, para o 'PVI fazer sentido, 
precisamos que I contenha o ponto inicial 0. Concluímos assim: 


se yo > 0 então I = ]— oo, l/yo[> 
se yo < 0 então I = ]l/yo, +oo[; 
se yo = 0 então I = IR. 




Há um exercício variante da verificação de soluções: Forme equações 
diferenciais, eliminando os parâmetros constantes, cujas soluções sejam as 
funções y(x) indicadas a seguir. 


71: y = Cyfx. 

Solução: Temos 


C _ C y/x _ y 
2y/x 2yfx y/x 2x ’ 


de onde obtemos 2xy' = y com o parâmetro C eliminado. 

72: y = Cx. 

Atenção: Derivando os dois lados, temos y' = (7; essa é uma equação 
diferencial, mas há dois inconvenientes: (i) produziu-se uma equação para 
cada valor do parâmetro C, em vez de uma única equação, enquanto C 
continua a aparecer; (ii) as funções y = Cx + D também satisfazem essa 
equação, assim como y" = 0. 


73: y = 8x + C. 

74 ax + b. 

75: y = ax 2 + bx + c. 


76: y = Cx 2 . 

77: y = Ce 3x . 

78: y = Asen(5x — 9). 
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79: x 2 + y 2 = R 2 (onde a função y é dada implicitamente, veja a seguir). 


Finalmente, abordamos o assunto de funções implícitas prometido: 

Mostre que, para cada equação diferencial abaixo, a relação indicada é 
“integral”, isto é, se a segunda igualdade define implicitamente uma função 
y(x), então essa função é solução da primeira equação. Já fizemos algo 
parecido em Cálculo básico: é a chamada derivação implícita, quando basta 
derivar os dois lados de uma igualdade. 


80: (3a: — 2 y)y' = 2x — 3 y, x 2 — 3xy + y 2 = C. 

Solução: Derivando a segunda igualdade com respeito a x, obtemos 
2a; - 3 (y + xy') + 2yy' = 0, ou seja, 2a; - 3 y = 3 xy' - 2yy' = (3a; - 2 y)y', 
como necessário. 

81: x = yy', x 2 - y 2 = C. 

82: (5a: — 3 y 2 )y' = 3a; 2 — 5 y, x 3 + y 3 = 5a :y. 


81: x = yy ', x 2 — y 2 = C. 


82: (5a: — 3 y 2 )y' = 3a; 2 — 5 y, x 3 + y 3 = 5a :y. 


83: x 2 y" + 3a :y' + y = 0, x = 2e xy . 


84: Chamamos atenção para esses problemas porque, em geral, expressar 
a solução de uma equação diferencial expliçitamente, embora desejável, é 
impraticável. Por exemplo, expressar y como função de x em sen(a;y) = 
cos(a :/y) ou em x = y 3 — e v é tarefa complexa, que não temos como estudar 
enquanto nos preocupamos com as equações diferenciais em si. 

Veja a Seção 2.3, pont.oH:32, como exemplo de soluções dadas nessa 
forma. Alguns textos, como [Dem], fornecem várias respostas em forma 
implícita. 

85: Entretanto, os Teoremas da Função Implícita e da Função Inversa per- 
mitem-nos, respectivamente, estudar o comportamento local das soluções 
dadas por essas expressões. Novamente na Seção 2.3, por exemplo, ao pro¬ 
varmos um TKl’ no ifcm 139, faremos uso desses teoremas. 

Enunciaremos versões univariáveis desses resultados tradicionalmente 
expostos em cursos finais de Cálculo (como [Rud]), mas importantes não só 
em Análise, como também em outras disciplinas, como Mecânica e Geome¬ 
tria Diferencial. 

86 Função Inversa: Sejam V. intervalo aberto e /: V. IR de classe C k . 
Se Xo G H é tal que f'(x o) / 0, então existem intervalos abertos U,V de 
modo que Xo € U, f(x o) € V, f\u: U —> V é uma bijeção e sua inversa é 
dé classe C k . 
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87 Função Implícita: Sejam fl região de IR 2 e / : íl —> IR de classe C k . Se 
(xo, 2/o) G O é tal que f(x 0 ,y 0 ) = 0e ^(x 0 ,y 0 ) ± 0, então existem região 
W Ç IR 2 e intervalo aberto U de modo que (xq, yo) G W, Xq € U e, para todo 
x G U existe único g(x) G IR tal que (a;, g(x)) GWe f(x, g(x)) = 0; a função 
g- U -> IR assim definida é de classe C k e g'(y 0 ) = -f| (x 0 ,yo 0^(zo, Jo)- 

88: Dê contra-exemplos de funções cujas derivadas não respeitem as condi¬ 
ções de cada teorema. 

89: Dê contra-exemplos para enunciados “globais” de ambos os teoremas. 


•.O 
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Algumas respostas e 
indicações 


Apresentamos as respostas para uma parcela representativa dos exercí¬ 
cios. Em geral, abreviamos as demonstrações e as soluções finais. 

9: Sim. 10: Não. 11: Sim. 12: Não. 13: Sim. 14: Sim. 23: y 2 ( 1 + 

(y')~ 2 ) = 1- 24: (y') 2 x 2 = y 2 . 25: xy' = y. 28: x = —2.10 _ 5 a;. 
29: Separe: ,J f = -2.10 ~ 5 dt. Integre: \nx = -2.10 ~ 5 t + C, donde 
x{t) = K exp(—2.10 _ 5 t) com K > 0. (Note x ^ 0.) Determine K com 
x(0) = 36 e resolva x(t) = 12. 30: Sendo s(a;) = ±1 o sinal de x, temos 
rnx = —s (x)GMm/x 2 . (Se x > 0 ou x < 0 então F < 0 ou F > 0, resp., 
porque a força gravitacional é atratora.) Simplificando: xx\x\ igual a uma 
constante estrit. negativa. ^Não vale no centro porque se assumiu x ^ 0 
para escrevê-la; a força atratora aí será infinita. 31: xx\x\ = — Gi? 3 ú47r/3 
para x > R ou x < —R e x = —xG5Ait/?> para —R < x < R. (Usa¬ 
mos x\x\ em vez de a: 2 (s(a:)) =tl .) Ao passar pelo centro ou pela superfície, 
a sonda continuará seu movimento sem perturbação. No centro, a equa¬ 
ção permanecerá válida porque a força atratora será contínua e nula. Na 
superfície, a lei do movimento (solução da equação) mudará porque a equa¬ 
ção mudará. 34: s 0 + 7 Íq /2 — V 0 t 0 e B = V 0 — 7 1 0 ; fatorando-se, 

s(t) = 7 (t—fo) 2 /2+Vo(f — ío)+So- 38: Trata-se de fixar um ponto pelo qual 
a curva passa. Sol. para 20: y = xoyo/x. 39: Quinze minutos: Seja x(t) o 
tamanho da população em função do tempo. O enunciado informa x = ax 
com cte. 0. Então x(t) = ke at com cte. k > 0 (porque x > 0). Dados 
numéricos implicam 1 = ke a0 e 8 = ke al \ determine k,a e resolva 2 = ke at 
para obter t. Note também: a população dobra três vezes em lh, então do¬ 
bra uma vez a cada |h; não há proporção direta. 41: y = —3 cosa:. 42: Não 
existe função. 43: y = (2xoyo+x%—x 2 )/2x. 44: y = sen(x— 7 r/ 2 ) = — cosa:. 
45: y = 2 sen a:. 46: y = — e 2x + 2xe~ x . 56: A = 4, B = 3, K = 21n3. 
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Capítulo 2 

Primeira ordem 


As equações de primeira ordem são aquelas que envolvem apenas a variá¬ 
vel independente, a função incógnita e sua primeira derivada. Estudaremos 
diversos tipos dessas equações cujas soluções podem ser escritas explicita¬ 
mente, começando pelas de variáveis separáveis, que já utilizamos no pri¬ 
meiro capítulo. Cada forma específica comporta uma demonstração parti¬ 
cularizada de um TEU, mas apresentaremos também um estudo qualitativo 
geral e abstrato para a forma normal das EDOs de primeira ordem. 


2.1 Separação de variáveis 

90: Uma equação de variáveis separáveis tem a forma 
(VS) A = /(a:) • g(y), y(x 0 ) = y 0 , 
onde já incluímos um PVI. 

São casos particulares as equações autônomas (que veremos em 121, 
ainda nesta seção, com / = 1) e o cálculo de primitivas (com g = 1). 

91: Para resolvê-la, basta desmembrar as diferenciais e passar todas as 
ocorrências de cada variável (x,y) para cada lado da igualdade: 


dy 

g{y) 


f{x) dx. 


Desse modo, podemos integrar ambos os lados usando técnicas de Cálculo 
básico e, então, isolar y em termos de x. A constante de integração terá 
seu valor determinado pelo PVI. 
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Esse “desmembramento” das diferenciais pode parecer misterioso ao es¬ 
tudante de Cálculo, para quem foi sempre dito que ^ deveria ser tratado 
como um único “bloco” e não como uma fração. Ao demonstrar o TEU, 
veremos como esse procedimento se justifica. 

92: Para integrarmos, é fundamental que cada lado da equação contenha 
apenas uma variável. Assim, nem toda equação tem variáveis separáveis e, 
caso tenha, isso pode não ser aparente: Procure sempre fatorar a expressão 
dada para verificar a possibilidade de separar as variáveis! 

(Nota: Não conhecemos a referência bibliográfica para o próximo exem¬ 
plo, mas preferimos adotá-lo por sua simplicidade e pela discussão lógica e 
gráfica que permite.) 

93 Exemplo: x{l + y 2 ) dx — y(l + x 2 ) dy = 0, //(•'•(>),-' yo- 
Rearranjando a equação dada, obtemos 

dy _ x 1 +y 2 
dx 1 + x 2 y 

de variáveis separáveis, mas essa operação requer que a solução y{x) não se 
anule. Suponha, então, que y{x) = 0 para algum valor de x: substituindo 
na equação original, vem a^l + O 2 )*— 0(í+x 2 )y' = 0 , ou simplesmente x = 0 . 
Assim, toda solução que se anula só poderá fazê-lo na origem e doravante, 
portanto, assumimos j/o / 0 ! e que y não se anula. 

Agora, isolamos e integramos: 

'=> \ ln(l + y 2 ) — \ ln(l + x 2 ) + C => 

1 + y 2 = k( 1 + x 2 ) onde k > 0 . 

Substituindo a condição inicial y(x o) = yo, obtemos k = ou seja, 

Análise gráfica: E fácil rearranjar a identidade obtida na forma tradi¬ 
cional da equação de hipérbole: y 2 — kx 2 = k — 1 . 

Se k > 1, podemos calcular y(x) = ±\/kx 2 + k — 1; cada sinal especifica 
uma função que descreve um ramo inteiro e está definida em todo o IR. 
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Quando 0 < k < 1, há um intervalo onde kx 2 + k — 1 é negativo, de modo 
que os ramos correspondem somente a funções definidas em semi-retas. 
Como devem ser funções, cada uma descreve apenas metade de um ramo. 
Não há solução nos vértices dos ramos, de acordo com nossa discussão 
inicial. 

Para k = 1, temos y = ±.x e podemos distinguir entre essas possibilida¬ 
des usando a condição inicial, já que assumimos yo ^ 0 : se yo = Xq então 
y(x) = x; se y 0 = —x 0 então y(x) = —x. Como essas retas passam pela 
origem, podemos tomá-las como soluções definidas em todo o IR. 

Concluímos com o seguinte esquema: 



Atenção: As soluções com k = 1 são retas, ou seja, são hipérboles 
degeneradas. Contudo, não são as assintotas das demais hipérboles que 
particionam o plano. Mostre, com as fórmulas básicas de Geometria Analí¬ 
tica, que as assintotas da hipérbole 1 + y 2 = /c(l + x 2 ) são retas passando 
pela origem ç com coeficiente angular ±\Z~k. Desse modo, as hipérboles 
não têm as mesmas assintotas e “fecham-se” progressivamente conforme 
são mais afastadas da origem. 

94: Faça o mesmo estudo algébrico e gráfico para o PVI x dx + y dy = 0, 
y(x o) = yo- 

Resolva estas equações por separação de variáveis: 
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95: y' tg x = y. 99: (xy 2 + x)dx + ( x 2 y + y) dy = 0. 

96: xdy + ydx = 0. 100: y' — e~ v + 1 = 0. 

97: xyy' = 1 — x 2 . 101: x 2 y' — y = 1 — xy'. 

98: y' = e x ~ y . 102: xy' = y 2 + 1. 

Resolva os problemas de valor inicial: 

103: xdy + ydx = 0, y{ 2) = 4. 106: x 2 dy — y 2 dx = 0, y( 1) = 1. 

104: x dy + y dx = 0, y( 0) = 0. 107: y' = 2^/y liM+yQ) = 1. 

105: x 2 y' - y = 1 - xy'. l/í 1 ) = 3. 108: (:r 2 + l)y' = \J\ - y 2 , y( 1) = 1. 

Discutiremos a existência e a unicidade de soluções de (VS) em dois 
casos: (i) assumindo que g não se anula e, depois, (ii) assumindo que as 
raízes de g são isoladas e trabalhando nelas. 

109 TEU: Sejam /: ]o, b[ —> IR contínua e y: ]c, d[ —> IR contínua que não 
se anula. Então, para cada ponto (xo,yo) do retângulo R = ]a, b[ x ]c, d[, 
(VS) tem uma única solução y: / —» ]p, d[, onde / é um intervalo aberto 
contido em ]a, b[ e contendo Xq. 



Observe, em 93, os gráficos e os domínios das soluções com 0 < k < 1: 
elas não estão definidas em todo IR e são um exemplo de que, aqui, não 
podemos pretender que / seja todo o intervalo ]a, b[. 

A demonstração desse TEU será feita em dois pontos: 

110 Existência: Sejam F,G primitivas de /, i, respectivamente; isso é 
garantido pela continuidade destas duas funções. A função G é estritamente 
monótona pois sua derivada 1 não se anula e então, por continuidade, tem 
sempre o mesmo sinal. Assim, sua imagem é um intervalo aberto J e existe 
a função inversa G _1 : J —> ]c, d[. 
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Como derivada de constante é zero, podemos transladar F verticamente 
de modo que F(xq) = G(yo). Porque F é contínua, existe então um inter¬ 
valo aberto I tal que Xç, £ I Ç ]a, b[ e ainda F[I] Ç J. 

Assim, podemos definir y: I —> ]c,d[ simplesmente como a composição 
y = G ^ 1 o F. Note que obtemos y(x o) = yo- 

Para verificar a solução, derivamos a identidade F = G o y utilizando a 
Regra da Cadeia: 

/ = F' = (Goy)' = (C oy) ■ y' = — • y', 

g°y 

donde yf = / • {g o y), como desejado. 

111 Unicidade: Continuemos trabalhando com as mesmas primitivas F. G 
acima, mas agora suponhamos que temos uma solução y qualquer. Mostra¬ 
remos que, em I, essa função y é igual à mesma solução G _1 o F que já 
calculamos. 

Nesse intervalo, podemos calcular (Goy) —\F e derivá-la. Vem 

(GoyY-F' = — V —/ = 0 
g^y 

porque y é solução de (VS). Desse modo, (Goy) — F é constante e seu valor 
é, por substituição, G(y(x o)) — RÁXq) = 0. Concluímos que G o y = F e 
basta calcular G -1 em ambos os lados dessa igualdade. 

Destacamos o seguinte eScólio, que explica a abundância de soluções 
dadas implicitamente: 

112 : Se G' = l/y, obtemos a solução implícita G(y) = f(t) dt + G(yo). 

113: Agora, consideramos o caso em que g tem raízes e, portanto, não se 
pode calcular l/y; assumiremos ainda que essas raízes são todas isoladas. 
(A maioria das questões práticas satisfaz essa condição.) 

Pode-se aplicar o TEU anterior nestas restrições do domínio de y: 

(i) em um intervalo aberto entre quaisquer duas raízes consecutivas de y; 

(ii) se y(y 0 ) ^ 0 , em um intervalo aberto contendo yo onde y não se anule. 

114: Se y(y 0 ) = 0, então a função constante yo é solução do PVI (VS); 
note qu&.sua derivada é zero em Xq (e sempre). Qualquer outra solução y 
também terá y'(x o) = 0 em vista da equação: 

y'(x o) = f(xo) ■ g{y(x o)) = f(x 0 ) • y(y 0 ) = f(x 0 ) -0 = 0 . 
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Assim, como veremos nos exemplos-exercício abaixo, nas raízes de g pode 
ocorrer (ou não) o seguinte fenômeno: há soluções cujos gráficos tangenciam 
a reta horizontal y = y 0 e, desse modo, são tangentes também entre si. (Não 
há soluções que se cruzam.) 

Já que vale o TEU entre as raízes de g, são únicas as soluções entre as 
retas horizontais definidas por essas raízes, e são obtidas pelo método que 
discutimos acima. Assim, caso alguma solução conflua com, ou bifurque de, 
uma das soluções constantes, deverá assumir a forma dada pelo método. 

Conclusão: Para g com raízes isoladas, (VS) sempre tem ao menos 
uma solução e as únicas soluções possíveis da equação dada são 

(i) as funções constantes que são raízes de g (são as soluções singulares, 
também chamadas equilíbrios da equação) e 

(ii) aquelas dadas pelo método de separação de variáveis (são as regulares), 

embora (VS), como um PVI específico, possa ter não apenas uma solução. 

Conhecendo-se as soluções regulares, devemos inspecioná-las para veri¬ 
ficar se incluem, em sua forma, as singulares, mas isso é diferente da pos¬ 
sibilidade de aquelas tangenciarem estas. Se isto não acontecer, conclui-se 
por um TEU irrestrito. 

Resolva cada uma das equações abaixo com a condição y(x o) = yo e 
identifique qual das situações descritas acima ocorre. 

115: yf = —xy- 1 . 118: y' = 1 + y 2 . 

116: yf = 2xy 2 . 119: y' = 3y 2 / 3 . 

117: yf = x~ 1 y. 120: y' = x~ 2 . 

No restante desta seção, concentraremos nossa atenção em um caso par¬ 
ticular muito útil da forma (VS) e em uma modalidade típica de aplicações 
desse tipo de equações. 

121: Uma equação autônoma é aquela em que a variável independente não 
consta explicitamente e tem a forma y' = g{y). Note que essa equação tem 
variáveis separáveis! 

122: Suponha que g: I — > IR é contínua e não se anula, sendo / um intervalo 
aberto. De acordo com o TEU para (VS), fixado (xo,yo) € IR x I, o PVI 
y' = g(y), y(x o) = yo tem uma única solução definida em um intervalo 
aberto contendo a^. 

123: Mostre que podemos assumir xq = 0. 
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124: As equações autônomas surgem especialmente na descrição de fenô¬ 
menos naturais que descrevem a evolução de alguma quantidade ao longo 
do tempo, mas em que a marcação do tempo, em si, é irrelevante. 

Por exemplo, nos pontos 26 e 27, conhecemos a equação de decaimento 
radioativo 

Q = -—<? 

T 

em termos da meia-vida r da substância em estudo. Evidentemente, tanto 
a condição de valor inicial Q( 0) = Qo como a solução Q = Qoe~ t ^ n2 ' , ^ T 
envolvem a variável t que mede o tempo, mas a equação diferencial indica 
a taxa de variação de Q apenas em termos de condições contemporâneas 
(no caso, o próprio valor de Q). Em outras palavras, o fenômeno estudado 
e a equação obtida não dependem de um “calendário” específico. 

O ponto anterior sugere que, nesse caso, não há perda de generalidade 
em supor-se sempre o instante inicial sendo 0. 

Como exemplos, na próxima seção, daremos destaque a modelos popula¬ 
cionais e epidemiológicos, mas voltaremos a encontrar equações autônomas 
em outras aplicações. 

125: A título explanatório, digamos que g está definida em toda a reta 
real e tem três raízes reais yi,r/ 2 , 2 / 3 - Sabemos que as funções constantes 
correspondentes são soluções de yf = g(y) e já as chamamos de equilíbrios 
dessa equação. Outras soluções podem ou não tangenciar com, confluir 
com, ou birfucar-se de cada equilíbrio. 

126: Em caso negativo, quando nenhuma solução corta as retas y = jq, 
V = 2/2 e y = j/ 3 , a situação é descrita pelo diagrama abaixo. Ele assume 
que as soluções estão definidas e são contínuas em toda a reta real e ilustra 
estes fatos: 

• já que yf = g(y ), a solução é crescente ou decrescente conforme g é 
positiva ou negativa', 

• já que y" = g , (y) ■ y', os pontos de máximo e mínimo de g indicam 
alguns dos valores de inflexão das soluções. 
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127: Essa análise sugere a seguinte classificação dos equilíbrios: 

• estável se g > 0 antes (à esquerda, ou abaixo) do equilíbrio e g < 0 
depois (à direita, ou acima) dele; 

• instável se g < 0 antes eg >0 depois; 

• semiestável em outros casos. 

Ou seja, o equilíbrio é estável (i / 2 no exemplo) ou instável Mh, 2/3 no exem¬ 
plo) quando as soluções em seu entorno aproximam-se ou afastam-se, res¬ 
pectivamente. A terceira possibilidade não é contemplada no exemplo e 
ocorre quando as soluções em um lado do equilíbrio aproximam-se dele, 
enquanto as do outro lado afastam-se. 

Identifique e classifique os equilíbrios destas equações: 

128: y' = 2y{l - y ). lSOix^sWa: - 2). 

129: y' = —3y(y — l)(y + 2) 2 . 131: * = sen( 2 a;) cos x. 

132: Mostre que todas as soluções entre os mesmos dois equilíbrios conse¬ 
cutivos (ou abaixo do menor equilíbrio, ou acima do maior equilíbrio) são 
todas translações horizontais umas das outras. 

2.2 Dinâmica populacional 

Os modelos de dinâmica populacional são excelentes exemplos de equa¬ 
ções autônomas justamente por não dependerem de calendários. Aqui, estu¬ 
daremos a formulação de alguns desses modelos e ilustraremos o estudo de 
equações autônomas através desses casos específicos. Qual modelo usar ou 
montar depende, afinal, de diversas hipóteses particulares de cada estudo e 
cada ciência. 

133: Embora populações sejam intuitivamente medidas por variáveis dis¬ 
cretas (cujos valores são números naturais), elas podem ser representadas 
por variáveis contínuas x(t) que, na perspectiva deste assunto, são mais 
suscetíveis às ferramentas do Cálculo. Partes fracionárias, especialmente 
de números grandes, podem ser simplesmente ignoradas. 

Outras quantidades frequentemente representadas assim são: propor¬ 
ções de indivíduos com certas características, dentre toda a população; 
probabilidades de eventos; concentrações ou densidades populacionais; bio- 
massa (soma das massas, ou outra medida qualquer, de todos os indivíduos 
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da população); totais de capital e montante; até massas d’água ou tamanhos 
de ambientes ou corpos isolados. 

Trabalharemos sempre com a condição inicial a:(0) = Xq', para outros 
valores do instante inicial, procede-se por translação horizontal da solução, 
como vimos na seção anterior, ou adota-se uma convenção manual como 
“1990 é o ano base”, de modo que 2015 corresponde a t = 25. 

134 Malthus: A velocidade de crescimento proporcional à própria popula¬ 
ção é o postulado mais simples desse estudo. Ele corresponde à equação 

x = ax 

com uma constante de proporção a > 0. A solução 

x(t) = Xoe at 

descreve o famoso “modelo de crescimento exponencial”. 

A mesma expressão é válida quando a hipótese de interesse em 

casos como o decaimento radioativo ou a autofagia. 

Experimentaremos algumas variações desse postulado e do modelo: 

135: Se, em uma população de x indivíduos, há sempre x/2 machos e x/2 
fêmeas, podemos supor que a taxa de crescimento x seja proporcional ao 
número de encontros entre machos é; fêmeas que resultam em cópula. Esse 
número, por sua vez, é proporcional aos números de machos e fêmeas, então 
x deve ser proporcional a (x/2)(.r/2). ou simplesmente 

x = ax 2 

para a > 0. (Resolva essa equação.) 

136: Um modelo cuidados incorporará não somente o aumento populacional 
por nascimentos, mas também as perdas por falecimentos e até os fluxos 
migratórios. Diversas expressões podem ser usadas para representá-los, de 
acordo com as hipóteses feitas e a ciência estudada. 

Por exemplo, se a taxa de mortalidade (proporcionalmente à população) 
for constante, então o termo —bx (onde b > 0 ) é somado ao termo de 
crescimentá', assim: os dois modelos anteriores ficam 

x = ax — bx = (a — b)x, 
x = ax 2 — bx. 
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Outro caso são as situações de grande violência, como guerras civis ou 
criminalidade. O mesmo raciocínio do ponto anterior pode ser usado para 
justificar o termo —bx 2 para as perdas por morte, proporcionalmente aos 
encontros entre dois indivíduos da população. Obtemos 

x = ax — bx 2 , 
x = ax 2 — bx 2 = (a — b)x 2 . 

Resolva essas quatro equações e identifique se há mudanças estruturais 
em suas soluções, com relação às anteriores. Estude ainda a conseqüência 
das hipóteses a > b e a < b. 

137 Verhulst: Populações muito grandes enfrentam escassez de alimentos 
e espaço e têm seu crescimento dificultado e freado. Um modo de incorporar 
esse fato ao modelo malthusiano é impor um redutor M — x que diminuirá 
a taxa de crescimento quando x é (menor que, mas) muito próximo de um 
limitante M: 

x = a(M — x)x. 

Essa é a chamada “equação do crescimento logístico” ou, simplesmente, 
equação logística. 

Ela é aplicada com sucesso em diversos estudos demográficos e também 
na modelagem de epidemias, em, : que x mede a proporção de infectados e a 
disseminação cresce proporcionalmente ao número ( 100 %— x)x de encontros 
entre não-infectados e já infectados. Pode ainda ser interpretada segundo 
nossa experimentação anterior, quando posta na forma 

x = aMx — ax 2 = ax — (3x 2 

(tomando-se a = /3eM = a//3) que representa crescimento malthusiano 
com mortalidade por violência. 

138: Resolva a equação de Verhulst por separação de variáveis, obtendo 


1V1 

1 — Cexp(— Mat)' 

Use a condição inicial x(0) = Xq para determinar C e reformule x como 




x(t) = 


xo + (M — xo) exp (—Mat)' 


Essa forma da solução inclui ambos os equilíbrios 0 e M; mostre que, para 
outros valores de xq, ela nunca iguala nem 0 nem M , provando “existência 
e unicidade” sempre. 


B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1- quad. 2012 


139: Podemos estudar o comportamento dessa função e esboçar seu grá¬ 
fico utilizando as técnicas de Cálculo elementar. Por exemplo, assumindo 
a,M > 0 e 0 < Xo < M, mostre que x tem domínio IR e sempre 0 < x < M 
e calcule 

lim x{t) = 0, 

t—t — O O 

lim x(t) = M. 

í-H- oo 

Entretanto, como veremos no próximo ponto, não é preciso manipular 
a fórmula explícita e complicada de x para efetuar a análise das primeira e 
segunda derivadas. 

140: Já conhecemos a primeira derivada x, dada pela própria equação: 

x = a(M — x)x > 0 => x é estritamenfe^sêscente 
(para a,M >0). De modo análogo, calculamos 

* = ít ±= íM M ~ x ) x \ = 

= a (0 — x)x + a(M -^ç)x = 

= a[—a(M — x)x] + a(M — x)[a(M — x)x\ = 

= a 2 (M — x)(M — 2xfx. 

(Os colchetes indicam as substituições de x por a(M — x)x, para a derivação 
e para eliminar as duas ocorrências após a derivação.) Então temos x = 0 
quando x = M/2 (e nos casos x = 0 ,M que descartamos), x > 0 quando 
0<x<M/2ex<0 quando M/2 < x < M. 

Observe que não determinamos os pontos t onde as primeira e segunda 
derivadas anulam-se, mas os valores x nesses pontos. Assim, contrariamente 
ao hábito dos cursos de Cálculo, marcamos esses valores de interesse no eixo 
das ordenadas, não no das abscissas. 

141: O valor M/2 pode ser interpretado assim, no contexto de contágio 
epidêmico mencionado na introdução à equação: ao atingir mais de 50% da 
população suscetível, o contágio deve desacelerar. 

Inicialmente, há poucos portadores da doença e a disseminação é lenta, 
embora se acelere com o aumento dos portadores, porque há mais e mais 
portadores contaminando indivíduos sãos. Quando 50% da população está 
contaminada, a disseminação tem velocidade máxima, mas deve iniciar sua 
desaceleração porque, com a diminuição da parcela sã, encontros entre por¬ 
tadores e não-portadores tornam-se progressivamente menos frequentes. 
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142: A análise acima fez uso fundamental da hipótese 0 < Xq < M , que im¬ 
plicou que 0 < x(t) < M para todo instante t. Ela foi o que, em particular, 
garantiu valor sempre positivo para o denominador de x(t) e a continuidade 
dessa solução em toda a reta real. 

Em geral, devemos considerar quando o denominador anula-se: 


(M — xo) exp (—Mat) = 0. 


Resolvendo em t, obtemos 

T _ -1 ln £o 
M a Xq — M ’ 

chamado ponto de turbulência. Sua existência requh? (xq — M) > 0 e, 
por sua vez, isso equivale às opções xq > M ou Xo <vlHporque assumimos 
M> 0). 

143: O gráfico da solução x(t), quando há turbulência, tem dois ramos, à 
maneira das hipérboles: um para t < T e outro para t > T; um está abaixo 
do equilíbrio 0 e outro acima do equilíbrio Ü/fr ■ (A correspondência entre as 
duas dicotomias depende do sinal de Xo-) 

A solução somente está continuamente definida no intervalo que contiver 
o instante inicial 0, seja ]— oo,T[ ou }T, ÍK)[. Distinguem-se, portanto, duas 
situações: 

(a) Xo > M, donde xq/{xo — M) ?1 e T < 0; a solução está definida con¬ 
tinuamente em ]T, oo [ e tem o gráfico superior no diagrama a seguir; 
mostre que 

lim x(t) = oo, 
t->T+ 

lim x(t) = M. 

í—>•+ OO 


(b) x 0 < M, donde x 0 /{x 0 — M ) < 1 e T > 0; a solução está definida 
continuamente em ]— oo,T[ e tem o gráfico inferior no diagrama a 


seguir; mostre que 

/-CV 




lim x(t) = 0 , 

t — y —oo 

lim x(t) = —oo. 

t — >T~ 


144: A análise restante é similar à anterior, notando-se que x < 0 e calcu- 
lando-se o sinal de x. Concluímos que o equilíbrio 0 é instável e o equilíbrio 
M é estável. 
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145: O que muda se tomarmos a < 0 e/ou M <0? Mostre que, então, o 
equilíbrio M identifica um “limiar” a partir do qual uma população pode 
prosperar ou extinguir-se. 

Há duas variações do modelo de Verhulst que incorporam caça preda¬ 
tória ou extração sustentável de recursos renováveis. Lembramos que essa 
equação já tem uma taxa de mortalidade intrínseca (/3x 2 , acima), então 
identificamos explicitamente um termo adicional: 

146: Se a extração ocorrer a um número fixo de indivíduos por unidade de 
tempo, vem 

x = a(M — x)x — b 
para alguma constante b > 0. 

Determine e classifique os equilíbrios dessa equação. Mostre que um dos 
equilíbrios está entre 0 e M e, para esse valor x, determine b em termos dos 
coeficientes intrínsecos a,M. O que é necessário para que a extração seja 
sustentável a essa taxa? 

147 Schaefer: Se a extração ocorrer com facilidade proporcional ao tama¬ 
nho da população, ou seja, “quanto mais houver, mais fácil caçar”, obtemos 

x — a(M — x)x — bx 

onde novamente b > 0. 

Determine e classifique seus equilíbrios. Qual valor de b permite caça 
renovável, sem alteração da população? 

Em muitos casos, se o número de indivíduos estiver abaixo de um “limiar 
inferior” ainda positivo, então a reprodução e o crescimento são inviabiliza¬ 
dos e a população extingue-se. Vejamos: 

148: Suponha M > L >S1$ onde L é o limiar. Mostre que o modelo 

x = a(M — x) (x — L), 

com equilíbrios L, M, reduz-se à equação de Verhulst. 

149 Volterra: Se acrescentarmos o redutor x — L, que corresponde à difi¬ 
culdade de crescimento quando x é (maior que, mas) muito próximo de L, 
obtemos 

x = a(M — x)(x — L)x. 

Determine e classifique os equilíbrios dessa equação e verifique o que 
ocorre quando xo < L. Resolva-a explicitamente e analise sua solução com 
as mesmas técnicas que usamos no estudo da equação de Verhulst. 


35 

B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1° quad. 2012 



Finalizamos esta seção com diversos problemas que têm o mesmo espí¬ 
rito do estudo de dinâmica populacional: 

150: Um cometa esférico e homogêneo é vaporizado pelo calor solar. Ele 
perde massa proporcionalmente a sua exposição à radiação solar, ou seja, 
seu volume decresce em uma velocidade proporcional à área dé^sua super¬ 
fície. Escreva e resolva a equação diferencial que rege o raio do cometa em 
função do tempo, com o PVI que associa Rq ao raio na origem dos tempos. 

151: Mostre que o mesmo modelo descreve a evaporação da água contida 
em uma taça cônica a uma temperatura constante. 

2.3 Equações lineares 

152: Uma equação de primeira ordem da forma 

V = P(x) ■ y +^4) 

é chamada linear. Quando q = 0, trata-se de uma equação linear homogênea 
que tem variáveis separáveis. 

Comecemos imediatamente por um TEU. Sua demonstração é construtí- 
vel, isto é, apresenta uma fórmula para a função procurada. Porém, não há 
necessidade de memorizá-la: aprecie o “método da variação das constantes” 
utilizado para deduzi-la e reproduza-o nos exercícios propostos. 

153 TEU: Suponha que p,q: I —> IR são contínuas sobre um intervalo 
aberto I. Dados xq £ I e yo £ IR, existe única y: I —> IR de classe C 1 tal 
que y' = py + q e //(;%.) ' Vo- 

Solução: Inicialmente, considere a parte homogênea y' = py, cuja solu¬ 
ção por separação de variáveis é 

y{x) = y 0 exp(J x xo p(s)ds). 

Aqui e abaixo, usamos o TEU para (VS) com a hipótese de p contínua; essa 
expressão inclui a solução constante 0 (pondo-se yo = 0 ) e nenhuma outra 
solução se anula. 

Neste momento, substituímos yo pelo símbolo C: agora, procuraremos 
soluções da equação completa, na forma 

y{x) = C(x) exp(f x xo p(s) ds ); 
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este é o chamado método da variação das constantes em que supomos C fun¬ 
ção de x, isto é, U C varia com x (Pospomos para o Capítulo 3, ponto 346, 
uma motivação e um contexto para essa idéia.) 

Calculamos a derivada do produto: 

y'(x) = C'(x) exp(f xo p(s) ds) + C(x) exp (f* g p(s) ds)p(x). 

Também temos, por substituição, 

p{x)y(x) + q{x) = p(x)C(x) exp(f xo p(s) ds) + g(a(£ 

Assim, conhecemos os dois lados da equação y' = py + q. Cancelando-se os 
termos com p, sobra 

q(x) = C'(x ) exp (J xo p(s) ds), 

cuja solução é 


C(x) = C{xo) + J X Xo q(r) exp(- f Xg p(s) ds) dr. 

Retornando C(xo) = C = yo, obtemos a solução procurada 

y(x) = y 0 exp(f xo p(s)ds) + f* o q(r) exp(f x p(s) ds) dr. 

Para provarmos a unicidade, suponha que tenhamos soluções y\, y -2 para 
o mesmo PVI. Observe que y = yi — z /2 satisfaz y(x o) = 0 e y' = py. Como 
no início, y{x) = y(xo) exp(f xg p(s$xls) e substituindo y(x o) = 0 informa 
V = 0, ou seja, yi = y- 2 - 

154 Exemplo: y' = 2 x~ l y — 5. 

I 2 passo: A parte homogênea é y' = 2 x~ 1 y, que pode ser resolvida por 
separação de variáveis: 


— =^ 1 ^— hi \y\ = 2 ln \x\ + C\ => y = Cx 2 . 
y x 

2 - passo: Fazemos a “variação da constante”, ou seja, tomamos a função 
y(x) = C(x) ■ x 2 . Substituindo-a em y' = 2 x~ l y — 5, obtemos 

C(x) ■ x 2 + C(x) ■ 2x = 2x~ x ■ C(x) ■ x 2 — 5. 

Note que b somando 2xC(x) aparece em ambos os lados dessa igualdade. 
Após eliminá-lo, determinamos C(x): 

C'{x) ■ x 2 = — 5 => C(x) = —5x~ 2 => C(x) = 5a ; -1 + K. 
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3- passo: Finalmente, concluímos que 

y(x) = C(x) ■ x 2 = 5x + Kx 2 . 

Verifique diretamente que y = 5x + Kx 2 é solução de y' = 2x~ x y — 5. 

(Não temos referência para o próximo exemplo.) 

155 Exemplo: y' = (1 + sen x)y + cosx, y(n) = e. 

I 2 passo: Resolva primeiro apenas a parte homogênea y' = (l+senx)y, 
de variáveis separáveis, obtendo y{x) = C exp(x — cosa:). 

2° passo: Escreva y(x) = C(x) exp(x — cosa:) e derive, substituindo na 
equação original: 


C'{x) exp(x — cosa;) + C(x) exp(x — cosx)(l + sen a:) = 

^ ^ ^ 
compare isto. . . 

p(x) / y(x) q(x) 

= (1 + sen a:) [C(x) exp(x — cosa:)] + [cosa;]. 


Cancelando, obtemos 


C'(x) exp(x — cosa:) = cosx 


e então 


r 

C(x) = / cossexp(coss — s) ds + K. 


Essa é uma integral mais difícil que não é, agora, nossa preocupação; deixe¬ 
mo-la indicada. Note que utilizamos o ponto inicial tt como limite inferior 
de integração, o que facilitará o último cálculo da resolução. 

3 2 passo: Substitua C(x): 

y{x) = Q{x) exp(x — cosa;) = 

[X 

= Kexp(x — cosa;) + exp(x — cosa;) / cos sexp(cos s — s) ds. 

J 7T 

Agora y(n) = e implica e = Ke K+1 + e 7r+1 0, donde K = e - '"’ e obtemos 

, V) f x 

y[x) = exp(x — 7 r — cos x) + exp(x — cos x) / cos s exp(cos s — s) ds. 

J 7T 

Resolva estas equações lineares: 
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156: y' = y tg x + 4 sen x. 
157: y' = 2y + e 3x . 


158: y' = x x y + x. 

159: s sen t + s cos t = 1. 


2.4 Equações exatas e fatores integrantes 

160: Uma equação 

p{x,y) dx + q(x, y) dy = 0 

(ou p + qy' = 0), com p, q definidas em uma região U Ç IR 2 , -é chamada 
exata se existir /: U —> IR diferenciável tal que 

df 9f , rr 

— - = p e 7A = q em todo U . 
ax oy 

161 Note bem: Em alguns momentos, y é a variável com respeito à qual 
se toma uma derivada parcial; em outros, y é uma função de uma única 
variável x e tem sua derivada usual. Quando derivamos / quanto a x, 
tratamos apenas da derivada parcial quanto a essa variável, não quanto à 
parte implícita em y, exceto em 164, onde calcularemos a derivada total. 

162: Suponha que p , q são de classe C 11 em U. Se p dx + q dy = 0 é exata 
com / como acima, então 

dp 

dy V/Oy dx' 

Se U é simplesmente conexo (como IR 2 e qualquer retângulo ou disco), então 
a identidade é suficiente para exatidão, como se mostra no Cálculo 

de várias variáveis. 

163: Portanto, caso verifiquemos que 

dp dq 
dy dx 

para uma equação qualquer, vale a pena procurarmos a função / (que pode 
existir ou não). Como veremos nos próximos exemplos, para determiná- 
-la, calculamos' f = f q dx ou f = f p dy, tratando a outra variável como 
constante. 

164: A equação f(x,y) = K define soluções implícitas da equação exata, 
sendo / como acima, porque 

í+w' = i + !£•£ = £/<*.»<*» = <>• 
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Ov 


d\[ 



Note que, dada a condição inicial y(x o) = yo, temos K = f(xo,yo)- 

165 Exemplo: (3a; 2 + 4a :y) dx + (2 y + 2a: 2 ) dy = 0, y( 1) = 0. 

Identificamos p = 3x 2 + 4a :y e q = 2y + 2x 2 e verificamos que e 
são iguais (ambos valem 4a:), então tentaremos encontrar /. 

Já que jJ- = p, integramos 


/ = J pdx = x 3 + 2x 2 y + A(y). 


Note que integramos p com respeito a x, então aparece uma constante de 
integração A, mas A é constante apenas com relação a x, podendo depender 
de y, de modo que escrevemos A = A(y). Também poderíamos experimen¬ 
tar / = Jqdy + B(x). 

Agora utilizamos a outra condição: ^ = q reescreve-se 
2a ; 2 + A' (y) =2 y + 2x 2 , 

donde A'(y) = 2y e então A(y) = y 2 + C. Assiffe, obtemos 

/( x, y) = x 3 + 2x 2 y + A(y) = a^-l- 2x 2 y + y 2 + C. 

Verificamos que, realmente, ^ = p e ^ = q\ afinal, fizemo-lo para que 
funcionasse! 

Fazendo / = K, obtemos a : 3 -Ç%x 2 y + y 2 = D, onde absorvemos as 
constantes. Com o ponto inicial (1,0), calculamos D = 1, de modo que 

» 3 + 2x 2 y + y 2 = 1 

e tem-se y = —x 2 ± y/.x A — a : 3 + 1 pela fórmula quadrática. Para determinar 
o sinal, conferimos novamente com o valor inicial y(l) = 0 , satisfeito apenas 
por 

y = —x 2 + \/ x A — x 3 + 1 . 

166 Exemplo: (g^% -|- 3x 2 y 2 ) dx + (2a : 4 + 2x 3 y) dy = 0. 

Temos p = S>x 3 y + 3x 2 y 2 e q = 2x A + 2x 3 y, com e || iguais (são 
8 a : 3 + 6x 2 y). Integramos 

/ = J pdx = 2 x 4 y + x 3 y + A{y); 


impomos 


= q => 2a ; 4 + a ; 3 + A'(y) = 2x 4 + 2x 3 y 
=> A\y) = 2 x 3 y — x 3 => 

=> A(y) = x 3 y 2 - x 3 y + C; 
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obtemos 


/( x, y) = 2a A y + x 3 y + A(y) = 2x 4 y + x 3 y 2 + C; 


escrevemos / = K , ou seja, 

2 a: 4 y + x 3 y 2 = D , 

de onde resolvemos y = — x ± x~ 2 y/x e — Dx. 

No próximo exemplo, escolhemos a integral de q quanto a ú, que é mais 
simples: 

167 Exemplo: (x 2 + 1 )y' = x 2 — 2xy. 

Temos p = 2 xy — x 2 e q = x 2 + 1 com e iguais (são 2a:). Para 
achar /, fazemos 


/ = qdy = {x 2 + 1 )y + 

J 


já que p = , obtemos 


2a;y — x 2 = 2 xy + B'( x), 


donde B'(x) = —x 2 e B{x) = —x 3 /3/\C. Então 


/( x, y) = {x 2 + % - x 3 /3 + C 
e a solução é y = (D + x 3 )/J^m 2 +]$■ 

L C* V 


Verifique a condição necessária para exatidão e resolva pelo método 


exposto acima: 

168: ydx + xdy = 0. 

169: y cos x dx + sen %dy = 0 . 
170: 3 e 3x y dx + e 3x dy = 0. 


171: 2xy 2 dx + 2x 2 ydy = 0. 
172: y 1 = Sa:- 1 ^. 

173: y’ = —(y + l)/(a: + 1). 


O método dos exemplos e exercícios pode ser aplicado sempre; demons¬ 
traremos um TEU de modo a ver como a solução é dada diretamente pelo 
Teorema da Função Implícita: 

174 TEU: Sejam p, q: R —1 IR de classe C 1 , onde Ré o retângulo dado por 
a<x<bec<y<d, tais que Suponha ainda que q nunca se 

anule em R. Então por cada ponto (xo,yo) G R passa uma única solução 
de pdx + qdy = 0. 
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Solução: Como R é simplesmente conexo, então a equação é exata em 
R, ou seja, existe /: R — > IR de classe C 2 tal que || = P e = q- Veja 
que ^{xg,yg) = q{xg 1 yg) 0: pelo Teorema da Função Implícita, existe 
um intervalo / contendo xg e uma função y: I —> IR tal que f{x,y{x)) — 
f(x 0 ,y 0 ) para todo x £ I. 

Já vimos que essa função deve ser solução de p + qy' = &T Quanto ao 
valor inicial, temos f{xg,y{xg)) = f(xg,yg)', como jd = q, Contínua sem 
zeros, a função f(x o, •) é estritamente monótona em yg e então y{x o) = yg. 

Para demonstrar a unicidade, suponha que y\ , y% sejam soluções do 
mesmo PVI. Então q(yi — y 2 )' = 0 subtraindo-se as equações. Porque 
q não se anula, temos (y± — 2 / 2 / = 0 e então y\ — yi é copstante, ao longo 
de um intervalo, e igual a y-\ (xg) — 1/2 (^ 0 ) = 0 . 

Agora, passamos ao segundo tópico desta seção: 

175: Diz-se que y: U — > IR é fator integrante de pdx + qdy = 0 se 

(yp) dx + (yq) dy = 0 for exata em U. 

Como multiplicar uma equação não a altera, qualquer solução dessa equação 
exata será solução da original. 

176: Se y é fator integrante, então 

<Km>) _ d(yq) 
dy dx 

O que se basta saber sobre U, y, p, q para isso e para a recíproca? 

177: Mostre que se 


dp _ dq_ 

dv dx 



só depende de x r) então a equação admite um fator integrante y que só 
depende de x e tem-se também que 

dy 

— = yip. 
dx 

Sugestão: Usando-se a regra de derivação do produto em 176, chega-se 
a uma equação diferencial parcial de y. 

178: No caso de tp e, portanto, y dependerem apenas de x, resolvemos a 
equação de variáveis separáveis y' = yip obtendo uma fórmula explícita 
para o fator de integração: y(x) = exp(— f ip{x) dx). 
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Escolha a primitiva (ou a constante de integração) mais simples; por 
que qualquer uma serve? 


179: Analogamente, se 


dg _ dp 

dx dy 


depende apenas de y, um fator integrante fj, pode ser obtido como solução 
de 

dy a 

— = yd. 

dy 

180 Exemplo: (8 xy + 3 y 2 ) dx + ( 2x 2 + 2 xy) dy = 0, y{^>y 1. 

Solução: Temos p = 8 xy + 3 y 2 , donde = 8 aj»+ 6 y, é q = 2x 2 + 2 xy, 

da 

donde -rff = 4x + 2 y; então a equação não é exata. Porém, 

% ~ !f _ 8 x + 6y - 4x - 2/^ 2 
q 2x 2 + 2 xy x 

e a equação tem um fator de integração y que só depende de x. Para 
encontrá-lo, pomos 

dy 2 

d V ^ x ’ 

donde y(x) = Cx 2 por separação de variáveis; escolhamos (7=1. 

A nova equação é 

(8x 3 y + 3 x 2 y 2 ) dx + ( 2x 4 + 2x 3 y) dy — 0 

e é exata: nós a resolvemos em 166. Lá, determinamos f(x, y) = 2x 4 y+x 3 y 2 
com jd = p e ^ = q. Já que, no ponto inicial, /(2,1) = 40 (ponto inicial), 
a solução é dada implicitamente por 

2 x A y + x 3 y 2 = 40; 


resolvendo em jKe escolhendo a solução que satisfaça y(2) = 1 , obtemos 
y(x) = — x + x% 3 Vx 8 + 40.t 3 . 

181: ( x 2 + y 2 ) dx + ( x 3 + 3xy 2 + 2xy) dy = 0 . 

Indicação: Com p, q correspondentes, verifique que essa equação não é 
exata, que 

S - a 

q x 3 + 3 xy 2 + 2 xy 
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depende de ambos x , y , mas que 


dg _ dp 
dx dy _ g 

V 

e isso independe de x, ou seja, apesar de constante, “depende apenas de y”. 
Como solução de ^ = 3/x, tome /x = e 3y . Confirme que 

e 3y {x 2 + y 2 ) dx + e 3y {x 3 + 3 xy 2 + 2 xy) dy = 0 
é exata e resolva-a. 

Determine um fator integrante e resolva: 

182: (3y 3 — x) dx + 3 xy 2 dy = 0. 186: x = ( x 2 — lj/tx. 

183: (3 xy — Ay) + (2x 2 — Ax)y' = 0. 187: (aijÁsh 5) dx + 3 x 2 y dy = 0. 


184: y sec x dx + sen xdy = 0 . 


o 2 2 

185: y' = 3y ~ x . 


188: (x + y + 4) + (3x + 1 )y' = 0. 
189: y’ = =^= 2 . 


2.5 Substituições 

Estudaremos, nesta seção, outras formas de equações que podem ser 
transformadas em, ou estudadas eom, equações das formas básicas que já 
conhecemos. Note que mudaremos o significado de algumas letras. 

Começamos com as equações “homogêneas” e apresentamos, a título de 
exemplo para as demais substituições, um estudo teórico de redução aos 
resultados obtidos para equações de variáveis separáveis. 

190: Uma equação, jry= y(x,y) é homogênea se vale g{ax,ay) = g{x,y) 
para qualquer a ^ 0 . 

Atenção: Conhecemos na Seção 2.3 as equações lineares homogêneas, 
que são bem diferentes. De fato, o adjetivo “homogêneo” costuma-se referir 
mais a elas do que a estas que vemos agora. 

191: Mostre que toda equação homogênea pode ser reescrita na forma y' = 
f(y/xkl e reciprocamente. 

Cuidado: O que acontece com x = 0 ? 

192: Mostre que a substituição intuitiva z ~ y/x transforma a equação 
y’ = f(y/x) na de variáveis separáveis z 1 = x~ 1 (f(z) — z), onde z{x) será a 
função incógnita. 
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193 Exemplo: xyy 1 = x 2 — y 2 . 

Se x = 0 ou y = 0, vemos que então ambos x, y = 0; isolemos y'. 

Veja que 

, _ x 2 - y 2 _ 1 - (y/x) 2 ' 
xy {y/x) 

o primeiro quociente mostra que a equação dada é homogênea pela definição 
e o segundo, que coloca y' em função de y/x , pode ser obtido dividindo-se 
numerador e denominador por x 2 . 

Efetuamos a substituição indicada, seja y = xz no primeiro quociente 
ou y/x = z no segundo, e separamos as variáveis: 

. 1 — z 2 . 1 — 2 z 2 zdz 

Z + xz = - => XZ = - --—^7 = dx. 

z z 1 — 2 \^_ J 

Integramos e obtemos 2 : 

— \ ln(l — 2 z 2 ) = x + C => z = ke~ 4x ). 

Não nos esquecemos de determinar y: 

y = xz = ±^(C- ke~ Ax ). 

194 Exemplo: y' + 1 = A 1 + y/x. 

Essa equação quase já está na forma y' = f{y/x) e requer 1 / 0, isto é, 
x < 0 ou x > 0 . 

Pondo z = y/x, obtemos zM*, : xzf+ 1 = \A + z, ou seja, 

dz dx 

— (1 z) X 

Isso pode ser integrado çpm a substituição u = \/l + z (ou z = u 2 — 1), 
resultando em 

—2 ln 1 1 — xx | = ln \x\ + Cq, 

donde 

i-JL 

\/N' 

K 2 2 K 

M VW 

í K 2 — 2 Kyfx se x > 0 , 

1 — K 2 + 2 Ky/^x se x < 0 . 
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Relembre no ponto 69 como concluir essa resolução quando x > 0. No 
caso x < 0, vem analogamente y{x) = —C — 2\/—Cx para x < 0 ou 
y(x) = C + 2 y/Cx para x ^ C. 

Mostre que estas equações são homogêneas e resolva-as: 

195: xyy' = x 2 + y 2 . 196: y' + 1 = — \/l + y/x. 

197: Suponhamos / definida em ]a, b[. Então a função f{y/x)% stá definida 
entre as retas y = ax e y = bx. Esse domínio G deverá conter o gráfico de 
qualquer solução de y' = f(y/x): 



198: Deduza um TEU específico a partir daquele de separação de variáveis: 
Suponha que (£ 0 , 2 / 0 ) G Gr-^que / é contínua com f(u) ^ u para todo 
u € ]a, ò[; então existe uma única solução de y' = f(y/x) com gráfico em G 
passando por esse ponto. 

Eis outras substituições interessantes: 

199: Considere a equação y' — f(ax + by) com a,b constantes reais. Mos¬ 
tre que a substituição intuitiva z = ax + by transforma essa equação na 
autônoma z' = a + bf(z). 

Resolva estas equações por esse método: 

200 : y' — ^Sen(4ai — 2 y )) -1 + 2 . 

Solução: Com z — 4x — 2y temos 

z' = 4 — 2y' = 4— 2[(senz ) -1 + 2] = (senz)^ 1 , 

donde sen zdz = —2 dx. Então — cos z = —2x + C e concluímos que 

4x — z , . . 

y = —-— = 2x — 2 arccos( 2 x — C). 
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201: y' = 2 + (3x - 5y) 2 . 

202 Bernoulli: Dada a equação 

y' = f(x)y + g{x)y k , 

onde k ^ 1 e as funções f,g são contínuas, mostre que a substituição 
z = y l ~ k transforma-a na equação linear 

z' = (1 - k)f(x)z + (1 - k)g(x). 

Isso foi observado por Leibniz. 

Resolva as seguintes equações, servindo-se da observação de Leibniz: 

203: y' = ay — by 2 (essa é uma forma da equação de Verhulst que já 
resolvemos por separação de variáveis). 

Indicação: Temos k = 2 e obtemos z' = — az -f b, que é linear (mas, 
como a, b são constantes, também é autônoma). Após os cálculos devidos 
para sua solução, obtemos 


b - Ce ~ ax ' 

204: y' = x~ x y + y 3 . 205: y' = xy + 2y K . 

206 Riccati: Suponha que y p é uma solução particular já conhecida de 

y' =-ftÁ+ii x )y + h ( x )y 2 , 

por exemplo, se / = 0 então y p = 0. Mostre que a substituição y = y p + z -1 
transforma-a na equação linear 

z' =<[-g(x) - 2y p (x)h{x)]z - h{x). 

207 Ortogonalidade: Suponha dada uma equação de primeira ordem na 
forma geral 

F(x,y,y') = 0. 

Cada solução dessa equação descreve, com seu gráfico, uma curva no plano; 
assumindo um TEU, essas curvas são todas paralelas e particionam o plano. 
Forme agor^á equação 

F(x,y,-í/y') = 0 

(para a mesma F), ou seja, simplesmente substitua cada ocorrência de yf 
literalmente por —l/y'. As curvas descritas por essa equação, sob certas 
condições, também foliam o plano. 
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Mostre que essas duas famílias de curvas são ortogonais, isto é, qual¬ 
quer curva da primeira família e qualquer curva da segunda família, nos 
pontos onde se encontrarem, têm tangentes perpendiculares. Por exemplo* 
as linhas de força de um campo elétrico e as linhas equipotenciais, em um 
plano contendo cargas dispersas, devem ser ortogonais. 

O procedimento de substituir y' por — 1 /y' é chamado determinação das 
trajetórias ortogonais. 

Finalmente, apresentamos um método de substituição intuitivamente 
óbvio quando um problema pode ser melhor formulado, entendido ou resol¬ 
vido em outro sistema de coordenadas: 

208: Dadas as fórmulas usuais para coordenadas polares versus retangula¬ 
res, 

x = r cos 9, y = r sen 9 , 

9 = arctg (y/x), r = \J x 2 + y 2 , 

compute 


dx = cos 9 dr — r sen 9 d9, dy = sen 9 dr + r cos 9 d9 

, _ xdy — ydx u xdx + ydy 
dO = — . , dr = — , 1 ' . 

\J x 2 + y 2 \Jx 2 + y 2 

Indicação: Pela Regra da Cadeia, por exemplo, dx = ^ dr + ^| d9. 

209: Converter a equação y' = f(x,y), de que se buscam soluções y(x), 
para coordenadas polares significa escrever 


sen 9 dr + r cos 9 d9 
cos 9 dr — r sen 9 d9 


= f(rcos9,rsen9) 


por substituição, em que se buscam soluções r(9). (Podemos sempre sim¬ 
plificar |r| =<»*'.: por quê?) 

210 Exemplà: ([Dem]) (a; 2 + y 2 ) dx — xy dy = 0. 

Por substuição, temos 

r 2 [cos 9 dr — r sen 9 d9] — r 2 sen 9 cos 0[sen 9 dr + r cos 9 d9] =0. 

Isso se simplifica como 

r 2 cos 3 9 dr — r 3 sen 9(1 + cos 2 9) d9 = 0, 
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que é de variáves separáveis: 


= sen 9 


' 1 1 
v cos 3 9 cos 9 


Sua solução é r(9) — C , exp(l/2cos 2 9)/\ cos0| com C > 0. Para voltar às 
coordenadas retangulares, substituímos algumas expressões: 

r = Cexp(i(tg 2 6» + 1)) • -f- 


donde 


ln \x\ - Ki + - 02 + x ) - k 2 + 


Assim, a solução é y = ±\/2a: 2 (ln |x| + K). 


2.6 Outros assuntos 

Nesta seção, damos destaque a aplicações, problemas e exercícios que 
envolvem equações de primeira ordem; aqui os agrupamos por assunto, 
sem os dispersar nas diversas seções de métodos de resolução. Também 
conheceremos alguns métodos numéricos e gráficos para aproximação ou 
esboço de soluções quando sua resolução formal é difícil. 

Relembramos que, em Física, as constantes de proporção (como a de 
Avogrado ou a de Planck, sejam adimensionais ou não) são sempre positivas 
e a multiplicação por — 1, quando ocorre, sempre é indicada explicitamente 
com o sinal negativo. 

211 Esfriamento: Newton formulou a seguinte lei sobre o esfriamento 
natural de um corpo previamente aquecido: A variação da temperatura do 
corpo é proporcional à diferença entre essa temperatura e a do ambiente; 
seja a o valor absoluto dessa constante de proporcionalidade. 

Com atenção aos sinais, determine a equação diferencial que rege a 
temperatura-ST do corpo, em função do tempo t, sendo a temperatura do 
ambiente igual a A(t). 

212: Suponha que a situação descrita tenha início em t = 0, com T( 0) = To 
(PVI). .,Npte que, portanto, a equação e suas soluções somente são válidas 
para>£ p 0. O corpo, no “passado” t < 0, não esquenta indefinidamente 
porque nem a equação tem validade, nem a situação era a mesma, nem a 
lei de esfriamento aplica-se. 
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213: Suponha agora que o corpo é originalmente mais frio que o ambiente 
e, então, que ele se aquece gradualmente; assuma uma lei análoga à do 
esfriamento, com a mesma constante de proporcionalidade. Mostre que 
a mesma equação vale também para esta segunda situação, estudando os 
sinais de ambos os lados. 

Desse modo, a equação é válida qualquer que seja a comparaçãosentre 
T e A ou mesmo caso seus gráficos se cruzem. 

214: Identifique a forma dessa equação e resolva o PVI T(0) = T 0 . 

215: Simplifique a solução sob a hipótese de A ser constante (por exemplo, 
em um frigorífico) e observe que a convergência é exponencial¬ 

mente rápida, mas somente se dá com t —> oo, isto é, o equilíbrio térmico 
nunca é atingido em tempo finito ou, aproximadamente, só ocorre após 
certo tempo. 

216: A lei de Newton é geralmente formulada apenas para A constante. An¬ 
teriormente, extrapolamos sua validade usando a equação diferencial para 
representar somente uma “fotografia” da dinâmica envolvida: a derivada é 
uma velocidade instantânea, não média. Porém, a adequação dessa extra¬ 
polação à realidade, ou seja, a confirmação dessa teoria, requer experimen¬ 
tação em ambientes de temperatura variável! 

217: Suponha que dois corpos, com temperaturas T\(t),T 2 {t), são colocados 
em contato, formando um sistema termicamente isolado, no instante t = 0 ; 
sejam Tio,T 2 o os valores iniciais. Podemos prosseguir com a extrapolação, 
assim, supondo que cada çórpo atua como “ambiente” para o outro. Escreva 
as duas equações diferenciais que regem Tj e Tj, assuma o mesmo valor de 
a para ambas as equações e mostre que 

(a) a soma Tj + Tj é constante com valor Tio + T 20 ; 

(b) não há mudança de temperatura (ou seja, as variações são zero) pre¬ 
cisamente quando T\ = T 2 ; 

(c) essa temperatura de equilíbrio vale |(Tio + Tjo). 

218: Refaça o ponto anterior usando constantes diferentes aq,a 2 para as 
equações è mostrando que 0 : 2/1 + oiTj é constante. Qual é a temperatura 
de equilíbrio? 

Mostre ainda que se toi,Ci são a massa e o calor específico do corpo 
de temperatura Ti e, analogamente, « 1 - 2,02 referem-se a Tj, então devemos 
tomar a 2 = «iiCi e oq = « 1202 - 
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219 Permeabilidade: Assumiremos esta lei: A quantidade de radiação 
absorvida ao passar por uma lâmina fina de chumbo é proporcional à quan¬ 
tidade de radiação incidente e à espessura da lâmina. 

As razões que levam a essa formulação são as seguintes: Se cada partí¬ 
cula radioativa tiver uma probabilidade independente e de ser absorvida, 
então a absorvição é proporcional à incidência. Porém, e é muito pequena 
para uma lâmina fina, porque as partículas têm grande energia. Se supu¬ 
sermos uma segunda lâmina idêntica atrás da primeira, a proporção das 
partículas que incidem sobre essa segunda lâmina é 1 — e e a daquelas que 
a atravessam é, novamente, 1 — e desse número, ou seja, 

(1 — e)(l — e) = 1 — 2 e + e 2 

das partículas originais. Como e 2 é muitíssimo pequeno, podemos des- 
prezá-lo: ao dobrar a espessura da lâmina fina, simplesmente dobramos a 
probabilidade de absorção. 

A mesma lei descreve diversos fenômenos de permeabilidade. 

220: Em termos formais, a quantidade de radiação absorvida é alh, onde 
a é uma constante positiva de proporção, / mede a incidência radioativa e 
h é a espessura fina da lâmina de chumbo. 

221: Se metade das partículas radioativas incidentes são absorvidas por 
uma certa chapa padronizada, quanto de radiação atravessará uma barreira 
com três vezes sua espessura? 

Veja: Não faça uma simples regra de três, porque as chapas considera¬ 
das não são lâminas finas! De fato, a barreira não poderá absorver 150% 
das partículas incidentes... 

Sugestão: Divida a barreira em três chapas idênticas à padronizada e 
use a fração que atravessa cada chapa como total incidente para a próxima 
chapa. 

222: Esse método sugerido permite responder a mesma pergunta para uma 
barreira de espessura qualquer? 

223: Para uma barreira de chumbo de espessura arbitrária, seja I(x) a 
quantidade de radiação que atravessa essa chapa até a profundidade x. 
(Pomos /(0)‘ = Iq, que é o total incidente.) Então, para h pequeno, a 
quantidade que chega a x + h é aquela I(x) menos a quantidade absorvida 
na região entre x e x + Z, que corresponde a uma lâmina fina. 

Deduza a equação diferencial da permeabilidade e resolva-a; verifique 
sua concordância com o cálculo em 221 . 
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Nos próximos pontos, estudaremos alguns circuitos elétricos simples, 
formando e resolvendo equações diferenciais que regem a corrente elétrica 
nesses circuitos. Embora tentemos motivar algumas dessas equações, não 
daremos ênfase nem às leis físicas subjacentes, nem à análise física posterior 
— como o cálculo da potência desenvolvida ou da energia dissipada pelo 
efeito Joule, ambos muito curiosos —, porque esses estudos cabem melhor 
em um texto especializado como [Tip]. 

224 RC: Este circuito liga, em série, uma chave interruptora, um resistor 
com resistência R e um capacitor de capacitância C , ambas R, C constantes 
positivas. Quando o circuito fechar, a corrente elétrica sempre satisfará 
I = U/R = Q/RC: é a carga Q armazenada no capacitor (±Q em cada 
lado) que gera a diferença de potencial U no circuito. Suponha que, no 
instante t = 0, tenhamos Q(0) = Qo e então a chave seja fechada: a corrente 
inicial é /o = Qo/RC e, conforme t aumenta, o capacitor descarrega-se e 
Q diminui. Como corrente significa carga deslocada por unidade de tempo, 
temos / = —Q, onde o sinal indica a diminuição de Q. 

Assim, 

Q = - Q/RC , Q(0) = Qo. 

Resolva esse PVI, obtendo a carga e (por derivação) a corrente como funções 
do tempo. Graficamente, como se comporta a descarga do capacitor? 

225: Suponha agora que o capacitor esteja descarregado (Qo = 0) e que se 
insira nesse circuito, em série, uma fonte de força eletromotriz, que origina 
uma diferença de potencial, constante E. Assim, o circuito formado é de 
corrente contínua quando se fecha a chave em t = 0 e inicia-se a carga do 
capacitor. Agora, I = Q k,E = Ur + Uq = RI + Q/C, de modo que 

Q = E/R - Q/RC. 

Resolva o PVI dessa equação com Q(0) = 0, determinando a carga e a 
corrente como funções do tempo e traçando seus gráficos. Qual é a corrente 
inicial? 

226: Suponha que o capacitor esteja carregado, ou seja, Q 0 ^ 0: o sinal de 
Qo depende da polaridade da conexão; conforme for positivo ou negativo, 
a carga aumentará sempre ou primeiro diminuirá (em valor absoluto), para 
anular-se, e então crescerá com polaridade oposta. 

Resolva esse novo PVI e determine graficamente o comportamento das 
soluções dependendo do sinal de Q 0 . 

227: Qual será a carga completa do capacitor? Em quanto tempo obteremos 
90% dessa carga? 
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228: Estude o circuito RC com uma fonte de força eletromotriz alternada, 
isto é, que origina a diferença de potencial E(t) = Eq cos(ojot): 

RQ + Q/C = E 0 cos(woí). 

229 LR: Ligamos, também em série, uma chave interruptora, uma fonte 
de força eletromotriz E, um resistor com resistência R e um indutor de 
indutância L, todas E,R,L constantes positivas. 

O indutor é uma espira ou solenoide que, sob a passagem de corrente 
elétrica, gera um campo magnético. Se a corrente variar, a variação cor¬ 
respondente do campo opõe-se à variação da corrente, gerando uma força 
contra-eletromotriz LI. (Se I cresce, então o indutor atuàppmo um resistor 
e a queda de potencial é LI > 0; se I diminui, então o indutor atua como 
um gerador e a “queda” é LI < 0.) 

Então E = Ur + Ur ~ RI + LÍ. Resolva o PVI 

LÍ + RI = E, /(0) = 0, 

determinando a corrente em função do tempo. 

230: Determine o valor de equilíbrio Jfi na i da corrente no circuito, quando 
a força contra-eletromotriz é nula. Essa é a corrente final do circuito. 

231: Estude o circuito LR com fonte de força eletromotriz alternada: 

LÍ + RI = E 0 cos(woí). 

Retornaremos a esse tópico na Seção 3.7, a respeito dos circuitos LC e 
LCR, quando pudermos estudar as equações de segunda ordem pertinentes. 
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57: A = 3, B = —3/2. 72: y = y'x: substitua y' = C. 73: y' = 8. 
74: y" = 0. 75: y'" = 0. 76: y' = 2 y/x: multiplique y' = 2 Cx por x e subs¬ 
titua y. 77: y' = 3y. 78: y" = —25 y. 79: y' = —x/y. 83: Derive uma vez* 
substitua 2e* y por x e derive novamente. 88: Para o Teorema da Inversa, 
tome f(x) = x 2 : vem /'(O) = 0; por menor que seja U, não há bijeção. 
94: Sol.: y = ±\/xq + y 3 — x 2 . 95: y = Csena;. 96: y = G/x. 97: y = 
± \J2 ln \x\ — x 2 + C. 98: y = ln \e x + C|. 99: y = A\JC/(x 2 + 1) — 1. 
100: y = ln11 — C/e x \. 101: y = (Cx— l)/(x + l). 102: y = tg(ln \x\ + C). 
103: y = 8/x. 104: y = 0. 105: y = (7x - l)/(x + 1). 106: y = x. 
107: y = (x ln x — x + 2) 2 . 108: y = (x + l)/-\/2 + 2a; 2 X115: g(y ) = y" 1 
definida em ]—oo,0[ e ]0,oo[, onde nunca se anula; então o TEU aplica-se 
em cada intervalo. Sol.: y = ± \/x\ + y^ — x 2 . 116: g(y) = y 2 definida em 
IR; não se anula em ]—oo,0[ e ]0,oo[, onde se aplica»gJTEU; sol. regulares 
y = —(x 2 — Xq — Vq 1 ) 1 não incluem sol. singular y <= y 0 = 0 e nunca a 
interceptam (porque nunca se anulam). Isso prova um TEU geral. (Note: 
a forma y = yo/(1 — y^x 2 ) inclui todas as soluções.) 117: g(y) = y definida 
em IR; não se anula em ]— oo,0[ e ]0,oo[, onde se aplica o TEU; soluções 
regulares y = xyo/xç, incluem sol. singular y = y 0 = 0 e nunca a intercep¬ 
tam quando j/o / 0 (porque x Z 0). Isso prova um TEU geral (mas para 
x < 0 ou x > 0). 118: g(y ) = 1 + y 2 definida em IR e nunca se anula; 
o TEU aplica-se sempre. Sol. y = (yo A tg(a; — Xq))/(1 — y 3 tg(a; — Xo )). 
119: g/y) = y 2 / 3 definida em IR; não se anula em ]—oo,0[ e ]0,oo[, onde se 
aplica o TEU; sol. regulares ,y = yj'’ 1 + (x — a;o) 3 não incluem sol. singular 
V = Vo = 0, mas a interceptam em x = Xq, de modo que não há um TEU ge¬ 
ral. Note ainda: as sol. podem bifurcar, como y = 0 e y = { ^ 

ou confluir. 120: g(y) = 1 definida em IR e nunca se anula; o TEU apli¬ 
ca-se sempre (mas naras j; </ () ou x > 0). Sol. y = y 3 + x^ 1 — a; -1 . 
123: Se y 1 satisfaz y' = g(y) e y( 0) = y 0 , então y 2 (x) = y\(x — x 0 ) satisfaz 
y' = d(y) e ^(^o) = yg- 128: 0 instável e 1 estável. 130: 2 semiestável. 
132: Use o TEU e 123. 135: x = (a:/ -1 — aí) -1 . 138: Dica: obtenha 
(x — M)x~ x = IV- Mj -1 e isole esse x. 139: Dica: como exp (—Mat) > 0, 
o denominador é sempre > Xq > 0, donde 0 < x < M. 146: E preciso, 
primeiramente, trazer a população a x, aguardando seu crescimento sem 
extração caso Wo < x. 147: Imponha x = 0. 148: Basta estudar o “saldo 
viabilizanto®’,’ y = x — L, donde y = a(M — L — y)y. 150: r(t) = Ro — kt. 
156: y= (K — cos 2x)/ cos x. 157: y = e 3x + Ke 2x . 158: y = Kx + x 2 . 
159: 5 — (í T K )/ sení. 169: y = k/senx. 171: y = ±kx. 172: y = Cx 3 . 
173: y = (C — x)/(x + 1). 182: Fator a; 2 ; sol.: y = [(x A — K )/4x 3 ] 1 / 3 . 
183: Fator A'|4 — 2a;| -1 / 2 : simplifique com K = \/2 e escreva-o |aa — 2| 1 / 2 . 
Integre df/dy = q “novo”. Sol.: y = C/x\J\x — 2|. 184: Fator seccc; sol.: 
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y = C/ tgx. 185: Fator a; -4 ; sol.: y = ±\/x 2 — Cx 3 . 186: Fator f -3 ; sol.: 
x = ±\/l — 2 Ct 2 . 187: Fator a: -4 / 3 . 188: Fator (3a; + l) -2 / 3 . 189: Fa¬ 
tor (a: + l) -2 ; sol.: y = (x + l)(ln|a: + 1| — C) — 1. 191: Dada g , tome 
f(s) = g(l,s) e use a = a: -1 ; dada /, tome g(x,y) = f/y/x). 192: Derive 
y = xz e iguale a f(z). 195: y = ± x^2 lna : + C. 196: Analog. a 194^ 
y = C — 2\JCx para 0 < x < C e y = —C + 2'J—Cx se — C < x < 0. 
201 : y = 3a;/5 — \fl tg (—xVih + C)/202 : Calcule z' e isole y'\ subs¬ 
titua y k por y/z e multiplique toda a equação por z/y. 204: y - '/C/x 2 - 
2a;/3) -1 / 2 . 205: y = [Cexp(a;(l - 7r)) - 2a; - 2/(1 - tt)] 1 /^ - ^. 211: f = 
—a(T — A). 214: Linear: por variação da constante, T(t ,) = e~ at To + 
e ~ at J* ae as A(s) ds. 215: T/f) = A+e~ at (To—A). 218: A energia térmica 
do sistema miCiT! + TO 2 C 2 T 2 é constante. 221: Um oitavo. 223: Temos 
I(x + h) = I(x) — al(x)h. Divida por h e tome h —> 0: I/j= —ai. Solução: 
I(x) = Ioe~ ax . Assumindo 1(1) = Iq/2, vem a =^^2 e 1(3) = Iq/8. 
224: Q(t) = Qoexp(—t/RC) e I(t) = Ioexp(—t/RC). 225: Q(t) = 
EC( 1 - exp (-t/RC)) e I/f) = E exp(-t/RC)/R. 226: Q(f) =EC- 
(EC —Qo) exp(—t/RC) e I(t) = (EC-Q 0 ) exp(-t/RC)/RC. 229: I(t) = 
E( 1 — exp (—Rt/L))/R. 230: /fi na i = E/R. 260: O PVI de segunda 
ordem especifica os valores da solução % de sua derivada em um ponto. 
Se duas soluções têm o mesmo valor-em um ponto, deverão ter deriva¬ 
das distintas aí. 271: L.D.: 1(0) + 0(5a^f + 0(—2a; 3 ) = 0. 272: L.D.: 
3(2sen(a; 2 - 1)) + 2(-3sen(a: 2 - 1)) = 0. 273: L.I. 274: L.D.: a; = 
|(a; + 5) + §(a: — 4). 275: L.I. 276: L.I. 277: L.I. 278: L.I. 279: L.I. 
280: L.I. 283: Suponha L.D.: uma das funções escreve-se como combina¬ 
ção linear das demais. Então também uma coluna da matriz é combinação 
linear das demais e o determinante anula-se. 317: y = Cie 2x + C^e 3 ^. 
318: y = C/e 2 * + C 2 xe 2x . 319: y = Cie 3 * + C 2 e~ 3x . 320: y = 

Ci cos3a; + C 2 sen3a;. 321: y = C\e 2x cosx + C 2 e 2x senx. 322: y = C\e x A 
C 2 e 2x . 323: y = C\e~ x + C 2 xe~ x . 324: y = C\e x cos2x + c 2 e x se\\2x. 
325: y = Cie x cos \/3x + C 2 e x sen \/3x. 326: y = C x + C 2 e* + C 3 e~ x . 
327: y = Ci + + C 3 e~ 3x . 328: Note P( 1) = 0 e fatore. Sol.: 

y = Cie* + C 2 e^ + C 3 e~ 2x . 329: Note —1 raiz simples: fatore. Sol. 
y = Cie - * + C 2 e»(j 2 cos(\/3x/2) + C 3 e x ^ 2 sen(i/3x/2). 330: y = Cie* + 
C 2 xe x + C 3 pQy 331: y = Cie^* + C^e^ 3 * + C 3 e -v/3 * + C 4 a;e -v/3 *. 
332: y = C\cos\/3x + C 2 sen\/3x + C 3 xcos\/3x + C<iXsen\/3x. 333: y = 
Ci + C 2 .r f- C 3 x 2 + C 4 a; 3 . 334: y = Ci + C 2 x + Cae^* + C 4 e - ' /2 *. 
335: y = Vi + C 2 cos x + C 3 sen x + C 4 a: cos x + C 3 x sen x. 336: Mostre 
P( A> = (A + 1)”. Sol. y = C i+ la;^ - *. 338: y = -4e* + 2e 4 *. 

339: y = 5e 2 * - 10a;e 2 *. 340: y = 3cos2a;. 341: y = 0. 342: y = 
—£ 2x cosx. 343: y = —^e x —xe x +^x 2 e x . 350: y = Cie*+C 2 e 2 *+e -5 */42. 
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Capítulo 3 


Ordens superiores 


Uma mesma teoria dará conta das equações de segunda, terceira e outras 
ordens, embora seja um pouco distinta do corpo de conhecimento que já 
desenvolvemos para a primeira ordem e refira-se a uma forma mais restrita 
de equações. Poderemos, também, esclarecer alguns dos métodos que já 
vimos e que podiam parecer “mágicos”, como a variação de constantes. 
Precisamos enunciar dois resultados para iniciar nosso estudo: 

260 TEU: Sejam f(x, y 0 , yi ■ ■ •, y n - 1 ) e suas derivadas ..., con “ 
tínuas em uma vizinhança de um ponto (xq, Oo, ai,..., a n -\). Então existe 
ó > 0 e existe uma única função y: ]:ro — â, Xq + á[ —> IR que é solução do 
PVI na forma normal 

y(x o) = ao, y\x o) = ai,..., y (n ~ l \x 0 ) = a n _i. 

261 Superposição: Suponha que y r e y s são soluções, respectivamente, 
das equações lineares Yh^oPíV^ = r e E"=oí ) *2/^ = s onde Po, ■ ■ ■ ,Pn,r,s 
são funções apenas da variável independente x. Mostre que y r ± y s é solução 

de J2i=oPiV$ = r±s - 

E também importante reforçar as diferenças de interpretação em relação 
à teoria de primeira ordem: 

262: Nas condições do TEU, os gráficos de duas soluções da mesma equação 
de segunda ordem podem cruzar-se, mas não podem ser tangentes. Por 
quê? 
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3.1 Dependência linear de funções 

Faremos bom uso dos conceitos de dependência e independência linear 
de funções, de modo que já os apresentamos aqui. A linguagem teórica 
adequada para a formulação do que aprenderemos é aquela básica de Álge¬ 
bra Linear sobre espaços de funções, como em [CoLo]; não há necessidade, 
porém, de aprendê-la antecipadamente. 

Suponha que /i,..., f n são funções sobre o mesmo domínio X e com 
valores reais: 

263: Diz-se que f\..... f n são linearmente independentes caso 

(Vci,..., c„ G IR) | [(Vx G X) J2i=i Cifi(x) = o] =^Si = ... = c n = 0 j, 

isto é, a única possibilidade de termos C\f\ + .+ Qff n = 0 é a trivial. 
264: Diz-se que f„ são linearmente dependentes caso 

(3ci,..., Cn G IR) | [(Vx G X) Ya=i c íM x ) = O] e algum a ^ 0 j, 

isto é, existe um modo não-trivial de termos c\f\ + ... + c n f n = 0. 

265: As funções f \,..., /„ são linearmente dependentes se e somente se pelo 
menos uma delas escreve-se como combinação linear das demais. (Este é o 
tradicional “divida pelo coeficiente.não-nulo e passe para o outro lado”.) 

266 Exemplo: As funções T-ç cos x são linearmente independentes. 

Para mostrá-lo, suponhamos que C \, C 2 G IR são tais que 

ci(x) + C 2 (cos x) = 0 

para qualquer valor de x. Em particular, com x = 0 ou x = 7r/2 obtemos, 
respectivamente, 

Oci + lc 2 = 0, (7 t/2)ci + 0c 2 = 0. 

Esse é um sistema de equações com incógnitas ci,c 2 cuja única solução 
é (0,0). Portanto, pela definição, {x,cosx} é um conjunto linearmente 
independente. 

Note: Foi preciso encontrar valores de x, em número igual ao das fun¬ 
ções apresentadas, de modo que se obtivesse um sistema com única solução 
(0,0), nesse caso em que as funções são linearmente independentes. Não 
há receita para encontrar tais valores, além de explorar os zeros das fun¬ 
ções. Caso as funções fossem linearmente dependentes, o sistema teria uma 
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solução não-trivial, mas seria preciso verificar que a combinação linear cor¬ 
respondente sempre se anula para qualquer valor de x. (Quando tratarmos 
do “wronskiano”, apresentaremos um cálculo simples para decidir sobre a 
dependência linear de funções, que lamentamos que não seja ensinado nos 
cursos elementares de Álgebra Linear.) 

267 Exemplo: As funções 5 (constante), 6 sen 2 x e 2 cos 2 x são linearmente 
dependentes. 

Para tanto, devemos procurar escalares Ci,c 2 ,C 3 , não todos nulos, tais 

que 

Ci (5) + c 2 (6sen 2 x) + c 3 (2 cos 2 x) = 0 

para qualquer valor de x. Tomando x = 0 ou x = | ou x = |, obtemos o 
sistema 

5d + 2c 3 = 0, 5ci + 6c 2 = 0, 5ci • 3c 2 + c 3 = 0, 

que é indeterminado; uma possível solução é dada por c 3 = 6, c 2 = —5 e 
c 3 = —15. Agora, verificamos que realménte 

6(5) — 5(6 sen 2 x) — 15j(á cos 2 x) = 0 

para qualquer x, o que conclui o raciocínio. Também podemos, a partir 
dessa equação, determinar uma função como combinação linear das demais, 
por exemplo, 

(6 sen 2 x) = |(5) — 3(2 cos 2 x). 


Decida se cada conjunto de funções é linearmente dependente ou não: 


268: 0, 5x e -2x 3 . 


273: xee 1 


269: 2sen(x 2 — 1) e —3sen(x 2 — 1). 274: e x e e 2x . 


270: x — 4 e x + 5. 
271: dS- 4, x e x + 5. 


275: senx e cosx. 


276: x e senx. 


272: x e x . 


277: x, e x e senx. 


O resultado a seguir é o mais importante em nosso estudo. Muito mais 
que saber demonstrá-lo, é importante entendê-lo: 
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278 : Considere estas funções: 


e XlX , xe XlX , 


x kl e XlX , 


x km e Xn 


e aiX cos(f3\x), xe aiX cos(/3ix),..., x ll e aiX cos(/3ix), 
e aiX sen(/3ix), xe aiX sen(/3ix ),..., x ll e aiX sen(fiix(), 


e ariX cos(f3 n x), xe anX cos((3 n x ),..., x ln e ar,x còp((3 n x), 
e“" x sen(/3 n x), xe anX sen (j3 n x),..., x ln e ariX sen(/3 n x). 

Assuma que k\,..., k m . h,... ,l n são números naturais, que Ai,..., A m são 
reais distintos e que a\ + (3\i ,..., a n + /3 n i são pomplexos distintos e não 
conjugados, sendo aq, fii ,..., a n , (3 n £ IR. Sob tais condições, o conjunto 
dessas funções é linearmente independente. 

Precisaremos de mais estas definições, para as quais suponha fixadas as 
funções .., /„: X -> IR: 

279: Se ci,..., c n são números reais, então podemos formar a função 

.9 “ Ci/i A A c n f n ] 

diz-se que as funções originais geram (linearmente) g ou que g é combinação 
linear delas. (Note que ci,..., c n são constantes, embora arbitrários.) 

280: O conjunto V de todas as funções geradas por /i,..., /„, para quais¬ 
quer escolhas dos valores reais de ci,...,c n , também é dito gerado por 

elas. Ele pode ser gerado por vários outros conjuntos de funções. Note que 

A. V 

sempre: 

(i) a função constante nula pertence a V ; 

(ii) /fv. . .</„£ V\ 

yjy 

(iii/lÇfeomas e múltiplos de funções em V também pertencem a V ; 

(iv) ou seja, V é fechado sob combinações lineares. 

Assim, V é um exemplo de espaço vetorial. 
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281: Caso sejam linearmente independentes, diz-se que o con¬ 

junto {/i,..., f n } é uma base de V. Seu número n é chamado dimensão 
de V e denotado dim V. (Para isso, mostra-se que n independe da escolha 
das funções geradoras, enquanto independentes.) 

282: Suponha que /iconstituam uma base do espaço V e que 
/: X —> IR seja uma outra função (pertencente, ou não, a V). Então o 
conjunto das funções 

g = Cl/l + • • • + C n fn + /, 

para quaisquer escolhas dos valores reais de ci,..., c n , é simplesmente indi¬ 
cado V + /. Esse é um exemplo de espaço afim e é uma translação de V, 
com a mesma dimensão. (Se / ^ V, ele não é um espaço vetorial.) 

283: Mostre que, nessas condições, toda função / pode ser escrita 

de maneira única como g = c±fi + ... + c n f n + jS ou seja, se supusermos 
também g = b±fi + ■ ■ ■ + b n f n + f então deveremos ter cada bi = Cj. 

Essas definições e a teoria de espaços vetoriais são inspiradas pelos con¬ 
juntos de vetores geométricos. Utilizamos essa interpretação para resumi- 
las, com dimensão n — 2, neste diagrama: 



No restante desta seção, teceremos alguns comentários sobre espaços 
de funções para o*)éstudante que já teve contato básico com a linguagem 
necessária, a título de firmar sua interpretação. 

284 : Em Álgebra Linear, trabalha-se com somas finitas; não há séries infi¬ 
nitas porque não há noção de convergência. Assim, toda combinação linear 
tem um número finito de termos. 

285 : Um conjunto de vetores S é linearmente independente se, para quais¬ 
quer n ^ 1 e iq,..., v n £ S distintos, os únicos escalares C\, ..., c n que 
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satisfazem y',"_ 1 CíVí = 0 são c\ = ... = c n = 0. Observe que n é finito, 
mas arbitrariamente grande. 

286: Mostre que um conjunto de vetores é linearmente independente seve 
somente se todo subconjunto finito seu é linearmente independente. 

287: Mostre que, no caso de S ser finito, basta verificarmos a condição para 
n= |Sj. 

288: O conjunto das funções J\T —>- IR é um espaço vetorial sobre IR, porque 
sabemos somar essas funções e multiplicá-las por escalares (números reais); 
essas operações são feitas ponto a ponto e resultam em funções também da 
forma X —> IR. Esse espaço é indicado IR A . 

Se X é um subconjunto aberto de IR, então cada conjunto C n (X), das 
funções de classe C n com domínio X, é um subespaço vetorial de IR V , 
porque as combinações lineares de funções contínuas, ou n vezes deriváveis 
são também contínuas ou n vezes deriváveis. 

Pelo mesmo motivo, a operação de derivação é uma transformação linear 

C n+1 (X) -A C n (X). 

3.2 Wronskiano 

289: Define-se, para f-\,, ff 1 . b[ —> IR de classe C n ~ 1 , o wronskiano 
de ordem n 

fi(x) ■■■ f n {x) 

fí{ x ) ■■■ fni x ) 

W(f : 

fí n ~ 2 \x) ••• À n ~ 2 \x) 

/i (n_1) (*) ••• â^Hx) 

Note que W(fi ,..., /„) é uma função contínua ]<x, b[ —> IR. 

290: Mostre que, se {/i,...,/„} é um conjunto linearmente dependente, 
então é identicamente nulo. Equivalentemente, se existe 

x £ ]a,D[ tal que W(fi ,..., f n )(x) i=- 0, então {/i,..., /„} é um conjunto 
linearmente independente. 

Veremos em 291 que a recíproca não é verdadeira. Porém, para um 
çpnjunto de soluções de uma certa equação diferencial, mostraremos que 
vale a equivalência “L.D. <t=> W = 0”. 

291: Sejam fi(x) = x 2 e / 2 (x) = x |x|. Mostre que: 
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(a) /i e / 2 são linearmente dependentes no intervalo [0,1]; 

(b) fi e / 2 são linearmente independentes no intervalo [—1,1]; 

(c) o wronskiano de /i e / 2 é identicamente nulo, calculando-o separada¬ 
mente para x ^ 0 e x < 0; 

(d) /i e / 2 não podem ser soluções de y" + g(x)y' + h(x)y = 0 para 
— 1 < x < 1 se ambas g , h forem funções contínuas nesse intervalo. 

292: Assuma que y\,...,y n sejam soluções de 

Pn{x)y (n) + p n -i{x)y {n ~ l) + ... + po(x)y = 0, 

com po, ■ ■ ■ ,p n ■ ]ffl, b[ —> IR contínuas e p n sem zeros. Mostre que, se essas 
soluções são linearmente independentes, então W (iq,, ...,y n )[x) ^ 0 para 
todo x £ ]a, b[. 

Solução: Suponhamos W(yi,... ,y n )(xo) = 0 para algum Xq £ ]a, b[. 
Consideremos este sistema linear algébrico: 

ciyi(xo) + ■ ■ ■ +/min{x o) = 0 
ciy[(xo) + ... + c n y' n (xo) = 0 

ciyí" 1} (x 0 ) + •• ^Fc„?/i n 1) (x o )=0 

Seu determinante é W{y\i... ,y n )(xo) = 0; logo, ele tem pelo menos uma 
solução ci,..., c„ não-trivialW^Çi^ esses valores, definimos 

y = cij/i + ... + c n y n . 

Por superposição (261), y é uma solução da equação do enunciado. Também 
temos 

iMlx o) = cij/^^o) + ... + c n yf£\x o) = 0 

para todo 1 n — 1. Pelo TEU, y = 0, de modo que {yi ,..., y n } é 

linearmente dependente. 

293: Conclua que, com a hipótese em 292 e o fato em 290, se o wronskiano 
anula-se;em um ponto então ele se anula em todo o intervalo de definição e 
as funções são linearmente dependentes. 

294: Novamente na situação de 292, mostre que o wronskiano das soluções 
por sua vez é solução de p n w' + p n _iW = 0. 
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Para 2- ordem: (Temos n = 2, equação p^y" + piy 1 + poU = 0 e 
soluções 2 /i, 2 / 2 -) Escrevemos 

W = W(y 1 ,y 2 ) = V ) V ] = 2 / 12/2 ~ 2/í 2 / 2 , 
l/l 1/2 

de modo que 

W' = i/i i /2 + 1 / 11/2 - i/i yi - i/i i /2 = 1 / 12/2 - 2 /i(ar 

e então 


P 2 W" = 2/1 (P22/2) - 1/2 (P2 2/í) = 

= 2 /l(—Pll/2 - POl/2) - y2(-PíP%-PoVl) = 

= —Pi (2/12/2 - 2 /í 2/2] = -PilE 

Em geral: Mostraremos esse fato nos próximos pontos. 

295 Derivação de determinantes: Seja (/. (J -) uma matriz quadrada de 
ordem n, onde as entradas fij são funções deriváveis em algum intervalo. 
Seja D(x) = det (fij(x)). Então D é derivável e 

D'{x) + ... + D n {x), 

onde Dj é o determinante da matriz obtida de (fij) derivando-se a i-ésima 
linha. 

Solução: O determinante' é uma certa soma de produtos, cada qual, 
quando derivado, torna-se outra soma de produtos, bastando então reorga¬ 
nizar os termos. Simbolicamente, isso pode ser feito a partir da definição de 
determinante e deixamos a cargo do estudante de Álgebra verificá-lo. Eis 
uma alternativá>i 

Consideraremos as entradas /, ? também como variáveis. Seja V t] o 
cofator da entrada f t j: então D = fijVij quando se desenvolve pela 

linha i. Note que V t j não depende de fy, donde dD/dfij = Vij. Agora, 
usamos a Regra da Cadeia: 


=E §- • %=E = t(tm)=t A. 


i=l x i=l 


296 Derivação do wronskiano: Dado W = W(fi ,..., f n ), calculamos 


68 


B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1- quad. 2012 


W' pelo ponto anterior: 


fí ■ 

f 

J n 



h 

fn 


h 

fn 

fí ■ 

f 

J n 

+. 

.+ 

fí ■ 

■ fn 

+ 

fí ■ 

■ fn 

f(n- 2 ) 

An- 2 ) 
Jn 



j(n~ 1 ) 

An- 1 ) 
Jn 


f(n- 2 ) 

f ( 7 — 2 ) 

/í"- 1 ) 

An- 1 ) 

Jn 



/(«-!) 

An- 1) 
Jn 


/í n) A 

An) 

Jn 


Todos esses determinantes, exceto o último, são nulos porque têm linhas 
repetidas. Então a derivada do wronskiano é obtida derivando-se apenas e 
simplesmente a última linha: 

h ■■■ fn 

fí ••• 

W = \ • • : 

An— 2 ) r(n—2) 

J 1 * * * Jn 

An) An) 

J 1 - - - Jn 

297: Verifique que o análogo não vale para derivadas de ordens mais altas. 


298: Para concluir, retomando-se a situação de 292, temos: 



l/l 

2 /n 


yi 

l/n 

PnW = 

2 ) 

* (ri.—2) 

= 




l/i 

’ 2 /n 


2/1 

2 /n 


(n) 

Pn2/l 

(n) 

•, Vnl/n 


\—\n —1 (i) 

- Ej=0 

\—\n —1 (i) 

• -Ei= 0 fti/A 



..l/r • 

Vn 


2/1 

2/n 

= Pn : 

*{n-2) 

An- 2) 

Pn—2 

„(n-2) 

„(n—2) 


Hl 

’ Vn 


2/l 

’ Vn 


(n—1) 

Hl 

(n—1) 

■ 2/n 


(n—2) 

2/l 

(n—2) 

■ Vn 


O segundo termo e os seguintes são todos nulos porque têm repetição de 
linhas. Logo, p n W + p n -iW = 0. 

Resolvamos essa equação explicitamente: 
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299 Abel—Liouville: Fixe xq £ ]a, b[. Por separação de variáveis, a solução 
de p n w' + p n -iw = 0 que passa pelo ponto (xq,w(xq)) é dada por 


j{x) = w(xo) expí — / 

' J X( 


Pn-l(s) 


Utilizando diretamente o TEU de equações de variáveis separáveis, vemos 
que essa é a única solução possível no caso de p n ,p n ^\ contínuas e p n sem 
raízes. 


300: Isso fornece outra demonstração de que o wronskiano das soluções da 
equação em 292 ou é zero em todo ]a, b[ ou nunca se anula em ]a, ò[. 


3.3 Espaço de soluções da equação linear 

Esta seção apresenta a equação diferencial linear de ordem superior e 
caracteriza o conjunto de suas soluções como um espaço afim — ou seja, 
uma translação de um espaço vetorial — _cui a dimensão é igual à ordem da 
equação. Tal resultado teórico é essencial para, posteriormente, podermos 
resolver ou mesmo compreender a resolução dessas equações. 

301: Uma equação da forma 

Pn{x)y (n) +p n - i(a:)2/ ( r S ^- • • ■ + Pi(x)y' + p 0 (x)y = r{x) 


é chamada linear. Assumiremos que as funções-coeficiente Pq, .... p n e a 
função r são contínuas em um intervalo aberto I Ç IR. 

302: Assumiremos que p n nunca se anula em I. Então podemos escrever a 
equação na forma normal 


»<"> = /(W.»<"-')) = 7 1) ~- 

U Pn{x) 


■ --Po(x)y 


Nesse caso, 

, NU í = _P0 df_ = _P]_ df = Pn— 1 

dy pn dy' p n ’dy^ n ~^ p n 

são contínuas em I x IR n . 

303 TEU: Nessas condições, fixados Xq £ I e cq, , c n _i £ IR, existe 6 > 0 
tal que ]xq — 5, xq + <$[ Ç I e existe uma única solução y do PVI 


= r, y(x 0 ) = c 0 ,..., y (n 1} (a;o) = Cn-i 

j =0 
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definida em ]x’o — S,x o + <5[. 

Isso é um caso particular de 260. Mostraremos posteriormenteque y está 
definida sobre todo I. 

304: Quando r = 0, dizemos que a equação linear é homogênea. Caso 
contrário, a equação é não-homogênea. A parte homogênea da equação é 
obtida substituindo-se r por 0. 

Pesquisaremos esses dois tipos em sequência: 

305: Defina L n : C n {I ) -A C(I), 

n 

3 =0 

Mostre que L n “preserva combinações lineares”, ou seja, é uma transforma¬ 
ção linear. Mostre que sua imagem está contida, de fato, em C(I). 

306: O núcleo de L n é 

Nuc(L„) = {f£ (Al) | L n {f) = 0}, 

precisamente o conjunto das soluções (definidas em todo I) da parte ho¬ 
mogênea da equação linear. A superposição, nesse caso, expressa o fato de 
Nuc(L„) ser subespaço vetorial de Q n (I) (isto é, ser fechado sob combina¬ 
ções lineares). No próximo ponto, mostraremos que dimNuc(T n ) = n. 

307: A equação diferencial L n (y) = 0 tem n soluções linearmente indepen¬ 
dentes i/i,.., ,y n sobre /. Mais ainda, se y é uma solução de L n (y) = 0 
sobre I, então existem escalares Cj ,..., C n univocamente determinados tais 
que y = Ciyi + ... + 

Solução: Fixe Xo £ I- Para cada 1 j £ n. consideremos o PVI 

/•,:(?/) 0, y ik \x 0 ) = ôj( k+ 1 ) para 0 < k < n- 1, 

onde S U v = { u J v & ° usual delta de Kronecker. Pelo TEU, existe uma 
única solução )jj : I -£ IR desse PVI. 

Use o wronskiano para mostrar que {yi, ■ ■ ■ ,y n } é linearmente indepen¬ 
dente. Alternativamente, suponhamos que ciyi(x) +... + c n y n (x) = 0 para 
todo® £ I. Calculando essa combinação linear em Xq, obtemos 

Cl -jM^oT+ c 2 ... + c n -jufxóT= 0, 
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de modo que C\ = 0 e, portanto, c 2 y 2 + ... + c n y n = 0. Derivando, obtemos 
também c 2 y' 2 (x) + ... + c n y' n {x) = 0 para todo x £ I. Novamente, 


C2 ■■■ + c n 0 

e obtemos C 2 = 0. Procedendo assim, concluímos que todo Cj = 0. 

Dada a solução y: I —> IR de L n {y ) = 0, tome Cj = y <:l ^ í ' > ('Á) para 
1 ^ j ^ n e defina z = C\y\ + ... + C n y n . Por superposição, z é uma 
solução de L n {y) = 0. Pelo modo como definimos cada y :l , temos 

n n 

Z {k) (x 0) = ^C j yf ) {x 0) Xy ' : li,- y {k) ( x 0) 

i=1 J=1 

para 0 ^ k ^ n — 1. Logo, j e 2 são soluções do mesmo PVI e o TEU 
implica y = z = Y,j Cjij-j. 

Porque a base {j/i,..., y n } é linearmente independente, são únicos os 
escalares Cj para cada solução y. 

308: Em resumo, se conhecermos uma base para Nuc(L n ), a solução geral 
de L n (y) = 0 é escrita como combinação linear dessa base. 

309: Os coeficientes C \,..., C n têm CÇwnesmo papel, aqui, que a constante 
de integração nas equações de primeira ordem. 

Por exemplo, digamos que uma base de soluções seja {1, x, cos x) e que 
se queira satisfazer os valores iniciais y{ ir) = 0, y’{n) = 2, y”( tt) = — 1 . 
Então escrevemos y{x) = C\ 1 + C 2 x + C 3 cos x, donde y'(x) = C 2 — C 3 sen x 
e y"{x) = —C 3 COSX. Por substituição, temos o sistema 

1 C 1 + nC 2 - 1C 3 = 0, C 2 + 0C 3 = 2, 1C 3 = -1. 


Sua solução consiste de C\ = — 1 — 27r, C 2 = 2, C 3 = —1 e obtemos a 

solução específicá^fx) = (—1 — 27t) + 2x — cosx. 

* 

Finalmente, voltamo-nos a equações não-homogêneas: 

310: Se y v é uma solução particular (ou específica, ou sem coeficientes a 
determinar) da equação linear L n (y) = r, então a solução geral tem a forma 

V = C±yi + .. • + C n y n + y p , 

onde Ci ,..., C n são constantes a determinar e {y\, ..., y n } é um base do 
espaço de soluções de L n {y) = 0. 
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De fato, se y também é solução de L n (x) = r, então a diferença y — y p 
é solução da parte homogênea (ou seja, L n (y — y p ) = 0) e pode ser escrita 
como a combinação linear Cij/i + ... + C n y n . 

Desse modo, o conjunto das soluções de L n (y ) = r é o espaço afim 

Nuc(L„) + yp, 

com dimensão n. 

311 Resumo: Para resolver a equação linear L n (y) = r de ordem n, de¬ 
vemos resolver a parte homogênea L n (y ) = 0, cuja solução geral i/h admite 
n constantes a determinar, e encontrar alguma solução particular y p da 
equação original. A solução geral será //h + y P - 

O termo //h pode ser escrito como combinação linear (nas constantes 
a determinar) de n soluções da parte homogênea que sejam linearmente 
independentes. 

Caso se especifique também um PVI, devemos determinar os valores das 
n constantes substituindo os valores iniciais em i/ h + y p e suas derivadas, 
não apenas em j/h (a não ser que a equação original seja homogênea, quando 
obviamente podemos tomar y p = 0). 

Resolver essas equações, para obter y\, e y p , é o que faremos nas próximas 
seções. 

_Cy 

3.4 Coeficientes constantes 

Nesta seção, trataremos das equações lineares com coeficientes constan¬ 
tes, para as quais podemos encontrar todas as soluções explicitamente. 

Lembramos que devemos encontrar (i) uma base de soluções da parte 
homogênea e (ii) uma solução particular da equação original, para então 
(iii) escrever a solução geral. 

312: Recorde que as soluções da equação de primeira ordem y' + ky = 0 
são da forma y(x) = Ce~ kx . 

313: A equação de segunda ordem ay" + by' + cy = 0, onde a,b,c são 
constantes reais^ê a ^ 0, tem soluções exponenciais? 

Solução: Escrevendo y(x) = e Xx , temos y'{x) = Xe Xx e y"(x) = \ 2 e Xx . 
Substituindo na equação, obtemos a\ 2 e Xx + bXe Xx + ce Xx = 0, ou seja, 
(aA 2 + bX + c)e Xx = 0. 

Como a exponencial nunca se anula, concluímos que y(x) = e Xx é so¬ 
lução da equação se e somente se A for uma raiz do chamado polinómio 
característico P(t ) = at 2 + bt + c. 
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Note que os coeficientes de P(t) = at 2 + bt 1 + ct° são os mesmos da 
equação ay ^ + by^ + cy^ = 0. 

As raízes desse polinómio podem ser encontradas pelo método de Bas- 
khara, que nos permite saber se elas são reais ou complexas. Digamos que 
sejam Ai e A 2 : 

314: Se Ai, A 2 são reais e distintas, obtemos duas soluções da equação: e XlX 
e e^ 2 *. Elas são linearmente independentes, porque uma não é múltipla 
da outra. Portanto, {e XlX , e X2X j é uma base para as soluções da equação 
proposta. 

Atenção: Isso só vale quando podemos aplicar nosso estudo anterior. 
Aqui, precisamos ter Ai, A 2 raízes reais para que y\ , y 2 sejam funções reais. 

315: Se Ai, A 2 são raízes reais, mas idênticas, isto é. se A = Ai = A 2 é uma 
raiz dupla do polinómio característico, então, do mesmo modo, sabemos 
que e Xx é solução. 

Como A é raiz dupla do polinómio, então não somente a \ 2 + b\ + c = 0, 
mas também 2aX + b = 0, como mostra o método de Bhaskara quando o 
discriminante é nulo. 

Use isso para verificar, por substituição direta, que também xe Xx é 
solução. (Aprenderemos em 437, quandoHratarmos de redução de ordem, 
como essa solução pôde ser encontrada.) 

Novamente, as duas soluções são linearmente independentes porque não 
são múltiplas entre si. Então {e Xx ,xe Xx } é uma base para as soluções da 
equação dada. 

Para lidar com as raízèsf complexas, precisamos considerar também fun¬ 
ções e equações com valores complexos; façamo-lo brevemente neste ponto: 

316: Suponha r, s,u,*n funções reais definidas em um mesmo intervalo. Se 
r+si é contínua e se vtyvi é uma solução complexa de ay"+ by'+cy = r+si, 
então u e v são soluções reais de ay" + by' + cy = r e ay" + by' + cy = s, 
respectivamente. 

Indicação: A função complexa u + vi pode ser vista como uma função 
IR —> IR 2 , ou seja, uma curva no plano à qual se aplica o Cálculo geométrico, 
com derivação linear e feita componente a componente. 

Note: Vale também o resultado mais geral, análogo, a respeito de uma 
equação linear homogênea qualquer y pj(x)y = 0. 

Voltamos às raízes do polinómio característico: 

317: Se Ai, A 2 são raízes complexas, elas são conjugadas e podemos escrever 
àa=a + /3ieX 2 =a- fdi com a, j3 G IR e f3 ^ 0. 
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318: Sabemos, por substituição, que é uma solução complexa de 

ay" + by' + cy = 0. Porém, essa exponencial é igual a 

e ax e ipx _ e ax ^ cos (^p x ^ _|_ jsen(/3a:)] = [ e ax cos(/3a:)] + [ e ax sen(/3a:)]í. 

Então, pelo que vimos, e ax cos(/3x) e e ax sen(/3;r) são soluções reais dessa 
equação. Usamos o wronskiano para determinar que elas são linearmente 
independentes: 


e ax cos((3x) 


e ax sen(/3x) 


ae ax cos(/3a;) — /3e ax sen(/?x) ae ax sen (fíx) + (3e ax cos (f3x) 


e 2ax /3 ^ 0. 


Logo, {e ax cos(füx), e ax sen(/3a;)} é uma base para o espaço das soluções 
reais da equação. 

319 Conclusão: Para resolver a equação ay" + by' + cy — 0, primeiro 
encontramos as raízes do polinómio característico P(t) = at 2 + bt + c. A 
solução da equação é então dada em termos delas, assim: 


raízes solução geral 

reais distintas Ai, A 2 Cie XlX + C 2 e X2X 

real dupla A C\e Xx + C 2 xe Xx 

complexas a ± ôr“ Ó$, ax cos((3x) + C 2 e ax sen (f3x) 


Não demonstraremos a conclusão análoga para ordens maiores que 2, 
mas convém entendê-la bem: 

320: Analogamente, dada 

dnV^ + • • • + a 0 y = 0 com a 0 ,..., a n € IR e a n ^ 0, 

formamos o polinómio característico P(t) = a n t n + ... + ait + ao- Liste as 
raízes reais de P como Ai,..., X p , com multiplicidades k\, ..., k p , respecti¬ 
vamente, e as raízes complexas em pares conjugados ui, fj ,*,..., n q , /x*, com 
multiplicidades l \,..., l q , respectivamente. Então 

P(í)^^„(í - Ai) fcl ... (í - A P ) kp (t - m) h (t - ui) 11 ... (t- fiq) lq {t - y* q ) lq . 
Escreva Hj = a j + Pji com £ IR e (3j ^ 0. Uma base para o espaço 
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de soluções da equação dada constitui-se das funções 

■Aix „_Aix fci-l Aix 

O ^ ^ * • ■ ■ O * 


A p a? X p x rpkp- 

D * tl/C- ^ ^ 10 

-lgApX 

cos(/lix), xe aiX cos(/3ix),..., 


sen(/?ix), xe aiX sen(/3ix),..., 

x ll ~ 4 e aiX sen^Sja;), 

cos(f3 q x), xe aqX cos(/3 q x ),..., 

x lq ~ 1 e a % x Cos(/3 q x), 

sen(/3 g x), xe aqX sen(/3 q x ),..., 

x lq ~l^ ,X! sen ((3 q x). 


Indicação: É preciso verificar que (i) todas essas funções são soluções 
da equação dada, que (ii) elas são linearmente independentes entre si e que 
(iii) há precisamente n delas. 

Determine a solução geral destas equações: 

321: 3 y" - 13?/ + 12y = 0. 

Solução: O polinómio característico é P(t) = 3í 2 — 13Í+12, cujas raízes 
são 3 e 4/3, distintas e simples. Então uma base do espaço de soluções é 
dada por {e 3x , e 4 x Z 3 } e a solução geral é y(x) = C\e 3x + C 2 e 4x ! 3 . 

322: y& + 8 y'" + 16 y' Co> 

Solução: O polinómio característico é 

P(t) = t 5 ^ 8 t 3 + 16 1 = t(t 4 + 8 t 2 + 16) = t(t 2 + 4) 2 . 

Suas raízes são 0 (simples) e 2i, — 2i (duplas), de modo que 

{1, cos 2 x, x cos 2 x, sen 2 x, x sen 2 x} 

é uma base do espaço de soluções. A solução geral é 

y(x) = Ci + C 2 cos 2x + C 3 sen 2x + C^x cos 2x + C 5 x sen 2x. 

32âi.i/' 1 ' - Ay'" + 13y" = 0. 

Solução: O polinómio característico é 

P(t) = t 4 - 4f 3 + 13< 2 = t 2 (f - At + 13), 
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cujas raízes são 0 (dupla) e 2 + 3i, 2 — 3i (simples). Então 


{1, x, e 2x cos 3x, e 2x sen 3x} 

é uma base do espaço de soluções e a solução geral é 

y(x) = Ci + C 2 x + C 3 e 2x cos 3x + Cie 2x sen 3x. 

324: y" - 5 y' + 6 y = 0. 334: y'" + y" - 6 y' = 0. 

325: y" - Ay' + Ay = 0. 335: y"' + 2 y" - y' - 2y - 0. 

326: y" - 9y = 0. 336: y"' +y = 0. 

327: y" + 9y = 0. 337: y"' - 3 y" + 3 y’ -y = 0. 

328: y" - Ay' + 5y = 0. 338: y^ 9 y = 0. 

329: y" - 3 y' + 2y = 0. 339: y W +"$ y" + 9y = 0. 

330: y" + 2y' + y = 0. 340: y&*± 0. 

331: y" - 2y' + 5y = 0. 341: y ( 4 ) - 2 y" = 0. 

332: 2 y" - Ay' + 8 y = 0. 342: y^ + 2 y"' + y' = 0. 

333: y"' - y'= 0. 343: YT j=0 (”)y (i) = 0- 

Resolva os problemas de valor inicial: 

344: 3 y" - 13 y' + 12 y = 0, y{ 0);= 1, y'(0) = 2. 

Solução: Conforme 321, a solução geral é y(x) = C\e 3x + C 2 e 4x Z 3 , 
de modo que y'{x) = 3C\e 3x + | C 2 e 4x ^ 3 . Substituindo, encontramos as 
equações lCi + 1©2 = 1 e 3Ci + |©2 = 2, que são resolvidas simultane¬ 
amente por Ci = 2ffi%yC 2 = 3/5. Assim, a solução específica do PVI é 
y{x) = §e 3x + \e 4x ! 3 . 

345: y" - 6 y' + 8y = 0, y(0) = 1, y'{ 0) = -2. 

346: y" - 4^ 4y = 0, y(0) = 5, y'(Q) = 0. 

347: 0, ?/(0) = 3, y'( 0) = 0. 

348: y" — Ay = 9, y{ 0) = 0, t/(0) = 0. 

- 4 y' + 5y = 0, y( tt/2) = 0, y'( tt/ 2) = e 7r . 

350: y"' - 3y" + 3y'- y = 0, y( 1) = 2, y'(l) = -1, y"{ 1) = 0. 
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Agora, continuando a considerar equações lineares com coeficientes cons¬ 
tantes, aprenderemos a determinar uma solução particular da equação não- 
-homogênea. Lembramos a forma em consideração: 

ay" + by' + cy = r(: r), ou em geral a n y + ... + aoj/ 0 -* = r(x) 

onde os coeficientes são constantes. 

351: Pelo princípio da superposição (261), se r é uma soma algébrica de vá¬ 
rias funções rj, então podemos quebrar a questão em vários “subproblemas”, 
procurando uma solução de cada equação ay" + by + cy = e tomando a 
soma algébrica dessas soluções como solução particular de ay" + by+cy = r. 

352: Note que se r é constante e se c / 0 (ou ao^jAO) então a função 
constante y = r /c (ou y = r/ao, respectivamente) resolve a equação. 

Estudaremos dois métodos diferentes para encontrar a solução particu¬ 
lar: A variação das constantes, de Lagrange, sempre conduz a uma solução 
mas requer integrações que, freqüentemente, podem ser complicadas. Já 
utilizar os coeficientes indeterminados requer, às vezes, estimativas iniciais 
ou “chute” para ser menos trabalhoso. 

353 Lagrange: Dada ay" + by 1 -fcy = r, seja {i/ 1 , 2 / 2 } uma base do espaço 
de soluções de sua parte homogênea; então escrevemos y h = C±yi + C 2 y 2 . O 
método da variação das constantes^,'ou dos parâmetros, propõe determinar 
uma solução da equação original com a forma 

y p (x) = C\{x)yi (x) + C 2 {x)y 2 {x) 

onde substituímos as constantes C \, C 2 por funções, assim como na Se¬ 
ção 2.3. Em termos técnicos, procuramos uma solução não apenas no su- 
bespaço vetorial gerado por yi,y 2 , mas em todo o ideal que elas geram no 
anel C 2 (I). 

Lagrange obsérvou que impor a condição adicional C[yi + C 2 y 2 = 0 
simplifica enormemente os cálculos: 


e então 

.07 


V 

0 


C[yi + Ciy[ + C 2 y 2 + C 2 y ' 2 — C\y[ + C 2 y 2 
y'p = C[yi + C' 2 y 2 + C\y'[ + C 2 y 2 . 


Como desejamos ay'/ + by' p + cy p = r, substituímos: 


a(C[ 2/í + C' 2 y 2 + C\y" + C 2 y 2 ) + b(Ciy[ + C 2 y 2 ) + c(C\y\ + C 2 y 2 ) — r. 


78 


B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1- quad. 2012 


Já que Ci(ay" + by[ + cyi) = 0 para ambos i = 1,2, simplificamos para 
obter a{C {//{ + C 2 y 2 ) = r. 

Assim, 

2/1 2/2 C[ _ 0 

Vi V'i_ Ç ' 2 _r/a_ 

Esse sistema linear funcional tem determinante W = W(yi,y 2 ), que nunca 
se anula. Sua solução, portanto, é dada por 


1 

0 

2/2 

:-^ec; = 

1 

2/1 

0 

w 

r/a 

2/2 

aW 2 

" w 

2/í 

r/a 


(Ou seja, calculando as funções em cada x, temos um sistema determinado 
cujas soluções, pelo método de Cramer, também dependem de x.) Desse 
modo, obtemos C \, C 2 integrando-se tais expressões: quaisquer primitivas 
bastarão, de modo que convirá tomar as mais simples (e sem constantes de 
integração). 

354 Exemplo: y" — 3y' + 2y = sen a:. 

Já sabemos determinar bases para as soluções da parte homogênea: a 
mais simples é {e x ,e 2x }, cujo wronskiano ef 


W = = e 3x . 

e x 2e 2x 

De acordo com Lagrange, temos 

e 2x sen x __ T „ e~ x . . 

C 1 = -—^— = —e senx => Ci = — (sen a; + cosa;), 

onde integramos por partes (duas vezes), por exemplo. Analogamente, 

e^sena; _ 2r _ —e~ 2x . 

C '2 = 3j = e sen a; => C 2 = —-—(2senx + cosa;), 

de modo que uma solução particular para a equação dada é 

e~ x — e~ 2x 

—— (sen x + cos x) ■ e x H--— (2 sen x + cos x) • e 2x = A sen x + ^ cos x 

z 0 

e a solução geral é 


y(x) = D x e x + D 2 e 2x + A sen;r + éL x 
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onde os parâmetros Di,D 2 são constantes arbitrárias. 

355: Generalize o método de Lagrange para equações de ordem arbitrária. 
356 Exemplo: y'" + y' = x 2 . 

A equação homogênea associada tem base {1, cos x, sena;} com wronski- 

ano 

1 cos x sen x 
W = 0 — sena: cosa; = 1. 

0 — cos x — sen x 

Procuramos uma solução particular 

y p (x) = C\{x) + C 2 (x) cos-x f-<.):s(.rj sena;. 

Lagrange impõe 

C}(a;)l + C' 2 {x) cosa; + C 3 (x) sena; = 0, 

donde y 1 (x) = — C 2 {x) sena; + C 3 {x) cosa:, e impõe também 

C[( a:)0 + C 2 {a^{— sena;) + C 3 (x) cosx = 0, 

donde y”(x) = — C 2 {x) cosx — C 3 (x) sena; e 

y'p{x) = —C 2 (x) cosx + C 2 (x) sena; — C 3 (x) senx — C 3 (x) cosx. 

Para que y”' + y p = x 2 , portanto, devemos ter 

—C 2 (x) cosx — C 3 (x) sena; = x 2 . 

Obtemos o sistema 

1 cosa; sena: C\ 0 

0 —sena; cosa: C 2 = 0 • 

0 —cosa; —sena; C 3 x 2 
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Como W = 1, vem 


cos x sen x 


C'i(x) = 0 —sena; cosa; = x 2 => Ci(a;) = —-, 


V 


x — cos x — sen x 


C' 2 (x) = 0 0 cosa; = -x 2 cos x 


O 


0 x 2 — sen x 


1 cos x 0 


Lr) 

2 

C 2 (x) = —x sen x — 2x cos x + 2 sen x , 


C' 3 (x) = 0 -senx 0 =-x 2 senx^ 


0 — cos x x 2 


AV 


C 3 (x) = x 2 cos x — 2x sen x — 2 cos x. 


Então 


z 3 2 v V 

y P (x) = — ■ 1 + {—x 2 sen x — 2 x cos x + 2 sen x) • cos x + 

3 / 

+ (x 2 cos x — 2x sen x — 2 cos x) ■ sen x = 


= -2 x 

3 



e a solução geral é 


y{x) = Di + D 2 cos x + D 3 sen x + —— 2x 

O 

com Di, D 2 , D 3 constantes arbitrárias. 

Determine a solução geral de cada equação: 


357: y" - 3;/ • 2// = e~ 5x . 

358: y" — 6y = sena;. 

359: y" +y’ -6y = 5e 2x . 

360: y H - 5 y' + 6 y= 12. 

361: y»' + 2 y' + y = 8e x . 

362: y" — 2y’ + y + 6 sen x = 0. 


363: y" — 4 y' + Ay = x 2 . 
364: y" — y = xe x . 

365: y" + y' — 2y = 8 sen 2a:. 
366: y" + 9 y = cos 2x. 

367: y" + y = 4x sen x. 

368: y'" -8 y = e 2x . 
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369: y" — 3 y' + 2y = 70(e~ 5x + sen x). 

370: y" + y' — 6 y = 100(e 2x + sen a; + 3). 

Resolva os problemas de valor inicial: 

371: y" -4 y = x 2 e 2x , y(0) = 1, y'(0) = -1. 

372: y" — 3 y' = x + cosx, y(ir/2) = 0, y'{ tt/2) = ir. 

373: O método dos coeficientes indeterminados afirma qtef”caso r tenha 
uma das formas 

r(x) = (ckX k + ... + c\x + co)e ax cos(/3aq?Òu 
r(x) = (ckX k + ... + cia; + co)e ax sen(/3x), 

que são bastante frequentes em aplicações, então uma solução particular da 
equação linear tem a forma semelhante 

y p (x ) = x l {bkX k + ... + bix + b 0 )e ax cos(fix) + 

+ x l (dkX k + ... + d\x + d 0 )e ax sen(/3x). 

Encontram-se l, bo, ■ ■ ■, bk, do, ■ ■ ■, (4 assim: 

(i) se k = 0, mantenha bo, do, qlie se tornam constantes multiplicativas; 

(ii) l será a multiplicidade (possivelmente 0) de a + pi como raiz do po¬ 
linómio característico; 

(iii) bo,. ■., bk, do, ■ ■ ■, dk 'serão os mais simples que satisfaçam o sistema 
induzido por substituição de y p . 

Referência: Uma discussão aprofundada, com argumentação formal e 
abstrata, é feita eqi [BoDi], Seção 3.5. 

Este método, ao contrário do de Lagrange, não requer integração, mas 
aplica-se apenas a algumas formas da função r. Pode-se, caso se intua uma 
forma mais restrita para y p , eliminar diversos coeficientes imediatamente e 
aliviar o trabalho, embora isso torne o método mais artesanal. 

374 Exemplo: y" + y = l.r sen x. 

O polinómio característico é t 2 + 1. Temos 

r(x) = 4x sen x = (4a; 1 + 0)e Ox sen(lx) 
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e 0 + li tem multiplicidade 1 como raiz do polinómio. Então procuramos 

y p (x) = a; 1 (òia; + bo)e 0x cos(lx) + a: 1 (dia; + do)e 0x sen(lx) = 

= (&ia; 2 + box) cosa; + (dia; 2 + dox) sen a:. 

Com essa expressão, calculamos 

y p + y p = (2b± + 6dix + 2 do) cosx + {—Ab\x — 2b 0 + 2d\) sen a; 

(verifique!), que desejamos que valha 4a; sen a;. Para tanto, comparamos os 
coeficientes de cada termo: 


2bi T 6dia; -t- 2do — 0 
—4b\X — 2 bo + 2di = 4x 


2òi T 2do — 0 
6di = ON, 

-4òi = 4 
— 2bo T 2di = 0 


Assim, obtemos b\ = — 1, b 0 = 0, di = 0, do = Te 
y p (x) = — a; 2 cos x sen x 


e a solução geral é 

y{x) = Di cos x + D% sen x — x 2 cos x + x sen x. 

375 Exemplo: y" — y = x 2 ^. ' 

O polinómio característico é t 2 — 1. Temos 

r(x) = x 2 e x = (la; 2 + 0a; + 0)e lx cos(0a;) 

e 1 + 0í tem multiplicidade 1 como raiz do polinómio. Então procuramos 

y p (x) = x 1 (b 2 X 2 + bix + b 0 )e lx cos(0a;) = (ò 2 x 3 + 6ia; 2 + box)e x , 

onde não escrevemos um termo com sen(0a;) porque este fator o anularia. 
Substituindo-se, calculamos 

y p {x) - y p (x) = e x [6b 2 x 2 + (6ò 2 + 4 b\)x + ( 26 i + 2b 0 )\, 

que devemos igualar a a; 2 e x . Isso é feito comparando-se o coeficiente de 
cada termo: 

6ò 2 = 1, 6ò 2 + 4&! = 0, 2b x + 2 bo = 0. 
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Esse sistema é resolvido por &o = 1/4, b\ = —1/4 e b 2 = 1/6, de modo que 
y p (x) = (la ; 3 - la ; 2 + \x)e x 

e a solução geral é 

y(x ) = D\e x + D 2 e~ x + e x (2a ; 3 — 3x 2 + 3x)/l2. 

376: Resolva os exercícios anteriores utilizando, agora, o método dos coefi¬ 
cientes indeterminados. 

3.5 Gravitação universal 

Aqui, nosso objetivo é descrever a órbita de um planeta, cometa ou ou¬ 
tro corpo qualquer em torno de um sol ou grande massa central situada na 
origem de um sistema de coordenadas, órbita essa regida exclusivamente 
pela força gravitacional newtoniana. Também deduziremos a equação dife¬ 
rencial de todas as órbitas de um corpo em um campo de forças central ou 
radial, de modo geral. 

Usaremos as letras m, M para designar tanto as massas desses corpos, 
como eles próprios, e consideraremos M grande o suficiente para que a força 
exercida por m seja desprezível e não o mova da origem. 

377: Em um espaço tridimensional, o movimento de uma partícula sujeita 
a uma força central é plano, isto é, a trajetória e o centro da força estão 
contidos em um único plano. Veja [FiNe], p. 156. 

378: Desse modo, podemos limitar-nos a um plano de trabalho descrito em 
coordenadas cartesianas Lx. ou nas polares associadas via x = r cos 9, 
y = rsen$; suporemos sempre que r > 0 porque M ocupa a origem. 

O que procuramos são funções vetoriais (x(t),y(t)) ou ( r(t),9(t )) que 
descrevem a posição de m em um instante t. 

379: Derivando quanto a t, obtemos 

x — f cos 9 — rd sen 9 , 
y = r sen 9 + rd cos 9 

e então 

x = r cos 9 — 2 rd sen 9 — rd sen 9 — r{9) 2 cos 9, 
y = r sen 9 + 2 rO cos 9 + rd cos 9 — r(9) 2 sen 9, 
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ou seja, 

x cos 9 + j/sen0 = f — r(9) 2 , 
x sen 9 — y cos 9 = —2 rd — rd. 

380 2- lei de Newton: Temos sempre 

F = m £i{ x ,y) = ( mx,my ). 

aceleração 

381: Quando o campo é centrífugo, isto é, a força atua no mesmo sentido 
do raio vetor, vale F = |F|(cos 9, sen 9); quando o campo é centrípeto, vale 
F = — |F|(cos0,sen0). 

Então x = ±^77 cos 6 * e y = ±^7 sen 9 em que o sinal é positivo quando 
o campo é centrífugo e negativo no caso centrípeto: substituindo, obtemos 

f — r{9 ) 2 = ±|F|/mv 
2 rè + rd = 0 . 

Desse sistema, derivamos a conservação do momento angular e a equação 
de Binet, que serão as únicas a serem utilizadas no restante do raciocínio. 

382 Conservação angular: O momento angular de m é 

L = mrV = mr 2 9 

(onde V é sua velocidade linear). Vemos que 

L = ^(?nr 2 0) = 2 rrtrfÒ + mr 2 0 = mr(2r9 + rd) = 0, 

de modo que L é constante. 

383: Há, portanto, duas possibilidades: 

Se L = 0, então 9 = 0 e m move-se em linha reta radial, afastando-se 
ou caindo em M. 

Se L ^ 0, então ou 9 > 0 ou 9 < 0 (já que mr 2 > 0), fazendo com que 
9 seja sempre estritamente crescente ou decrescente. Assim, um corpo em 
órbita não inverte o sentido de sua translação. 

384 2 - lei de Kepler: A área varrida pelo raio vetor de m entre os instantes 
ío ©*Çé 
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onde integramos em coordenadas polares. Então 


À(t) = 9(t) = L/2m 

e concluímos pela constância da velocidade areolar. Isso implica que m 
tem velocidade menor no ponto da órbita mais afastado de M e velocidade 
maior no ponto mais próximo. 

385 Binet: Note que 


l ) = -r- 2 % = -r- 2 %/% = -r~ 2 r/è = -r 


logo, 

&ÁV- 1 ) = M-™/L) = -Wi)(f/ f )<5 - mf/Lè . 

Substituindo do sistema, vem 


^(r- 1 ) = -m[r(9 ) 2 ± \F\/m] A/=y l/r =F \F\/LÒ. 


Então, para r em função de 9, temos 


d 2 ( 1\ ^ m|Ejr 2 

dO 2 Vr/ r / L 2 ’ 

em que o sinal é negativo no caso de campos centrífugos. 

Até aqui, trabalhamos com qualquer campo central: a equação de Binet 
e a 2- lei de Kepler são gerais. 

Agora, veremos como a lei da gravitação universal de Newton mantém 
os corpos em órbitas, que são secções cônicas, descritas pelas três leis de 
Kepler. 

386: No caso particular da força gravitacional F exercida por M sobre m 
em ( x,y ), seu valor absoluto é |E| = GMm/r 2 e o campo é centrípeto. 
Assim, 


387 1- lei de Kepler: Com z = l/r e |F| = GMm/r 2 , obtemos da 
equação de Binet esta equação linear de 2- ordem: 

d 2 z GMm 2 

dÕi +Z= L 2 ' 

A parte homogênea tem solução geral z\, = C\ cos 6 + C 2 sen 9, que pode 
ser reescrita como z h = —k cos (9 — 9 0 ) com k > 0; a escolha desse sinal será 


* 

rs/ 


EV= — |F|(cos 0 ,sen 0 ) = — 


GMm 
( x 2 + y 2 ) 3 / 2 


{x,y). 
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justificada a seguir. Por outro lado, uma solução particular da equação 
original é a constante z p = GMm 2 / L 2 . Logo, 

—= —kcos(9 — 9 0 ) + GMm 2 / L 2 , 


ou seja, 


„ (g) = E/k 

1 ; 1 — E cos(9 — 9 0 ) 


se definirmos E = kL 2 /GMm 2 > 0. 

Escreva uma solução análoga para um campo centrífugo. 

388 Cônicas: Mostre que isso é a equação de uma cônica de excentricidade 
E ^ 0 com um foco em M e sendo 9q o ângulo entre seu eixo principal e 
o das abscissas. (A cônica é uma parábola se E = J.,Mipse se E < 1, ou 
hipérbole se E > 1.) Conclua que planetas não caem em espiral. 

Indicação: Consideremos o caso da elipse, quando E < 1. Uma apli¬ 
cação bruta das transformações de coordenadas, ou uma espiada apressada 
em um formulário, pode induzir-nos a trabalhar com a origem no centro 
da elipse (a tradicional fórmula x 2 /a 2 + y 2 /b 2 = 1 ), mas em nosso caso 
a origem está no foco onde se encontra M. Outra dificuldade é que uma 
mudança de coordenadas envolve um cálculo que, começado do jeito errado, 
pode ser muito trabalhoso. 

Suponhamos, porém, que a elipse tenha semi-eixo maior a e semi-eixo 
menor b. Então sua excentricidade é E = V a 2 — b 2 /a, a distância entre 
seus focos é 2 \Ja 2 — b 2 , e a soma das distâncias de um ponto na curva aos 
focos é 2a. Uma delas é r, de modo que a outra (no lado oposto ao foco em 
M) é 2a — r. O ângulo entre o raio vetor de medida r e o eixo principal, 
dependendo das posições relativas, é 9 — 9q ou 9q — 9 ou 2n ~ 9 + 9o, todos 
com mesmo cosseno. Do triângulo formado pelos dois focos e o ponto da 
curva, obtemos 

(2a — r ) 2 = ^ 2 \Ja 2 — b + (r ) 2 — 2 ^ 2 \Ja 2 — b 2 ^ (r) cos (9 — 9 0 ). 
Simplifique e resolva em r para determinar 


1 — E cos (9 — 9o) 1 — E cos(0 — 9 0 ) ’ 

em que fazemos k = E/[a( 1 — E 2 )\. 

389: Quais são as possíveis trajetórias de um elétron defletido por uma 
fonte negativa? 
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390 Periodicidade: Quando E < 1 e a órbita é elíptica, mostre que r é 
periódico em função de t. (Abaixo, chamaremos T a esse período.) 

Solução: Quando a órbita é uma elipse (E < 1), a função r(9) está 
definida para todos os ângulos (seu denominador nunca se anula) e é ex¬ 
plicitamente periódica quanto a 9, com período 2tt. Porque a velocidade 
areolar é constante, o tempo necessário para o raio vetor varrer toda a elipse 
é sempre o mesmo independentemente do ângulo inicial. Isso corresponde 
a uma translação completa, ao longo da qual 6 aumenta ou diminui 2 ir e, 
assim, r é periódico também quanto a t. 

391 3- lei de Kepler: Mostremos que a razão T 2 /g? é a mesma para 
qualquer corpo em torno de M, sendo T o período dç translação e a o 
raio médio de uma órbita elíptica (ou seja, seu semi-eixo maior). (Sejam 
também, como acima, b seu semi-eixo menor e E sua excentricidade.) 

Vimos que a velocidade areolar é constante L/2m, de modo que nab/T = 
L/2m. Então T = 2'Kabm/L e 


T 2 47r 2 ?n 2 6 2 

a 3 L 2 a 


47r 2 m 2 

~ 1 Ã~ 


a(l - E 2 ) = 


^{E/k) = ^ 


L 2 

GAÍm 2 


47T 2 

GM 


é constante e independe de m è a. 

Reciprocamente, a lei da gravitação de Newton é uma conseqüência da 
astronomia de Kepler; tal fato tem grande valor científico e histórico, porque 
se deduz uma lei geral e teórica a partir de leis experimentais específicas: 

392 Dedução da gravitação universal: Postula-se que a força gravita- 
cional é atratora, ao longo da reta entre os dois corpos, ou seja, o campo 
gravitacional é (central) centrípeto. Calculamos A = L/2m e, assim, a 
2- lei de Kepler afirma que L é constante. 

A 1- lei informa que elipses 


r( 6 ) 


E/k 

1 — Ecos(9 — 9 0 ) 


onde E < 1 


satisfazem a equação de Binet para campos centrípetos. Porém, 


+ \ = (k/E)(Eco S (9 - 9o)) + (k/E)( 1 
donde k/E = m\F\r 2 /L 2 . 
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Enquanto mostrávamos a 3 â lei de Kepler, calculamos 


47T 2 77p 


«(1 - E 2 ) = 


47 r 2 m 2 


(E/k) 


exclusivamente a partir da velocidade areolar. Agora, assumindo essa lei, 
escrevamos T 2 /a 3 = Cm constante: obtemos 


C M = 


Alt 2 Til 

\F\r 2 ' 


Pelo princípio da ação e reação de Newton, m também atrai AI gravita- 
cionalmente com uma força de mesmo módulo |.Fj. Por simetria, então, 

47 t 2 m 47 t 2 Aí ’ Ay 

C M r 2 C m r 2 

em que C m é outra constante. Note que cada constante depende apenas 
do corpo correspondente. Portanto, tomando m = I, deduzimos que Cm = 
Ci/M; definindo G = 47 t 2 /Ci e retomando m arbitrário, obtemos 

47t 2 m ÇAÍm 
1 1 = Ci/M ' Ã 2 = r 2 ’ 

Encerramos nossos cálculos com duas extensões da teoria desenvolvida: 

393 Movimento de Aí: Supusemos sempre que M não se movia e per¬ 
manecia como centro de nosso sistema de coordenadas. Contudo, pode-se 
requerer o estudo de órbitas quando m é comparável com AI: por exemplo, 
nos casos Terra-Lua, Sol-Júpiter, estrelas binárias, etc. O que se pode 
supor, nesses casos, é que o centro de massa de AI e m tem movimento 
inercial, isto é, retilíneo e uniforme. 

Alterações: Qualquer referencial inercial, agora, é igualmente ade¬ 
quado. Deve-se descrevér o movimento de m, M por vetores V m , Vm, respec¬ 
tivamente; as velocidades e acelerações pertinentes são V rn . Vm e V m , Vm ■ 
Desse modo, pode-se escrever a segunda lei de Newton para m e Aí, com as 
forças -GAIm(V M -V m )/\V M -V m \ 3 e -GAÍm(V m -V M )/\V m - V M \ 3 , res¬ 
pectivamente. Subtraia essas equações e conclua que as duas primeiras leis 
de Kepler não sofrem modificação, enquanto que a terceira transforma-se 
em T 2 //tP = 47 T 2 /G(AI + m). 

394 Relatividade: Considerando-se efeitos relativísticos, a equação de 
Binet transforma-se na equação de Schwarzschild que não é linear: sua re¬ 
solução complica-se, mas pode-se fazer uma boa aproximação para explicar 


B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1° quad. 2012 



a discrepância da órbita de Mercúrio. Veja [FiNe], p. 174-176, e [0’Ne], 
p. 372-380. Nessa teoria, também, corpos que transladam dispendem ener¬ 
gia mecânica na forma de ondas gravitacionais e, portanto, planetas caem 
em espiral, embora seu tempo de queda seja muito superior à duração do 
Universo. 

3.6 Osciladores harmônicos 

Osciladores harmônicos correspondem a classes especiais de equações 
lineares de 2- ordem com coeficientes constantes, que descrevem diversos 
fenômenos físicos de interesse. 

Primeiramente, apresentaremos a teoria geral dos osciladores harmôni¬ 
cos, para depois contextualizar alguns exemplos. 

395: O oscilador harmônico simples é a família de fenômenos descritos pela 
equação 

ax + cx = 0, 

sendo as constantes a,c, em geral, estritamente positivas. 

396: Também se estuda o oscilador amortecido com equação 

ax + bx + cx = 0, 

onde ainda b é uma constante positiva. 

Em estudos aplicados, o termo bx pode ter outra forma, mais complexa, 
dependendo da modelagem feita. Veremos, porém, que a linearização (que 
consiste em substituir o termo real pelo linear) simplifica enormemente 
sistemas complexos, a um certo preço. 

397: O oscilador é forçado se houver parte não-homogênea em sua equação, 
geralmente uma expressão senoidal: 

ax + bx + cx = r(x). 

Vamos estudar todas essas classes, uma por vez: 

398: A equação ax + cx = 0 com a, c > 0 tem polinómio característico 
P(u ) = au 2 + c com raízes complexas ázi^/c/a, de modo que sua solução 
geral é 

x(t) = Ci cos + C 2 sen = 

= Acos(cot + 9), 
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onde uj = \fcfa e A, 9 são constantes a determinar, com A > 0. Chamam-se 
A de amplitude, 9 de fase e ca de frequência angular do oscilador. 

399: Note que x é periódica com período T determinado assim: 

oo(t + T) + 9 = (cot + 9) + 2tt <t=> T = 2n/io, 

ou seja, lo = 2ir/T = 2irf onde / = T -1 é a frequência do oscilador 
harmônico. 

400: Por simples derivação, obtêm-se estas expressões para a ‘Velocidade” 
e a “aceleração” do oscilador: 

x(t) = —lo Asen(wí + 9) = Alo cos(tot + 9 + 
x(t) = —oo 2 Acos(iot + 9) = Alo 2 cos (cot 7r). 

Desse modo, x tem fase atrasada 7r/2 em relação a x e amplitude Alo, 
enquanto x tem fase atrasada ir em relação & x e amplitude Alo 2 . 

401: Vemos também que x = —lo 2 x. 

402: Os valores iniciais tradicionais são A0i = Xq e i:(0) = vq. Então 

xo = A cos(w0 + 9) = A cos 9, 

v 0 = Aio cos(w0 4 - 9 + |) = — Aco sen 9, 

o que permite determinar A e 9, ou mostrar diretamente que 
x(t) = aio cos (cot) + — sen(wí) 

LU 

(caso lo = 0 então x = aio). 

403: A energia do oscilador harmônico é definida assim: 

E = ^a(x) 2 + \cx 2 . 

Os termos E(f*= a{x) 2 / 2 e Ep = cx 2 /2 são chamados energias cinética e 
potencial, respectivamente. 

No caso de alguns fenômenos mecânicos, como 417, essa definição mo 
tiva-se somando-se a energia cinética usual e uma energia potencial igual, 
em valor absoluto, ao trabalho efetuado sobre a massa a da origem 0 até o 
valor atual de x: 

nX nX pX 

7 = / “força” ds= (as) ds= (—cs) ds = —\cx 2 

J 0 J 0 J 0 
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já que as + cs = 0. Em geral, seja qual for o fenômeno estudado, podemos 
definir E e suas componentes cinética e potencial como acima. 

404: Temos 

È = ■ 2xx + ■ 2xx = x(ax + cx) = 0, 

de modo que E é constante. Calculando-a em t = —9/lo, por exemplo, 
obtemos 

E = ^a(x) 2 + \cx 2 = |a(0) 2 + |c(A) 2 = 

ou seja, a energia é proporcional ao quadrado da amplitude. 

405: Resolver a equação ax+bx+cx = 0, com b > 0, requer conhecer o sinal 
de b 2 — 4ac para que se determinem as raízes de au 2 + bu + c. Aqui, 

suporemos que b é muito pequeno, correspondente a um amortecimento 
ligeiro. Desse modo, ò 2 — 4ac < 0 e as raízes de P(u) são complexas com 
parte real —b/2a. A solução geral então é 

x(t) = Ae~ bt / 2a cos (cot + d), onde lo = X = 

2a 

ou seja, agora lo depende também de'T>, mas ainda é constante. 

406: Faça o gráfico de x(t) para t & 0, encapsulando-o entre os gráficos de 
Ae~ bt / 2a e — Ae~ bt l 2a . Note, agora, que A é a amplitude máxima válida 
apenas para t = 0 (realizada se 9 = 0). 

407: Se b 2 — 4 ac ^ 0, não há oscilação e o sistema apenas converge ao 
equilíbrio (a; = 0). 

408: Novamente se define a energia 

'•(7 

E = \a{x) 2 + \cx 2 . 

Oy 

Agora temos 

E = x{ax + cx) = x(—bx) = —b{x) 2 ^ 0, 
de modo que E é decrescente. 

409: Também para o oscilador amortecido convém considerar a quantidade 

^nat = \/ c/«. 
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chamada frequência angular natural do oscilador e obtida desprezando-se o 
amortecimento. 

410: A oscilação é forçada quando uma força excitadora senoidal com pe¬ 
ríodo próprio é aplicada ao oscilador: 

ax + bx + cx — Fq cos(u)qÍ), 

com Fq,loo constantes positivas. (Isso assume que a força excitadora tem 
valor máximo justamente em t = 0 ; em uma situação prática, basta trans¬ 
ladar horizontalmente a solução que obteremos.) 

411: Caso se deseje estudar uma força excitadora constante, basta pôr 
loq = 0, de modo que a solução particular é x = Fq/c constante. 

Doravante, assumiremos coo > 0: 

412: Defina os seguintes valores: 


Aq = 



se lo o 7 ^ 
se cüq = 



Verifique, então, que 

xjf íi) =/ Aq cos (w 0 í + 0 o ) 


Se LOq ülnat, 

se loq = w nat ; 
se loq < Lo na t* 


é uma solução particular dessa equação, conhecida como solução de estado 
permanente. 

413: Já a solução Xh(£) = Ae~ ht l 2a cos (u>t + 9) da parte homogênea é cha¬ 
mada transiente, porque seu valor é desprezível após um intervalo de tempo. 
Assim, a solução completa 


x(t) = x h (t) + x p (t) « x p (t) 

é inicialmente complicada, mas seu comportamento converge ao de x p (í); 
em particular, sua frequência converge à mesma da força excitadora. 

414: Quais são as alterações caso 5 = 0, ou seja, o oscilador é forçado, mas 
não amortecido? Trabalhe com a hipótese loq ^ w n at- 
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415 Ressonância: Quando wo = w na t, isto é, a frequência da força excita- 
dora é igual à frequência natural do oscilador, diz-se que ocorre ressonância 
e temos, por substituição e simplificação, 

x p {t) = p*- cos(w 0 f - §). 
ou o 

Nesse caso, a amplitude A 0 = F 0 /bui 0 é a máxima possível. A fase do 
oscilador está atrasada n/2 em relação a da força excitadora Fq cos(woí); 
porém, 

x P (t) = sen(w 0 í - f) = ~r cos(w 0 í), 

de modo que a velocidade do oscilador está em fase com a força excitadora, 
ou seja, o movimento nunca se opõe à força, possibilitando a absorção 
máxima de energia. 

416: Finalmente, consideramos a ressonância sem amortecimento, ou seja, 
ò = 0ewo=w = w nat : mostre que uma solução particular é 

. . FqLü . 

x p{t) = sen(wt), 

cujo comportamento leva a amplitudes cada vez maiores. Mostre que a 
energia desse oscilador é, a partir de certo valor para t, sempre crescente e 
ilimitada. 

Um mecanismo ou estrutura que esteja submetido a um tal oscilador não 
aguentará estresses crescentes e colapsará. (O mesmo efeito pode ser veri¬ 
ficado acima, com b muito pequeno e a solução evolui para uma amplitude 
A 0 muito alta.) 

Já resolvemos a equação de Binet (385) por uma mudança de variáveis 
que a torna um oscilador harmônico simples. Os seguintes exemplos apresen¬ 
tam diversas outras manifestações dos osciladores harmônicos. Complete 
os detalhes da formulação, se necessário: 

417: Suponha que um bloco de massa m desloca-se sobre um trilho hori¬ 
zontal, atado a uma mola com constante elástica k. Fixe como origem a 
posição do bloco quando a mola tem seu comprimento natural e assuma 
que o bloco apenas se desloca em um intervalo pequeno, de modo a sempre 
valer a lèj. de Hooke e nunca torcer ou inverter a mola, ou colidir com o 
ponto fixo da outra extremidade. 

Se o bloco mover-se exclusivamente sob ação da força elástica, temos, 
pela 2- lei de Newton e porque a força elástica opõe-se ao movimento, 

mx = — kx => mx + kx = 0, 
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que é um oscilador harmônico simples. A discussão sobre energia, acima, 
vale literalmente. 

Se o movimento do bloco sofrer atrito do trilho, supõe-se que esse atrito 
seja proporcional à velocidade do bloco, mas oposto: então 

mx = —kx — bx => mx + bx + kx = 0, 

que é um oscilador harmônico amortecido. 

Descreva um sistema massa- mola, como este, mas em que haja força 
excitadora. 

418: Suponha, agora, que o bloco esteja suspenso pela mola e que se mova 
verticalmente. Assuma que a mola não tenha massa. Sendo g a aceleração 
gravitacional e supondo que o eixo esteja orientado para cima, temos 

mx = —kx — mg. 

Determine a mudança de variável que transforma essa equação em uma de 
oscilador harmônico simples, sem excitação ou amortecimento. 

419 MCU: Suponha que um pino é fixado à borda de um toca-discos ou 
roda de ceramista e iluminado por uma fonte de luz de raios paralelos. 
Mostre que o movimento de sua sombra reta em uma parede é um oscilador 
harmônico simples. 

420: Na situação descrita no ponto 31, assuma que a sonda permanece 
sempre no interior do túnql, óiUseja, no interior do planeta. Identifique o 
oscilador harmônico associado. 

Na Seção 2.6, estudamos os circuitos elétricos RC e LR. Podemos com¬ 
pletar esse estudo, agora, com os circuitos LC e LCR, que são regidos por 
osciladores harmônicos. Relembre o funcionamento de capacitores e indu¬ 
tores, respectivamente, em 226 e 231. 

421 LC: Ligamos, em série, uma chave interruptora, um indutor de indu- 
tância L e um capacitor de capacitância C, ambas L , C constantes positivas. 
Então, sendo Q a) carga armazenada no capacitor e I a corrente no circuito, 
temos I = Q e LI + Q/C = 0, isto é, a diferença de potencial no indutor 
neutraliza aquela no capacitor, donde 

LQ + Q/C = 0. 

O PVI usual é dado por Q(0) = Qq, carga armazenada inicialmente no 
capacitor, e 1(0) = 0, quando a chave é fechada no instante t = 0. 
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Assim, Q é um oscilador harmônico com frequência angular w = 1/ y/LC 
e I é sua “velocidade”. Como I é senoidal, trata-se de um circuito de corrente 
alternada: o indutor “facilita tanto” a descarga do capacitor que termina 
por transferir a carga inteira para a outra folha do capacitor, espelhando 
a situação inicial; o capacitor carrega-se e descarrega-se, sempre oscilando 
entre +<5 0 e — Q 0 . 

As energias elétrica e magnética do circuito interpretam as energias 
potencial e cinética que discutimos acima. 

422 LCR: Este circuito conecta, também em série, um resistor de resistên¬ 
cia R ao circuito LC acima, sendo ainda R constante positiva. Então as 
quedas de potencial no resistor, no capacitor e no indutor são, respectiva¬ 
mente RI, Q/C e —LI, donde 

LQ + RQ + — = 0 
e Q é um oscilador harmônico amortecido com 



e I é sua “velocidade”. Trata-se de um circuito de corrente enfraquecendo-se, 
alternada para valores usuais de R e “contínua” quando R é muito grande. 

423: Agora, suponha que se insira nesse circuito, ainda em série, uma fonte 
de força eletromotriz alternada, isto é, que origina a diferença de potencial 
E(t) = E 0 cos(woí). Então RI + Q/C + LÍ = E, donde 

LQ + RQ + ^ = E 0 cos (w 0 t) 

O 

e Q é um oscilador harmônico amortecido e forçado e, mais uma vez, I é 
sua “velocidade”. (Lembramos que uma força eletromotriz contínua pode 
ser estudada com ujq = 0.) 

Finalmente, estudaremos o oscilador harmônico associado ao movimento 
pendular. Trabalharemos sempre com um pêndulo de haste rígida com com¬ 
primento R e toda a sua massa m concentrada pontualmente na extremi¬ 
dade livre; pêndulos reais são conhecidos em [Tip]. 

Deduziremos a equação do movimento de um pêndulo pela análise das 
energias envolvidas no processo, sendo desnecessário o estudo de configura¬ 
ções de forças. 
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424 Pêndulo simples: Aqui, supomos que a haste revolve sem qualquer 
atrito. Postulamos que o pêndulo move-se em um plano e que pode ficar em 
repouso em duas posições: a “estável” (em termos práticos) é a do pêndulo 
pendurado, e a “instável”, em que a haste sustenta a massa, é impossível no 
caso de um pêndulo de fio esticado. 

Note: No vocabulário das Seções 2.1 e ??, esses pontos não são singu¬ 
laridades atratoras ou repulsoras, mas do tipo elíptico: a massa do pêndulo 
oscila em torno da posição “estável”, ou no complemento de um entorno da 
“instável”. O pêndulo não só pode oscilar comumente, como também partir 
de perto da posição instável e quase dar a volta, em grande amplitude, ou 
girar sem inverter o sentido de rotação, com velocidade máxima no ponto 
do equilíbrio estável e mínima no outro. 

425: Seja a o ângulo da haste com a vertical (0 na posição estável), consi- 
derando-se também alguma orientação para a positivo ou negativo. Então 
a velocidade angular do pêndulo é à e a velocidade linear da massa m é 
V = àR. Consideremos as energias cinética Eç = mV 2 /2 e potencial Ep da 
massa m: se E é a energia mecânica total do sistema, portanto constante, 
temos E = Eç + Ep. (Estas não são as energias associadas ao oscilador 
harmônico como definimos acima.) 

426: Adotamos a convenção de que Ep = 0 para o ponto mais baixo da 
trajetória do pêndulo. Verifique, assim, que Ep = mgR( 1 — cosa). As pos¬ 
sibilidades que descrevemos ocorrem conforme E = 0 (equilíbrio estável), 
E = 2mgR (equilíbrio instável), 0 < E < 2 mgR (oscilações) ou E > 2 mgR 
(giros completos). 

427: Então 

\m{a) 2 R 2 + mgR( 1 — cosa) = E. 

Derivando quanto ao tempo e simplificando, obtemos 

Rá + g sen a = 0. 

Essa não é uma equação de oscilador harmônico, porque não é linear. 

428: Para pequenas amplitudes, sen a ss a, de modo que à equação do 
pêndulo simples corresponde o oscilador harmônico, ou equação linearizada, 

Rà + ga = 0 

que descreve o fenômeno físico aproximadamente. Sua solução é 

a(t) = A cos (cot + 9) 
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com ui = \Jg/R, ou seja, o período é T = 2i:\JRjg, como aprendemos em 
Física. 

Note: Nessa solução, a é um ângulo descrito em termos trigonométri¬ 
cos, cujos argumentos é que esperaríamos que fossem ângulos. Se quiser 
considerar, portanto, a posição vertical ou horizontal da extremidade com 
massa m, deveremos manipular 

cosa = cos(Acos(uit + 9)). 

Por mais estranho que seja, não há nada errado nisso. 

Embora a equação linear seja muito mais simples de resolver e estudar, 
enquanto sua solução seja uma aproximação muito boa ao que efetivamente 
ocorre, é importante lembrar que o sistema original não é linear e tem um 
comportamento muito diferente. Isso é relevante não apenas em Física e nas 
mais diversas engenharias, mas também em diversos estudos econômicos, 
demográficos e situacionais, onde a tomada de decisão governamental ou 
corporativa deve estar atenta a fatores de caos. 

429: Podemos também considerar a equação original, sem linearização: 



onde K = (2 E/mR 2 ) — (2 g/R) é constante. Essa é uma equação de 1- or¬ 
dem; com qual método podemo^resolvê-la? 

430 Pêndulo amortecido: No pêndulo simples, a energia total se man¬ 
tinha constante. Passemos a “supor que o pêndulo sofre algum atrito ou 
amortecimento e que isso implica em uma perda de energia, até o instante 
t, de Ea(í) ^ 0 (sendo £a(0) = 0). Temos, portanto, Eq + Ep + Ea = Eq, 
onde Eq é constante igual à energia inicial do sistema. 

431: Suponha)“por exemplo, que a força resistiva seja proporcional à velo¬ 
cidade da extremidade do pêndulo: \F\ = kV com k > 0. Então a energia 
perdida é igual ao trabalho que essa força realiza do instante inicial até o 
instante “considerado, digamos t\. 

£píi)= f \F\ds= ( 1 \F\-f t dt = f 1 \F\Vdt= í 1 kV 2 dt, 

■J trajetória J 0 J 0 J 0 

donde 

È A = kV 2 = k{à) 2 R 2 . 
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Como anteriormente, derivamos Eq + Ep + Ea = Eq e simplificamos para 
obter 

kR . 

Ra - a + g sen a = U, 

m 

cuja linearização é o oscilador amortecido 

kR . 

Ra -I- a + qa — ü. 

m 

3.7 Coeficientes variáveis 

Agora, trataremos das equações lineares em geral, cujos' coeficientes 
podem ser funções da variável independente: 

Pn(x)y (n) + ... +pi{x)y' +p 0 {x)y = r(¥), 

em particular as equações de segunda ordem 

p{x)y" + q(x)y' + u(x)y = r(x). 

Novamente, devemos encontrar (i) soluções básicas da parte homogênea e 
(ii) uma solução específica da equação dada, para depois (iii) obter a solução 
geral. 

Comecemos por determinar a solução particular: 

432: O método da variação das constantes de Lagrange (descrito em 353 e 
seguintes) aplica-se também às equações de coeficientes variáveis, bastando 
novamente que se conheça uma base de soluções da parte homogênea. A 
demonstração é idêntica à do caso de coeficientes constantes. 

433: Especificamente para uma equação de segunda ordem com a base 
{yi, 2 / 2 } e W = IF(yuJfoVtemos 

V 2 _ _ V2lL e c 1 - — Vl ° - — 

r /p y'i pW 2 w y' x r /p pW ’ 

onde agora também p depende de x. 

434: Também o método dos coeficientes indeterminados (373) aplica-se ao 
novo caso, embora os coeficientes da equação devam também ser considera¬ 
dos na proposta de uma forma para a solução particular. 

435 Exemplo: ( x 2 + x)y" — (x 2 — 2 )y' — (x + 2)y — x(x + l) 2 , conhecendo 
a base {e x ,x~ 1 } para a parte homogênea. 
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Deve-se verificar que realmente 2 / 1 ( 2 ) = e x e 2 / 2 ( 2 ) = x _1 são soluções de 
(x 2 +x)y" — (x 2 — 2)y' — (x + 2)y = 0, o que podemos fazer por substituição. 
Também se verifica, usando-se por exemplo o wronskiano, que 2 / 1 5 2/2 são 
linearmente independentes. De fato, temos 

W(x) = W(y u y 2 )(x) = 6 X = -e^aT 1 + x~ 2 ) ^ 0. 

e —x 

Cabe uma nota: p(x) = x 2 +x anula-se em x = — 1 e x = 0, de modo que 
devemos estudar a equação separadamente em ]— 00 , — 1, 0[ e ]0, 00 [. 

(Isso posto, poderíamos dividir toda a equação por p( .>")■) Observe que, 
nesses intervalos, 2 / 1 , 2/2 estão definidas, de classe e W realmente não 
se anula. 

Então, pelo método de Lagrange de variação das constantes, procuramos 
uma solução particular da equação original tendo a forma y p (x) = Cj (x)e x + 
C 2 (x)x~ 1 , para o que devemos ter 


C'(x) = X 1 -< x + l ) 2 = Xe -x 

lK} {x 2 +x)-l~e*(x~ 1 +x- 2 )} ’ 

donde C\(x) = — (x + l)e~ x (integrando-se por partes), e 

C> (x] = , ^x(x + l) 2 _ 2 

’ (x 2 Ct-®)-[-e*(*- 1 +a:- 2 )] 

donde C 2 (x) = —x 3 /3. (Aqui, ignoramos as constantes de integração.) 
Então a solução geral s é* 

y(x) = D\e x + D 2 x~ l + [— (x + l)e~ x ]e x + [—x 3 /3]x _1 = 

= D\e x + D 2 x~ x — (1 + x + x 2 /3), 


onde agora D\,D 2 são parâmetros arbitrários, mas constantes. 

436: Para usar o método dos coeficientes indeterminados, digamos que 
já suspeitássemos de um polinómio de grau 2 como solução da equação: 
escrevemos 

2 /p (x) = ax 2 + bx + c 

onde a,b,c são constantes a determinar. Como então 2 /p( 2 ) = 2 ax + b e 
2 /p ix) = 2a, substituindo na equação devemos ter 

(x 2 + x)(2 a) — (x 2 — 2) (2 ax + b) — (x + 2 )(ax 2 + bx + c) = x(x + l) 2 , 
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ou seja. 


(—3a)x 3 + (—2 b)x 2 + (6a — 2b — c)x + (2b — 2c) = x 3 + 2x 2 + x. 

Comparando-se esses polinómios, temos —3a = 1, —2 b = 2, 6a — 2b — c = 1 
e 2b— 2c = 0, de modo que a = —1/3, b = —1 e c = —1, tal qual deduzimos 
acima. 


Para determinar uma base de soluções da parte homogênea da equação 
linear de coeficientes variáveis, o estudo das raízes do polinômib caracte¬ 
rístico não é válido. E comum que a equação homogênea tenha, de fato, 
soluções exponenciais e trigonométricas, mas precisamos de outros cami¬ 
nhos para encontrá-las exatamente. Nossa abordagem será via redução de 
ordem: 

437 Redução da 2- ordem: Suponha já conhecida uma solução 2/1 de 
p(x)y" + 9 ( 2 ) 2 / + it(:r)y - 0: 

devemos determinar outra solução y 2 , linearmente independente de 2 / 1 , para 
formar a base { 2 / 1 , ' 2/2 } • 

Vimos em 294 que o wronskiano W = W(yi,y 2 ) deve satisfazer a equa¬ 
ção pw' + qw = 0 e, então, ser dado pela solução de Abel-Liouville (299). 
Explicitamente, temos 


2 / 1 ( 2 ) 2 / 2 ( 2 ) 

2 /í ( 2 ) 2 / 2 ( 2 ) 


= ^(^=W(x 0 )exp(^- j 


Expandindo o determinante e substituindo expressões em x para j/i, 2/í 
e a integral, obtemos uma equação de primeira ordem cuja incógnita é y 2 . 
Basta, portanto, resolvê-la segundo os métodos do Capítulo 2. 

Para que y 2 seja independente de 2 / 1 , devemos ter W(yi,y 2 ) ^ 0. Para 
isso, impomos tV(xid) = 1 ou outra constante não-nula que simplifique 
os cálculos, assim como um valor simples para Xq, desde que no domínio 
correspondente. 

438 Exemplo: (x + l)y" — (x + 2)y' + y = 0, conhecendo 1 / 1 ( 2 ) = e x . 
Substituindo, verificamos que 2/1 é realmente solução dessa equação. 
Note que p(x) — x + 1 anula-se em x = —1, de modo que devemos 
estudar a equação separadamente em ]— 00 , —1[ e ] —l,oo[. Trabalhemos 
primeiro com ] — 1, 00 [: 
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Fixamos Xq = 0 G ] — 1, oo[ e W(xq) = 1, de modo a anular a primitiva 
abaixo e facilitar os cálculos. A formulação acima torna-se 



Portanto, 

e x y' 2 - e x y 2 = exp([s + ln \s + l|]g) = exp(x + ln \x + 1|), 

ou seja, 

y' 2 -y 2 = e ln ^ =x+í 

já que trabalhamos com x > — 1. Essa é uma equação de primeira ordem: 
como aprendemos na Seção 2.3, sua solução é 

2 / 2 ( 2 ;) = -(x + 2) + Ce x . 

Podemos tomar C — 0, porque precisamos de uma única expressão para , 
e obtemos a solução geral 

y(x) = C 1 ef + C 2 (x + 2) 


da equação dada, onde absorvemd&.o sinal em C 2 . 

Para resolver no domínio ]— 00 * — rfá note (pelo Teorema do Valor Inter¬ 
mediário) que existe Xq < —1 tal que Xq + ln|xo + 1 | = 0 ; não é preciso 
explicitá-lo para conduzir cálculos como os acima e obter a equação 

2/2 — 2/2 = —x — 1, com solução 2 / 2 ( 2 ;) = x + 2 + Ce x . 

Novamente, a solução geral é y(x) = C\e x + 62 ( 2 : + 2). 

439 Exemplo: y" + 4 xy' + (4x 2 + 2 )y = 0, com y\{x) = exp(— x 2 ). 
Verifique que 2/1 é solução da equação dada. Desta vez, temos 


jm_ 2 . _ 2 _cy 2 

ou seja, jjT - ® y 2 + 2xe~ x 2/2 = e _ x , donde 

2/2 + 2 £ C £/2 = e - ® 2 . 

Uma solução dessa equação, também pela Seção 2.3, é 2 / 2 ( 2 ;) = xe - * , de 
modo que a solução geral da equação dada é 

y(x) = Cie~ x + C 2 xe~ x . 
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Dada uma solução das partes homogêneas, determine a solução geral 
destas equações: 

440: xy" — (x + 2 )y' + 2y = 0, y\ = e x . 

441: xy" + (x — 1 )y' — y = 0, t/i = x — 1. 

442: xy" + 3y' = 0, 2/i = 1. 

443: (x - l)j/" - xy' + y = 0, 2/1 = x. 

444: x 2 í/" + xy' — Ay = 0, y\ = x 2 . 

445: xy" — y' + 4x 3 y = 0, yi = sen(x 2 ). 

446: (x 2 + x)y” — (x 2 — 2)y' — (x + 2)y = x(x + l)^^y= e*. 

447: x 2 y" + xy' — 4y = 10x 3 , yi = x 2 . 

448: x 2 y" — 2y = —2x, 2/1 = x _1 . 

449: x 2 y" — xy' + y = x, yi = x. 

450: Quando a, b, c são constantes reais, sendo a / 0, e A é raiz dupla do 
polinómio P(t) = at 2 + bt + c, aprendemos que e Xx é solução da equação 
ay" + by' + cy = 0 e verificamos queSêfií' também é solução. Como podería¬ 
mos, inicialmente, encontrar esta solução, conhecendo apenas e Xx ? Mostre 
os detalhes do cálculo. 

Quanto a reduzir uma ordem arbitrária ou maior que 2, temos as se¬ 
guintes possibilidades: 

451: Conhecidas n —X, soluções linearmente independentes de uma equação 
linear homogênea de ordem n, o cálculo do wronskiano e o uso da solução 
de Abel-Liouville como acima induzem uma equação linear não-homogênea 
de ordem n — 1 para a n-ésima solução. 

452: Suponhamos conhecida uma solução yi da equação linear homogênea. 
Então as substituições y = y\Z e z = Z’ reduzem a ordem em uma unidade, 
preservando linearidade e homogeneidade. 

Solução: Por indução, mostra-se que 

y w = (yiZ) il) = É (]) y{ j) z {i ~ j) 

1=0 
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para cada 0 ^ i ^ n. Substituindo todos y(°\ ..., na equação original, 
obtemos uma equação da forma 

q„(x)ZW + .. . + qi(x)Z' + q 0 (x)Z = 0. 

Já que y = yi é solução da equação original, então Z = 1 é solução desta 
equação, de modo que qo(x) = 0. Logo, pondo z = Z r , obtemos 

q n {x)z { - n ~ 1) + ... + qi(x)z = 0 , 

como desejado. 

453: Vejamos isso explicitamente para o próprio caso n = 2: com y = y\Z, 
temos 

y' — y'\Z + y\Z' e y" = y”Z + 2 y' % Z' + y\Z". 
Substituindo-se em py" + qy' + uy = 0, vem 

p{y"Z + 2y[Z' + y\Z") + q{y[Z + y\Z') + u(yiZ) = 0, 

ou seja, 

\pyi}Z" + [2 py[ + qyi]Z' + [py” + qy[ + uyx]Z = 0, 

onde o terceiro coeficiente é nulo porque y\ é solução da equação original. 
Assim, com z = Z ', 

\pyi\z' + i 2 py'i + qyi\z = o. 

que é uma equação linear homogênea de ordem 1 . 

454: Resolva os exercícios anteriores utilizando, agora, este método para 
redução de ordem. 

Resta o problema de determinar y -[, isto é, uma primeira solução explí¬ 
cita (e não-nula) da parte homogênea da equação linear: 

455: Em alguns casos, isso pode ser feito com o próprio método dos coe¬ 
ficientes indeterminados, aplicado à parte homogênea. (O método de La- 
grange produz apenas a solução trivial!) Também se podem tentar outros 
métodos, como a determinação de série de potências, que estudaremos opor¬ 
tunamente. 
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3.8 Equações e métodos especiais 

Esta seção considera equações de diversas formas. Em particular, para 
algumas equações lineares homogêneas de coeficientes variáveis, apresenta 
ao menos uma solução, cumprindo com uma pendência da seção anterior. 
Ainda, apresentamos outros métodos de resolução que também podem ser 
aplicados a equações de primeira ordem. 

456 Euler: Sejam a,b,c constantes quaisquer, para as quais consideramos 
a equação linear homogênea 

ax 2 y" + bxy' + cy = 0. 

Forme o polinómio Q(t) = at 2 + (b — a)t + c (que não se chama “polinómio 
característico”) e mostre que: 

(a) x x é solução dessa equação se e somente se A é raiz de Q\ 

(b) se (b — a) 2 = 4ac e A é a raiz dupla de Q, então x x ln \x\ também é 
solução da equação; 

(c) se as raízes de Q são complexas a±BL com a, /? reais, então X a cos(/3 ln x) 
e x a sen(/j lna;) são soluções dessa equação quando x > 0 . 

457 Exemplo: x 2 y” — 6 xy' + 12 y = 0. 

Formamos o polinómio lí 2 + (—6 — 1 Jt + 12 = t 2 — 7t + 12, cujas raízes 
são 3 e 4. Assim, as funções linearmente independentes a : 3 e x 4 formam 
uma base para as soluções dessa equação. 

Determinar séries de potências por substituição na equação dada é um 
procedimento similar ao método dos coeficientes indeterminados. Eis o 
exemplo clássico: 

458 Potências: y" = — y, y{ 0) = 1, y'(0) = 0. 

Solução: Escrevemos uma série de potências com centro no valor inicial 
dado: neste caso, u,= a n{ x ~ 0) n - Derivamos a série termo a termo, 

formalmente e sem compromisso com convergência, como se fosse um po¬ 
linómio gigante: y' = na n x n ~ l e então y" = E^ 2 n ( n — 1 )a n x n ~ 2 . 

Substituímos ,as ! séries na equação e reindexamos as séries de modo que 
todas tenham as potências na mesma forma; como nossa equação envolve 
apenas y" e y , re-escrevemos y = Y ^=-2 a n- 2 X n ~ 2 . (Caso houvessem co¬ 
eficientes funções na equação envolvendo x, também desenvolveríamos em 
séries e multiplicaríamos dentro das somas.) Temos 

r oqj oo oo 

n(n — 1 )a n x n ~ 2 = — a n - 2 X n ~ 2 =>• ^^(n(n — 1 )a n + a„_ 2 )x n ~ 2 = 0 . 

n=fy n —2 n =2 
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Já que a série obtida é nula, anulamos cada coeficiente seu, formalmente 
e sem considerar convergência: n(n — l)a„ + a n - 2 = 0 para cada n ^ 2. 
Vemos que a n = —a n _ 2 /(n(n — 1)), de modo que se n é par então a n = 
(—l)”/ 2 ao/n! e se n é ímpar então a n = (—l)(” _1 )/ 2 ai/n!. Assim, 

OO OO OO 

y = J2 a nX n = ^(-l) fc a 0 cr 2fe /(2fc)! + ^(-l) fe ai x 2fc+1 /(2fc + 1)T = 

n—0 /c—0 fc—0 

= ao cos x + ai sen x, 

em que rearranjamos a série sem preocupação com convergência. Obtivemos 
uma família de funções candidatas a solução da equação dada; verificamos 
que é o que de fato ocorre, e procedemos normalmente para determinar 
soluções restantes; também podemos estudar o raio de convergência da 
série obtida na solução. Quanto ao PVI, substituímos os valores iniciais na 
série obtida e suas derivadas, obtendo 1 = aoO° + aiO e 0 = aoO + aiO°; 
logo, ao = 1 e ai = 0, com y = cosa:. 

Frequentemente, uma série de potências tem um número finito de coefi¬ 
cientes não-nulos, ou seja, é simplesmente um polinómio. Aqui estão alguns 
exemplos relacionados a equações diferenciais: 

459 Legendre: Verifique que a equação 

(1 — a: 2 )|/”V 2s wij' + n(n + 1 )y = 0, 

onde n € IN, possui uma solução polinomial de grau n chamada polinómio 
de Legendre. 

460 Hermite: Mostre que, se n £ IN, a equação 

y" — 2xy' + 2 ny = 0 

possui uma solução polinomial de grau n, o polinómio de Hermite. 

Há dois casos em que se pode reduzir a ordem de uma equação diferencial 
com a forma propícia. Lembre que a forma geral da equação de segunda 
ordem é F(x,y,y',y") = 0, para a qual se procuram soluções y(x). E com 
essa& letras que rotulamos os casos abaixo: 

461 Ausência de y: Suponha que y não conste na equação: diz-se que “a 
variável dependente está ausente”. Substitua z = y' e z* = y", obtendo a 
equação F(x,z,z') = 0, que é de primeira ordem. A solução procurada y 
será uma primitiva de z. 


106 

B. EDO ©2012 Vinícius Cifú Lopes. Versão preliminar: UFABC, 1- quad. 2012 


462 Ausência de a:: Este é o caso de equações “autônomas” que corres¬ 
pondem a “fenômenos independentes de calendário”. Se “a variável indepen¬ 
dente (, x ) estiver ausente”, fazemos z = y' e notamos que 

„ , dz dz dy dz 

^ dx dy dx dy 


Assim, a nova equação é F(y, z, z^) = 0, onde a variável independente 
agora é y, de primeira ordem. 

Atenção: Não se trata de achar 2 como função de x. Peljt) método 
proposto, após determinarmos z(y), devemos resolver a equação yf = z(y), 



e integrar ambos os lados quanto a x. (Enfim, vemos que há duas equações 
de primeira ordem a resolver.) 
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y = C/tgx. 186: Fator a: -4 ; sol.: y = ±\/ x 2 — Cx 3 . 187: Fator t~ 3 ; sol.: 
x = ±\/T — 2Uf 2 . 188: Fator a; -4 / 3 . 189: Fator (3a; + l) -2 / 3 . 190: Fa¬ 
tor (x + l) -2 ; sol.: y = (x + l)(ln|x + 1| — C) — 1. 193: Dada g, tome 
f(s) = g(l, s) e use a = a; -1 ; dada /, tome g(x,y) = f(y/x). 194: Derive 
y = xz e iguale a f(z). 197: y = ±xy/2 ln x + C. 198: Analog. a 196, 
y = C — 2 \JCx para 0 < x < C e y = —C + 2 \J—Cx se — C < x < 0. 
203: y = 3x/5 — \fí tg(— x\/3Ã) + C)/5\/ò. 204: Multiplique a equação 
por y~ k e calcule z' em separado. 206: y = (C/x 2 ~ 2x/3 )~ 1 ^ 2 . 207: y = 
[C exp(x(l - 7r)) - 2a; - 2/(1 - tt)] 1 /^" 1 ). 213: T = -a(T - A). 216: Li¬ 
near: por variação da constante, T(t) = e~ at Tn + e~ at fg ae as A(s) ds. 
217: T(t) = A + e~ at (To — A). 220: A energia térmica do sistema 

rri\C\Ti + m 2 c 2 T 2 é constante. 223: Um oitavo. 225: Temos I(x + 
h ) = I(x) — al(x)h. Divida por h e tome h —> 0: V = —ai. Solução: 
I(x) = Ioe~ ax . Assumindo /(l) = Iq/2, vem a =^2 e 1(3) = I 0 / 8 . 
226: Q(t) = Q 0 exp(-t/RC) e I(t) = I 0 exp(-t/RC). 227: Q(t) = 
EC( 1 - exp (-t/RC)) e I(t) = E exp(-t/RC)/R. 228: Q(i) = EC- 
(EC-Q 0 )exp(-t/RC) e I(t) = (EC-Q 0 )exp(~t/RC)/RC. 231: I(t) = 
E( 1 — exp (—Rt/L))/R. 232: /fi na i = E/R. 262: O PVI de segunda 
ordem especifica os valores da solução e de sua derivada em um ponto. 
Se duas soluções têm o mesmo valor em um ponto, deverão ter deriva¬ 
das distintas aí. 268: L.D.: 1(0) + 0(5a;) + 0(—2a; 3 ) = 0. 269: L.D.: 
3(2sen(a; 2 - 1)) + 2(-3sen(a; 2 - TK= 0. 270: L.I. 271: L.D.: x = 
|(x + 5) + §(x —4). 272: L.I. 273: L.I. 274: L.I. 275: L.I. 276: L.I. 
277: L.I. 290: Suponha L.D.: uma das funções escreve-se como combina¬ 
ção linear das demais. Então também uma coluna da matriz é combinação 
linear das demais e o determinante anula-se. 324: y = Cie 2x + C 2 e 3x . 
325: y = C x e 2x + C 2 xe 2x : 326: y = Cie 3x + C 2 e~ 3x . 327: y = 

C\ cos 3a; + C 2 sen 3a;. 328: y = C\e 2x cosx + C 2 e 2x senx. 329: y = C\e x + 
C 2 e 2x . 330: y = + C 2 xe~ x . 331: y = C\e x cos 2a; + c 2 e x sen 2a;. 

332: y = C\e x cos \/&ç + C 2 e x sen \/3x. 333: y = C\ + C 2 e x + C 3 e~ x . 
334: y = C 1 + Çg.e% + C 3 e~ 3x . 335: Note P( 1) = 0e fatore. Sol.: 
y = C\e x + C 2 c~ r + C 3 e~ 2x . 336: Note —1 raiz simples: fatore. Sol. 
y = C\e~ x + C 2 e x / 2 cos(v / 3a;/2) + C 3 e x ! 2 sen(v / 3a;/2). 337: y = C\e x + 
C 2 xe x + C 3 x&y 338: y = C \+ C 4 xe~' /3x . 
339: y = Ciço¥\/3x + C 2 sen\/3x + C 3 xcos\/3x + C 4 Xsen\CSx. 340: y = 
Ci + C 2 x b C 3 x 2 + C 4 X 3 . 341: y = C\ + C 2 x + C 3 e^ x + C^-^. 

342: y = Ci + C 2 cosa; + U 3 sena; + Ura; cos a: + Usa; sen a;. 343: Mos¬ 
tre F*($M= (t + í) n . Sol. y = C j+ ix j e~ x . 345: y = 3e 2x - 2e 4x . 

346: y = 5e 2x - \Ç>xe 2x . 347: y = 3 cos 2a;. 348: y = 0. 349: y = 
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— e 2 x cosx. 350: y — e x 1 + 5xe x 1 — 4x 2 e x 4 . 355: Dica: forme 


Vi(x) ••• y n (x) 



' 0 ' 

y[(x) ••• y'„(x) 


'C[ 


0 

y[ n ~ 2) (x) ■■■ y ( n~ 2 \x) 


C' n _ 


0 

.yi l l \ x ) ■■■ y { n~ 1 ] (x)_ 


0 ^, 

r/a 


357: y = Die x + D 2 e 2x + e~ 5x /42. 358: y = Die 2x + D^ 3x - (7 sen a; + 
cosa;)/50. 359: y = Die 2x + D 2 e~ 3x +xe 2x . (A sol. particular pode conter 
—e 2cc /5; absorva-o em Di.) 360: y = Die 2x + D 2 e 3x + 2. 361: y = 
Die~ x + D 2 xe~ x + 2e x . 362: y = Die x + D 2 xe x — 3 cos a;. 363: y = 
Die 2x + D 2 xe 2x + (2x 2 +4x + 3)/8. 364: y = Die x +D 2 e~ x + (x 2 ~x)e x /4. 
365: y = Die x + D 2 e~ 2x — | (3 sen 2a; + cbs/ír)*, 366: y = Z?i cos 3a; + 
D 2 sen3a: + |cos2a;. 367: y = ZAcosa; + D 2 senx + xsenx — a: 2 cosa;. 
368: y = Die 2x + D 2 e~ x cos \/3x + D.f** 1 T ^1 y/3x + xe 2 x /12. 369: y = 
Die x + D 2 e 2x + 7 sen a; + 21 cosa; + he~ bx /3. 370: y = Die 2x + D 2 e~ 3x + 
20xe 2x — 14 sen a; — 2 cos a; — 50. (A sol. particular pode conter —4e 2x \ 
absorva-o em D\.) 371: y — ^e 2x + jM e ~ 2x + e 2 x (x 3 /12 — a; 2 /16 + a;/32). 
372: y = (tt 2 /24 - tt/ 3 + 8/27) + (7?r/18 + l/270)e 3a: - 3 ’ r / 2 - ^(cosa; + 
3sena;) — x 2 /6 — x/9. 418: Dica: Translade x para que a origem seja onde 
o bloco permanece em equilíbrio, isto é, a mola é esticada apenas pelo peso 
do bloco. 440: y = Cie x ^?C 2 (x 2 +2x+2). 441: y = Ci(x—1)+C 2 e~ x . (Use 
a; 0 = 1; para integrar xe~ x r/>(x— l) 2 , faça e~ x /(x— \)+e~ x /(x —\) 2 e integre 
o 2° termo por partes, cancelando o 1°.) 446: Confira 435. Para encontrar 
y 2 = x~ 4 , a redução de ordem requer: para integrar (x 2 — 2 )/(x 2 + x), 
escreva-a 1 — (2/x) + l/(x+l); exponencie a primitiva obtendo e x (x+l)x~ 2 ; 
na eq. 1- ordem, integre e~ x (x + l)x~ 2 fazendo f e~ x xdx por partes que 
já indica toda a primitiva. 447: y = Uicc 2 + C 2 x ~ 2 + 2a; 3 . 448: y — 
Cix 2 + U 2 a; _1 -pfe. 449: y = Cix + U 2 a;lna; + (lna;) 2 /2. 


> 
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